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Abstract 

We study the limit of the hyperbolic-parabolic approximation 

C vt + A(v e , ev e x )v% = sB(v s )v s X!C V s e R N 

{ v E (0, x) = v . 

The function fi is defined in such a way to guarantee that the initial boundary value problem is well 
posed even if B is not invertible. The data g and Vo are constant. 
When B is invertible, the previous problem takes the simpler form 

f v! + A{v e , ev%)v e x = eB(v s )v s xx v s e 
I v e (t,0) = v b 
[ v e (0, x) = vo- 

Again, the data Vb and vo are constant. The conservative case is included in the previous formulations. 

It is assumed convergence of the v e , smallness of the total variation and other technical hypotheses 
and it is provided a complete characterization of the limit. 

The most interesting points are the following two. 

First, the boundary characteristic case is considered, i.e. one eigenvalue of A can be 0. 

Second, as pointed out before we take into account the possibility that B is not invertible. To deal 
with this case, we take as hypotheses conditions that were introduced by Kawashima and Shizuta relying 
on physically meaningful examples. We also introduce a new condition of block linear degeneracy. We 
prove that, if it is not satisfied, then pathological behaviours may occur. 
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1 Introduction 

The aim of this work is to describe the limit of the parabolic approximation 

vl + A(v e , sv%)v% = eB{v £ )v e xx v e e «. N 

Q(v s (t,0))=g (1.1) 
v £ (0, x) = v 

for e — * + . In the previous expression, the function fi is needed to have well posedness in the case the 
matrix B is not invertible. The function fi is defined in Section [21 When the matrix B is indeed invertible, 
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system takes the simpler form 

v e t + A(v e , ev%)v% = eB{v e )v e xx v e G R N 

v"(t,0) = v b (1.2) 

V S (0, x) EE Vq. 

Even if the equations in and (|1.2[) are not necessarily in conservation form, nevertheless the conservative 
case 

is included in the previous formulation. Indeed, one can define 

A(v e , S v s x ):=Df(v s )-e(B(v s )) x 

and obtain an equation like the one in (|1.1|> or (|1.2jl . 

In the present paper we assume that when £ — » + the solutions w e converge, in the sense that will be 
specified in the following, to a unique limit v. Since in both and (|1.2|) the initial and boundary data are 
constant, then the limit u solves a so called boundary Riemann problem, i.e. an hyperbolic initial boundary 
value problem with constant data. Results in 2 show that the study of boundary Riemann problems is a 
key point to determine the semigroup of solutions for an hyperbolic initial boundary value problem. We will 
came back to this point at the end of the introduction. 

The goal of the work is to determine the value of v{t, x) for a.e. point (t, x). In particular, we determine 
the value of the trace v of the limit on the axis x = 0. The reason why this is interesting is the following. 
Let us focus for simplicity on the case in which B is invertible, i.e. on (|1.2[) . It is known that, in general, 
the trace of the limit, which we denote by v, is different from Vb, the boundary datum imposed in (|1.2p . The 
relation between Vb and v is relevant for the study of hyperbolic initial boundary value problems and was first 
investigated (as far as we know) in [27] . Also, in [28] it was proved that the value of v in general depends on 
the choice of the matrx B. In other words, if in ljl.2|) one keeps fixed vq, Vf, and the function A and changes 
only the matrix B, then in general the value of v will change, even if the system is in conservation form. 

The most interesting points in the work are the following two. First, we cover the characteristic case, 
which occurs when an eigenvalue of the matrix A can attain the value 0. The non characteristic case occurs 
when none of the eigenvalues of A can attain the value 0. The characteristic case is more complicated to 
handle than the non characteristic one. Loosely speaking, the reason is the following. Suppose that the fc-th 
eigenvalue can assume the value 0. Then in the non linear case we do not know a priori if the waves of the 
fc-th family are entering or are leaving the domain. 

The second point is that we cover the case of a non invertible viscosity matrix B. To tackle this case we 
assume that the so called Kawashima Shizuta condition holds. We also introduce a new condition of block 
linear degeneracy. We provide a counterexample which shows that, if this condition is violated, then there 
may be pathological behaviours, in the sense that will be specified in the following. 

The exposition is organized as follows. In Section ll.il we give an overview of the paper, introducing the 
main ideas involved in the analysis. In Section [2] we discuss the hypotheses assumed in the work. In Section 
Owe give a characterization of the limit of the parabolic approximation ([1.20 . i.e. when the viscosity matrix 
B is invertible. Finally, in Section [4] we discuss the limit of (jl.lj) when the matrix B is singular. 

1.1 Overview of the paper 
1.1.1 Section [2} hypotheses 

Section [5] describes the hypotheses exploited in the work and it is divided into three parts. 

Section |2~T1 describes the hypotheses assumed in the case the matrix B is invertible. These hypotheses 
were already considered in several previous works and they are automatically satisfied when the system 
admits a dissipative entropy. 

In Section |2~21 we discuss the hypotheses assumed in the case the matrix B is singular. These hypotheses 
can be divided into two groups. 
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The first group is composed by conditions that were already exploited in several previous works (e.g. in 
plH [351 [Ml H3 El H3 EH)- In particular one assumes that there exists a regular and invertible change of 
variables u — u(v € ) such that the following holds. If v satisfies 



vl+A[f,ev%)vl=eB{v B )vl 



XX 1 




(1.3) 




(1.4) 



for a suitable b(u) € M rxr . Also, B and A satisfy suitable hypotheses that are reasonable from the physical 
point of view since they were introduced in 34, 35, 36, 37J relying on examples with a physical meaning. In 
particular, we assume that the so called Kawashima Shizuta condition is satisfied. 

Apart from these hypotheses, we introduce a new condition of block linear degeneracy, which is the 
following. Let 



be the block decomposition of A corresponding to Q1.4p . namely An and En belong to M( Ar_r ) x ( JV ~ r ) and 
A22 and £22 belong to M rxr . The condition of block linear degeneracy says that, for every given real 
number <r, the dimension of the kernel of [.An (it) — aEn(u)] is constant with respect to u. In other words, 
the dimension of the kernel may vary as a varies, but it cannot change when u varies. 

Block linear degeneracy is not just a technical condition. Indeed, in Section 12.2.21 we discuss counterx- 
amples which show how, when the block linear degeneracy is violated, one can have pathological behaviors. 
More precisely, we exhibit examples in which block linear degeneracy does not hold and there is a solution 
of Q which is not C 1 . These can be considered a pathological behaviour since one usually expects that 
the parabolic approximation has a regularizing effect. 

On the other side, block linear degeneracy is not an optimal condition, in the following sense. It is possible 
to show that block linear degeneracy is satisfied by the Navier Stokes equation written using Lagrangian 
coordinates, but it is not satisfied by the Navier Stokes equation written in Eulerian coordinates. On the 
other side the two formulations of the Navier Stokes equation are equivalent, provided that the density of 
the fluid is strictly positive. This remark was first proposed by Frederic Rousset in [UJ and it suggests that 
it is interesting to look for a condition strong enought to prevent pathological behaviours but at the same 
time sufficiently weak to be satisfied by the Navier Stokes equation written in Eulerian coordinates. This 
problem is tackled in the forthcoming paper [TB] , 

In Section [2.2. II it is defined the function fi which is used to define the boundary condition in JO}. The 
point here is the following. 

Consider an hyperbolic initial boundary value problem, and for simplicity let us focus on the conservative 
case: 



It is known that if one assigns a datum u € (i.e. if one assigns N boundary conditions), then the initial 
boundary value problem (|1.6p may be ill posed, in the following sense. In general, there is no function u 
which is a solution of 

Ut + f(u) x = 

in the sense of distributions, which assumes the initial datum for t — and satisfies 




(1.5) 




(1.6) 



lim u(t, x) 



= u 
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for almost every t. Also, it is know that a necessary condition to obtain a well posed problem is to assign 
a number of conditions on the boundary which is in general smaller then N. Assume the boundary is non 
characteristic, i.e. assume that none of the eigenvalues of the jacobian Df(u) can attain the value 0. In 
this case, one can impose on the boundary datum a number of conditions equal to the number of positive 
eigenvalues of Df{u). 

We can now came back to the parabolic equation 

u t + A{u, u x )u x = B(u)u xx . (1.7) 

Let us write u = (uj., u 2 ) T , where u\ £ M. N ~ r , u 2 6 M r . Here r is the rank of B, as in (|1.4|) . With this 
notations equation (|1.7p can be rewritten as 

f u u + A n u lx + A 2 iu 2x = , , 

\ u 2 t + A 12 u lx + A 22 u 2x = bu 2xx 

Roughly speaking, the reason why one has to introduce the function 13 is the following. Let tin be the 
number of strictly negative eigenvalues of the block An, which itself is a (N — r) x (N — r) matrix. Also, 
let q denote the dimension of the kernel of An. The second line in (II. 8|) contains a second order derivative 
of u 2 and hence u 2 can be seen as a parabolic component. On the other side, only first derivatives of u\ 
appear and hence u\ can be seen as an hyperbolic component. Actually, there is an interaction between the 
two components (this is ensured by the Kawashima-Shizuta condition). On the other side, because of the 
hyperbolic component one is not completely free to assign the boundary condition in (jl.ip . As pointed out 
in previous works (e.g in |43j ) the number of conditions one can impose on the boundary is N — nu — q. 
Indeed, one can impose r conditions on u 2 . On the other hand, one can impose on u\ a number of conditions 
equal to the number of positive eigenvalues of An, i.e. to N — r — rin — q. Summing up one obtains exactly 
N-nn-q. 

Thus, the function 6 in (jl.ip takes values in M. N ~ n ^~^ and g is a fixed vector in M. N ~ nil ~ q . The precise 
definition of 6 is given in section 12.2.11 and it is such that the initial boundary value problem (jl.ip is well 
posed. In Section 14.31 it is given a more general definition for the function fi. 

Section l2 . 2 . 21 discusses three examples. The first two show that, if the condition of block linear degeneracy 
is violated, then there may be solution of 

u t + A(u, u x )u x = B(u)u xx 

exhibiting pathological behaviors, in the sense explained before. More precisely, the first example deal with 
steady solutions 

A(u, u x )u x = B(u)u xx , (1.9) 
while the second one deals with travelling waves, 

[A(u, u) — aE{u)]u = Bu" . 

In the previous expression, a represents the speed of the wave and it is a real parameter. Finally, the third 
example in Section f2 . 2 . 21 shows that if the rank of the matrix B is not constant, then there may be solutions 
of (|1.9p exhibiting pathological behaviours of the same kind discussed before. 

Section 12.31 discusses the hypotheses that are assumed in both cases, when the matrix B in (jl.ip is 
invertible and when it is not. It is assumed that the system is strictly hyperbolic (see Section 12.31 for a 
definition of strict hypcrbolicty) . Also, it is assumed that when e — > + the solutions of v e of (jl.ip converge 
to a unique limit. Also, it is assumed that the approximation is stable with respect to the initial and the 
boundary data and that the limit has finite propagation speed. 

We refer to Section [2~3l for the exact statement of the hypotheses, here instead we underline another point. 
The proof of the convergence of v e in the case of a generic matrix B is still an open problem. However, 
there are results that provide a justification of our hypotheses. In particular, in [28] it is proved the local in 
time convergence in the case B is invertible, but in general different from the identity. Moreover, in [S] the 
authors proved the global in time convergence in the case of an artificial viscosity (B(v e ) is identically equal 
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to In)- The analysis in [5J exploits techniques that were introduced in [5J [51 QZ3 [II] to deal with the Cauchy 
problem. In [3J it is proved the same kind of convergence we assume in the present properties. Also, other 
properties we assume here (stability of the approximation, finite propagation speed of the limit) are as well 
proved in [3J. Analogous results were proved in [35] for a special class of problems with 2 boundaries. 

Also, we point out that there are several works that study the stability of the approximation in the case 
of a very general viscosity matrix B. Actually, the literature concerning this topic is very wide and hence 
we will quote only works that concern specifically initial boundary value problems: [371 1321 133"! 144) . 

1.1.2 Section [3J the characterization of the hyperbolic limit in the case of an invertible 
viscosity matrix 

Section [3J discusses the characterization of the limit of (jl.2j) when the matrix B is invertible. Actually, 
because of the hypotheses we assume in Section 12. 1[ we study the equivalent (in the sense specified therein) 
problem 

E{u e )u e t + A(u £ , eu%)u% = eB{u e )u% x u e G R N 

u £ (t,0) = u b (1.10) 
u £ (0, x) = uo. 

Also, Section [3J is divided into fours parts. 

Section 13.11 collects preliminary results that are needed in the following. 

Section f3 . 2 1 gi ves a quick review of some results concerning the characterization of the limit in the Riemann 
problem. These results were introduced in [7J. 

The Riemann problem is a Cauchy problem with a picccwisc constant inial datum with a single jump. 
Let us focus for simplicity on the conservative case: 

x =0 

(1.11) 

A solution of p. lip was first described in [35] assuming some technical hypotheses (i.e. that all the fields 
are either genuinely non linear or linearly degenerate). Since we will need in the following, here we briefly 
review the ideas exploited in [3J5] to obtain a solution of Ql.lljl . 
Denote by 

Ai(u) < • • ■ < Xn(u) 

the eigenvalues of the jacobian Df(u) and by r\(u) . . . rjy(u) the corresponding eigenvectors. For simplicity, 
we assume that all the fields are genuinely non linear. In this case, one can chose the orientation of r$ in such 
a way that VAi • r.i > 0. For every i = 1 ... N we denote by S z s .(u + ) a curve in M. N which is parameterized 
by Si and which enjoys the following property. For every value <S*. (u + ) there exists a speed A close to Aj(u + ) 
such that the Rankine Hugoniot condition is satisfied, i.e. 




/(«+)- /(£*,(«+)) =x[u + -S i Si (u+)]. (1.12) 



Also, let lZ l s .(u + ) be the integral curve of r^(tt) starting at u + , in other words TZ\ . (u + ) is the solution of the 
Cauchy problem 

4-K i = n(tf l (u + )) 



ds 

TV {u+) = u+ 



Finally, let T*.(u + ) be defined as follows: 



ls > [u > ~ \ Sl(u+) Si < 
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In it is proved that is a C 2 curve such that 



dT l St (u+) 



(is; 



r l (u + ). 



In the following, we will say that Tj (w + ) is the i-th curve of admissible states. Indeed, every state T*. 
can be connected to u + by either a rarefaction wave or a shock which is admissible in the sense of Liu. In 
other words, if Si > then the solution of the Riemann problem 



is 



u(t, x) 



Ut + f(u) x = 

«(0, ar) = | ^ (u+) a;<0 

,T > t\i(u + ) 



(1.13) 



s<tAi r. 



,(« + )) 



(1.14) 



tAi(l* («+))< a: <tAi(u + ), a: = tA; (l*(u+)) 



The meaning of the third line is that u(t, x) is equal to T*(u + ) (the value assumed by the curve at the point 
s) when x is exactly equal to Aj evaluated at the point T*(u + ). The value of u is well defined because of the 
condition VA^ • fj > 0. 

On the other side, if Sj < then the solution of the Riemann problem (|1.13|) is 



u(t, x) 




(1.15) 



The speed A satisfies the Rankine Hugoniot condition and it is close to A,(u + ). 

In this way one obtains N curve of admissible states (u + ), . . . T^ N (u + ). To define the solution of the 
Riemann problem (|l.lip one can proceed as follows. Consider the function 

V(si ...s N , u+) = T s \ o . . . T»Z\ ° T s N N (u+). (1.16) 

With the notation T^ n z] a T^ N {u + ) we mean that the starting point of T^ n z] is T^ N (u + ), i.e. that 



T, 



N-l 



N , 



It is proved in [39] that the map ip is locally invertible with respect to s% . . . s^. In other words, the values 
of si ... sn are uniquely determined if one imposes 

u~ = ip(si . . . s N , u + ), 

at least if u~ is a close enough to u + . One takes u~ as in (jl.lip and obtains the values s\ . . . s^. Indeed, 
we assume, here and in the following, \u + — u~\ << 1. Once s\ . . . sn are known, one can obtain the limit 
gluing together pieces like (|1.14[) and (|1.15[) . 

The construction in (33] was extended in [ID] to more general systems. Also, in [7] it was given a 
characterization of the limit of 



f + A(u s , eu%)u% = eB{u E )u E x 
u+ x>0 
u~ x < 



u 

u £ (0, x) 



(1.17) 



when \u + — u is sufficiently small (under hypotheses slightly different from the ones we consider here). 
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The construction works as follows. Consider a travelling profile for 

u t + A(u, u x )u x = B(u)u xx , 

i.e. u(x — at) such that 

B(u)u" = (A(u, u') - crE(u)ju'. 
In the previous, expression, the speed of the wave a is a real parameter. Then u solves 

u' — v 

B(u)v' = (A(u, v) - (xE{u)\v (1.18) 
a 1 = 0. 

The point (it + , 0, Ai(it + )) is an equilibrium for (|1.18p . Also, one can prove that a center manifold around 
0, Xi(u + )) has dimension N + 2. 
We recall here that every center manifold is invariant with respect to (I1.18[) and moreover satisfies the 
following property: let (u°, p°, cr°) belong to M c and denote by (u(x), p(x), <r(x) the orbit starting at 
(u°, p°, cr°). Then 

lim (u(x), p(x), a(x))e- cx / 2 = (0, 0, 0) lim f u(x), p(x), tr(xj) e- cx ' 2 = (0, 0, 0). 

x — >+oc V / x — y— oo V / 

The constant c in the previous expression is strictly positive, depends on the matrix E~ x A and it is defined 
in Section [2] (Hypothesis [3|) . 

Fix a center manifold M. c . If (u, v, a) is a solution to Ijl.l8|> laying on Ai c , then 

U 1 = Vifi{u, V i} (Ti) 

v[ = (f>i{u, Vi, (Ji)vi (1-19) 
a> = 0. 

The functions f{ and (pi are defined in Section 13.2.11 

The construction of T l s .u + works as follows. Fix Si > and consider the following fixed point problem, 
defined on a interval [0, Sj]: 

' u(t) = u+ + [ T h(u(0, «i(0, *i(0K 

^( T ) = /i( T 7 w , Vi, (Ti) - conc/i(r, u, Vi, (Tj) (1-20) 
CT *( T ) = ^tW^ conc /*( T > u ^ v ^^)- 
We have used the following notations: 



Jo 

where 

A»[iii, u,-, (Tj](^) = <fii(ui(£), Vi(g), cri(£) \ +c E (u Q )a. 

Also, conc/i denotes the concave envelope of the function /j: 

conc/j(r) = inf{/i(s) : /i is concave, h(y) > fi(y) Vi/6 [0, s 4 ]}. 

One can show that the fixed point problem p.20[) admits a unique solution. 

The link between (|1.20p and (|1.18[) is the following: let (u,, Vi, Oi) satisfy (|3.32[) . Assume that Vi < on 
]a, b[ and that Vk(a) = Vk(b) = 0. Define cti(r) as the solution of the Cauchy problem 

dca 1 



dr Vi(r) 
a t {a + b/2) = 
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then (ui o en, Vi o Q!i, (Tj o ai) is a solution to (|3.31[) satisfying 

lim Ui o on (x) = Ui (a) lim it, o olAx) = Ui(b). 

x — > — oo x — >+oo 

Thus, Ui(a) and Ui(b) are connected by a travelling wave profile. 

On the other side, if V{ = on the interval [c, d], then the following holds. Consider TZ\ u(c), the integral 
curve of ri(u) such that TZqu(c) = u{c). Then u(d) lays on TV s .u(c), thus u(c) and u(d) are connected by a 
rarefaction or by a contact discontinuity. 

If Si < 0, one considers a fixed problem like (|1.20p , but instead of the concave envelope of fi one takes 
the convex envelope: 

conv/j(r) = sup{/i(s) : h is convex, h(y) < fi(y) V y e [0, Si]}. 

Again, one can prove the existence of a unique fixed point (u^, Vi, <Ji). 
The curve Tg. is defined setting 

T s >+ := u(si). 

This curve contains states that are connected to u + by rarefaction waves and shocks with speed close to 
Xi(u+). 

If u~ — T].m + , then the limit of the approximation (jl.1711 is 

!U+ X < (Ti(0)t 

Ui(s) X — (Ti(s)t (1-21) 

Ui(si) = T^U + X > <Ti(Si)t 

In the previous expression, is given by (|I.20p and it is a monotone non increasing function. 

It can be shown that in the case of conservative systems with only genuinely non linear or linearly 
degenerate fields the i-th curve of admissible states T*.u + defined in [7] coincides with the one described 
in [39] . Once T*.u + is known, then one defines ip as in (|1.16|) and find the limit gluing together pieces like 



In Section 13.2.21 we give a characterization of the limit of 

E{u e )u e t + A(u e , eu%)u e x = eB(u £ )u% x u e £ R N 

u e (t,0)=u b (1.22) 
u £ (0, x) = uq. 

when the boundary is not characteristic, i.e. when none of the eigenvalues of E~ 1 A can attain the value 
0. The idea is to costruct a locally invertible map <fi which describes all the states that can be connected 
to Mo, in the sense that is specified in the following. Loosely speaking, the map 4> represents for the initial 
boundary value problem what the map ip defined in 1)1.1 6|) represents for the Cauchy problem. Once is 
defined, one takes Uf, as in p.22| and imposes 

Ub = <f>(ua, s i ■ ■ ■ s n)- 

If uq and Ub are sufficiently close, then this relation uniquely determines the values of s± ... s^. Once these 
values are known, then the limit u{t, x) is uniquely determined. More precisely, one can define the value of 
u(t, x) for a.e. (t, x). We will come to this point later. 

The construction of the map 4> works as follows. Denote by Xi(u) . . . \n(u) the eigenvalues of E~ 1 (u)A(u, 0) 
and by r\{u) . . . rjv(tt) the associated eigenvectors. Also, assume that for every u, 

Xi(u) < ... X n (u) < — c < < c < A„+i < . . . Ajv(m) 

In other words, n is the number of negative eigenvalues of E~ 1 (u)A(u, 0). These eigenvalues are real because 
this is one of the hypotheses listed in Section [5] (Hypothesis [3]). 
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For i = n + 1 . . . N consider the i-th curve of admissible states. Fix N — n parameters s n +i ■ ■ ■ sn and 
define 

As before, the notation o T^ N uo means that the starting point of is 

T^uq, T _1 = T^u . Thanks to the results in [7] we quoted before, u is connected to uq by a sequence of 
rarefaction waves and shocks with stricly positive speed. 

To complete the construction, one considers steady solutions of 

E{u)u t + A{u, u x )u x = B(u)u xx 

i.e. couples (U, p) such that 

U ' =P (1 23) 

B(u)p' = A(u, p)p 1 > 

The point (u, 0) is an equilibrium for (|1.23|) . As shown in section 13.2.21 the stable manifold around (u, 0) 
has dimension n, i.e. has dimension equal to the number of strictly negative eigenvalues of E~ 1 A. Also, the 
following holds. Let ip be a map that parameterizes the stable manifold, then ip takes values into M. N x M. N 
and it is defined on a space of dimension n. To underline the dependence on u we will write </>(u, s\ . . . s n ). 
Denote by 7r u the projection 

tt„ :R N xR N ^R N 
(u, p) i— » u 

Fix n parameters si . . . s n and consider tt u o ip(u, s\ . . . s n ). Consider the problem 

A(u, u x )u x = B{u)u xx 
u(0) = tt u o %p{u, si . . . s n ), 

then there exists a unique solution of this problem such that 

lim u{x) = u. 

x — > + oo 

Setting u e {x) :— u(x/e), one finds a solution of 

A(u e , eu%)u £ x = B(u e )u xx 
such that u e (0) — n u o i/j(u, si . . . s n ) and for every x > 0, 

lim u £ (x) — u. 

In this sense, we say that there is a loss of boundary condition when passing to the hyperbolic limit, because 
the boundary condition ir u o ip(u, si . . . s n ) disappears in the limit. We point out that the idea of studying 
steady solutions to take into account the loss of boundary condition was already exploited in many previous 
work, e.g in [25] . 

To complete the characterization of the limit, we define <p as 

4>{uq, si...s n )=it u o T s ™+1 o . . . T* N Z\ o T^ N u , si...s„y 

In section 13.2.21 we prove that the map <fi is locally invertible, i.e. the values of Si ... sat are completely 
determined is one sets 

Ub = 4>{u , si . . . s N ). 

We take the same Ub as in (|1.22p . As pointed out at the beginning of the paragraph, once Si . . . sjv are 
known, then the value of the self similar solution u(t, x) is determined for a.e. (t, x). One can indeed glue 
together pieces like (|1.21j) . In particular, it turns out that the trace of u(t, x) on the axis x = is 

u = T?+\ °...T s N N -_\°T s N N u . 
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As we have already underlined before, the relation between Ub (the boundary datum in (|1.22p ) and u is 
interesting for the study of hyperbolic initial boundary value problems. 



In Section 13.2.31 we give a characterization of the limit of the parabolic approximation 

E{u e )u e t + A(u £ , eu%)u% = eB(u £ )u e xx u e G R N 

u e (t,0) = u b (1.24) 

U 6 (0, x) = Uq. 

when the boundary is characteristic, i.e. one eigenvalue of E~ x A can attain the value 0. The characterization 
of the limit works as follows. We costruct a locally invertible map cj> which describes all the states that can 
be connected to uq. Once <f> is defined, one takes Ub as in (|1.24| and imposes 

u b = 4>{uq, si.. . s N ). 

If uq and Ub are sufficiently close, then this relation uniquely determines the values of s\ ... sn- Once these 
values are known, then the limit u(t, x) is uniquely determined. 

Formally, the idea is the same as in Section 13.2.21 However, the construction of the map 4> is definitely 
more complicated in the boundary chararacteristic case. 

Roughly speaking, the reason is the following. Let Ai (it) . . . Ajv(u) be the eigenvalues of E~ 1 (u)A(u, 0). 
They are real by Hypothesis [31 Assume 

Xx{u) < ... Afc_i(u) < -c < Afc(n) < c < A fe+ i < . . . Ajv(w), 

where c is a suitable positive constant. In other words, there are at least k — 1 strictly negative eigenvalues, 
N — k strictly positive eigenvalues and one eigenvalue close to 0. 
Define 

then Uk is connected to u,q by rarefaction waves and shocks with stricly positive speed. We now want to 
define the k-th curve of admissible states. To define Tjj Uk, we might try to consider the fixed point problem 

u{t) =u k + I r k {u{£), v k (Q, o-fc(O)^ 

Vk(r) = fk(r, u, v k , Ok) - conc/ fc (r, u, v k , a k ) (1-25) 
(Jk ^ r ) = ^tW)d7 C0nc ^ fc ( T ' u ' Wfe ' CTfc )' 

where ffc, fk and ce are the same as in ()1.20p . However, if we consider (|1.25|) we are not doing the right 
thing. Indeed, we might have that the speed Ok is negative at a certain point r. Since eventually we want 
to define the limit u(t, x) as in (|1.2ip . we want <Jk to be greater then 0. 

Another problem is the following. Consider the system satisfied by steady solutions of 

E(u)u t + A{u, u x )u x = B(u)u xx , 

i.e. consider 

U ' =P (126) 

B(u)p' = A{u, p)p [ ° j 

Also, consider the equilibrium point (u,k, 0). Let M. s be the stable manifold of (|1.26[) around (u,k, 0). For 
simplicity, assume \ k (uk) = 0. Then, there might be a solution (U, p) such that 



lim 

: — ^+oo 



U{x), p(x)) = (u k , O) 



but (U, p) does not belong to the stable manifold. However, this kind of solution should be taken into 
account when we study the loss of boundary condition. 
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To tackle these difficulties one can proceed as follows. Instead of the fixed point problem (|1.25p . one 
considers 

tt(r) = u k + / r k (u(£), v k {£), 

Vk(r) = fk(r, u, v k , cr fe ) - mon/ fc (r, u, v k , a k ) (1-27) 

Vkir) = mon/ fc (T, u, v k , a k ). 

In the previous expression, mon/t denotes the monotone concave envelope of the function f k , 

mon/fe(r) = inf |g(r) : g(s) > f k (s) Vs, g concave monotone non decreasing in [0, Sj\ j. 

Some properties of the monotone envelope are discussed in Section 13.1.21 here we stress that mon/fc is a 
concave and non decreasing function, thus the solution a k of ljl.27|) is always non negative. 

Also, the following holds. Denote by (u k , v k , a k ) the solution of (jl.27|) (existence and uniqueness are 
proved in Section 13.2. 3[) . Define 

8 = mm{s : a k (s) = 0} (1.28) 

and 

s = max{s : a k (s) = 0, v k (s) = 0}. 

Assume < s < s < s k . Then u k (s) is connected to u k by a sequence of rarefaction and shocks with positive 
speed. Also, one can show that there exists a steady solution U, 

A{U, U X )U X = B(U)U XX 

such that U(0) = u k (s k ) and 

lim U(x) = u k (s). 

X — > + oc 

However, in general this solution does not belong to the stable manifold of system (|1.26p . This means 
that considering system (|1.27[) we also manage to take into account the converging steady solutions we were 
missing considering just the stable manifold of (|1.26p . 

Heuristically, to complete the construction one should consider the stable manifold of (|1.26|) and hence 
take into account the steady solution that, for x that goes to +oo, converge to u k (s) with fast exponential 
decay, in the sense specified in Section f3.2.3l Actually, the situation is more complex. The reason, loosely 
speaking, is the following . There may be a solution U that converges to u k (s_) and such that some of 
its components converge with fast exponential decay, but other components converge more slowly. This 
possibility is not covered if we consider only the solutions laying on the stable manifold and those given 
by (|1.27p . To take into account this possibility some technical tools are introduced. More precisely, one 
considers suitable manifolds: center stable manifold and uniformly stable manifold. The existence of these 
manifolds is a consequence of results in [33j . but some of the most important properties are recalled in 
Section 

Eventually, one manages to define a locally invertible function <fi. One then takes Ub as in (|1.24p and 
imposes 

Ub = <t>(uo, si . . . s N ). 

If wo and Ub are sufficiently close, then this relation uniquely determines the values of s\ . . . sm- Once these 
values are known, then the limit u(t, x) is uniquely determined and can be obtained gluing together pieces 
like (|1.2ip . In particular, it turns out that the trace of u on the x axis is u k (s), where u k solves (|1.27p and 
s is given by (|1.28p . 

1.1.3 Section[4l the characterization of the hyperbolic limit in the case of a singular viscosity 
matrix 

In SectionH]we discuss the characterization of the limit of (11. ip when the matrix B is not invertible. Actually, 
because of the hypotheses introduced in Section [21 one studies the limit of 

uf + A(u E , eu x )u% = eB(u e )u e xx u e e 

6(u £ (i, 0))=g (1.29) 
u £ (0, x) = uq 
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This system is equivalent to (|1.1|> in the sense specified in Section [2.21 Also, fi = &ya(v £ )J: a precise 
definition is given Section [51 In particular, 6 ensures that the initial boundary value problem (|1.29p is well 
posed. Section [4] is divided into several parts. 

In Section 14.11 we introduce some preliminary results. The point here is the following. In Section [3] we 
give a characterization of the hyperbolic limit when the viscosity matrix is invertible. A key point in the 
analysis is the study of travelling waves 

[A(U, U') - aE(U)]U' = B(U)U" (1.30) 

and of steady solutions 

A(U, U')U' = B(U)U". (1.31) 

To give a characterization of the hyperbolic limit when the viscosity matrix is not invertible, we have to 
study again systems (11.30)) and (|1.3ip . However, being the viscosity matrix B singular, a technical difficulty 
arises. Let us focus, for simplicity, on the case of steady solutions. If B is invertible, we can write 

p' = B(u)- 1 A(u,p)p ( L32 ) 

In this way, we write system (ll.31[) in an explicit form. On the other side, if the matrix B is singular, 
additional work is required to reduce (|1.31|) in a form like (|1.32p . This is indeed done in Section 14.11 What 
we actually obtain is not a 2 A^-dimensional first order system like (|1.32p , but a lower dimensional first order 
system. The exact dimension depends on the structure of the matrix A, in the sense specified in Section [43] 
In Section [4. 2. II we review the characterization of the hyperbolic limit in the case of a Riemann problem, 
i.e. the limit of 

uf + A(u e 7 eu%)u% = sB{u e )u £ xx u £ e R N 
u+ x>0 
u~ x < 



u £ (0, x) 



As for the case of an invertible B (Section 13.2. ip . the key point in the analysis is the description of the 
i-th curve of admissible states T l s .u + . However, there are technical difficulties due to the fact that B is not 
invertible. Actually, in Section 14.2.11 we only give a sketch of the construction, and we refer to [7] for the 
complete analysis. 

In Section 14.2.21 we introduce a technical lemma. The problem is the following. Consider a steady 
solution (|1.3ip and assume that it is written in an explicit form like (|1.32p . This is possible thanks to the 
considerations carried on in Section [4. II Given an equilibrium point for this new system, consider the stable 
manifold around that equilibrium point. For reasons explained in Section [4.2l we need to know the dimension 
of this stable manifold. Lemma 14.71 ensures that the dimension of the stable manifold is equal to n — mi — q, 
where n is the number of strictly negative eigenvalues of A, mi is the number of strictly negative eigenvalues 
of the block An and q is the dimension of the kernel of the block An. The block An is defined by (|1.5p . 
Lemma 14.71 gives an answer to a question left open in [44j . 

In Section 14.21 we discuss the characterization of the hyperbolic limit of 



u% + A(u e , eu%)u% = eB{u E )u £ xx u e e R N 

6(u £ (i, 0)) =g (1.33) 
u s (0, x) = uq 

when the matrix B is singular, but the boundary is not characteristic, i.e. none of the eigenvalues of 
i? _1 (u)A(u, 0) can attain the value 0. 

To provide a characterization, we construct a map <P{uq, s nil + g +i . . . s^) which enjoys the following 
properties. It takes values in and when s nil+q+ \ . . . sm vary it describes all that states that can be 
connected to uq, in the sense specified in the following. Note that similar functions are constructed in 
Section [3.2.21 and Section [3 . 2 . 31 and they are used to describe the limit of the parabolic approximation when 
the viscosity is invertible. However, in those cases the functions depend on N variables, while in the case of 
a singular matrix B the function <fi depends only on (N — Tin — q) variables. 
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Also, in the case of a singular viscosity we have to compose <f> with the function 8, which is used to assign 
the boundary datum in (|1.33[) . It is possible to show that the map 8 o is locally invertible, in other words 
the values of s nil + g +i . . . sjv are uniquely determined if one imposes 

fi o ip(s nil+q+ x . . . s N ) = g, (1-34) 

provided that |fi(lfo) — g\ is small enough. We will plug the same g as in (|1.33p . Once the values of 
s„ 11+9+ i . . . sjy are known, the limit u{t, x) can be determined a.e. (t, x) in the same way as in Sections 
13.2.21 and 13.2.31 In particular, one can determine exactly the value of the trace of the limit u on the axis 
x = 0. As pointed out before, this is important for the study of hyperbolic initial boundary value problems. 

The construction of the map <f> works as follows. Denote as before by n the number of the eigenvalues 
of E~ 1 (u)A(u, 0) that are strictly negative, since the boundary is not characteristic the number of strictly 
positive eigenvalues is N — n. For i = n + 1 . . . N, let T s z . be the i-th curve of admissible states, whose 
construction is reviewed in Section [4.2. II Fix s n+ i . . . sjv and define 

u = T s IX\o...T s n n Z 1 1 oT s n n u . (1.35) 
As before, with the notation Tj^"* o T^uo we mean that the starting point of is T^uq, i.e. 

rrN— 1 rriN - rr,N - 

T o oT s N u o = T SN u . 

The value u is connected to uq by a sequence of rarefactions and shocks with strictly positive speed. 
To complete the construction, we consider steady solutions of 

E(u)u t + A(u, u x )u x = B(u)u xx , 

i.e. we consider solutions of 

A(u, u x )u x = B(u)u xx (1.36) 

In Section [4.1.21 we discuss this system and we explain how to write it as first order O.D.E.. Also, in 
Section 14.2.21 we study the stable manifold of (|1.36[) around an equilibrium point such that u = u and 
u x = 0. In particular, we prove that this manifold has dimension n — nn — q. Let if>(u, s nil + q . . . s n ) a map 
that parameterizes the stable manifold (we are also putting in evidence the dependence on It). For every 
s nil + q . . . s n there exists a solution u of (|1.36p such that u(0) = ip(u, s nil + q . . . s n ) and 

lim u(x) = u. 

x — > + oo 

setting u e (x) := u(x/e) one obtains a steady solution of 

E{u £ )u £ t + A(u e , eu%)u% = B{u £ )u £ xx , 

such that for every x > 

lim u e (x) — u. 

In other words, we experience again a possibile loss of boundary condition from the parabolic approximation 
to the hyperbolic limit. 

The map </> is then defined as follows: 

4>{uq, s nil+q+1 ...s N ) := ipir™- 1 oT* N u Q =T* N u , s„ 11+9+ i . . . s„). 

In Section [4.31 we prove that 8 o <jy is locally invertible with respect to s nil +g+i . . . sjv- 

Here, instead, we point out the following. The fact that the dimension of the stable manifold of II. 361 has 
dimension n — rin — q is proved in Section 14.2.21 and it is not, a priori, obvious. However, n — ri\\ — q is 
exactly the number that makes things works if one wants the function 8 o to be locally invertible, in the 
sense that 8 o </> takes values in R Ar ~™n~9 and hence to be locally invertible it should depend on N — n\\ — q 
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variables. Since we have to take into account N — n curve of admissible states, we are left with n — n\\ — q 
variables for the stable manifold. 



In Section [3.2.21 we provide a characterization of the limit of 

ttf + A(u s , eu x )u% = eB(u 6 )u e xx u e G l w 

fi(u £ (i, 0))=g (1.37) 
u £ (0, x) = uq 

when the matrix B is singular and the boundary is characteristic, i.e one of the eigenvalues of A can attain 
the value 0. Actually this case requires no new ideas, in the sense that one can combine the techniques that 
are described in Section l3.2.3l (to deal with the fact that the boundary is characteristic) and in Section l4.2.3l 
(to deal with the singularity of the matrix B) . 



Finally, in Section [4T3] we introduce a technical lemma which guarantees that the function fio<^> in (|1.34j) is 
locally invertible. We also provide a more general definition for the function 13. Such a definition still ensures 
that the initial boundary value problem (|1.37|) is well posed and that the map 6 o cj> is locally invertible. 

1.2 Unquoted references 

Existence results for hyperbolic initial value problems were obtained in [3T] and [JS] relying on an adaptation 
of the Glimm scheme of approximation introduced in |29j . 

These results were later improved by mean of wave front tracking techniques. These techniques were 
used in a series of papers ([HI [T5] [TH1 HH HOI H3 HH H3]) to establish the well posedness of the Cauchy 
problem. A comprehensive account of the stability and uniqueness results for the Cauchy problem for a 
system of conservation laws can be found in [16j . There are many references for a general introduction to 
system of conservation laws, for example [25] and to |46j . As concerns initial boundary value problems, in pQ 
the existence results in [31] and [45] were substantially improved, while well posedness results were obtained 
in [26] . 

In [5] it is studied the limit of the wave front tracking approximation. In particular, the authors extended 
to initial boundary value problems the definition of Standard Riemann Semigroup. Such a notion was 
introduced in [H] for the Cauchy problem. Roughly speaking, the analysis in [2] guarantees that to identify 
the semigroup of solutions it is enough to consider the behaviour of the semigroup in the case the initial 
and the boundary data are constant. It hence one of the most important motivations for the study of the 
parabolic approximation (|1.1|) . 



Acknowledgments: the authors express their gratitude to Denis Serre for having carefully read a previous 
version of the work and for his useful remarks. Also, the authors wish to thank Jeffrey Rauch for useful 
observations concerning the definition of boundary condition for the hyperbolic-parabolic system. Finally, 
we would like to thank Fabio Ancona for his useful remarks concerning the concave envelope of a function. 

2 Hypotheses 

This section introduces the hypotheses exploited in the work. 

The exposition is organized as follows. In Section [2. II we introduce the hypotheses we impose on system 
(|1.2p when the viscosity matrix B is invertible. In Section l2~2l we discuss the hypotheses exploited when the 
matrix B in (|1.1[) is singular. Finally, in Section |2~31 we introduce the hypotheses needed in both the cases, 
when the matrix B in (|l.ip is invertible and when it is not. 
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In Section [2.2.11 we also discuss three examples which show that, if some of the conditions introduced 
are not satisfied, then there may be pathological behaviors, in the sense specified there. Also, we define the 
function fi used to assign the boundary condition in (jl.ip . 

In [T3] it will be discussed a way to extend to a more general setting the condition of block linear 
degeneracy, which is the third condition in Hypothesis [2] 



2.1 The hypotheses assumed when the viscosity matrix is invertible 

In this section it is considered the system 

v t + A(v, v x )v x = B(v)v xx (2.1) 
in the case the viscosity matrix B(v) is invertible. We assume the following. 

Hypothesis 1. There exists an (invertible) smooth change of variables v = v(u) such that (|2.ip is equivalent 
to system 

E{u)u t + A{u, u x )u x = B{u)u xx , (2.2) 

where 

1. for any u, the matrix E{u) is real, symmetric and positive definite: there exists a constant ce{u) such 
that 

VCeR^, (E(u)£, 0>c E (u)\e- 

2. for any u, the matrix A(u, 0) is symmetric. 

3. for any u, the viscosity matrix B{u) is real and there exists a constant cb(u) > such that 

V(ei N , (B(u% >c B (u)\Z\ 2 . 

In particular, the matrix E{u) may be the differential Dv{u) of the change of coordinates. 
The initial boundary value problem p.2[) is equivalent to 



u 



\ + A(u £ , eu x )u% = eB(u £ )u xx u £ G 
u £ (t, 0) = u b (2.3) 
u £ (0, x) = U , 

where v(uq) = % and v(ub) — v>b- Indeed, one can verify that, thanks to the invertibility of the viscosity 
matrix, it is possible to assign a full boundary condition Ub € M. N . 
In the case system (|2.1[) is in conservation form 

«*+/(«)„ = (£(«)«*), (2.4) 

Hypothesis [T] is guaranteed by the presence of a dissipative and convex entropy. 

For completeness we recall that the entropy rj is dissipative for (|2.4p if for any v there exists a constant 
cd(v) such that 

V£ G R N , (D 2 V (v)£, B(v)0 > c D (vm 2 (2.5) 

It is known (see for example |30j ) that if system (|2.ip admits a convex entropy 77, then there exists a 
symmetrizing and invertible change of variable v = v(u). More precisely, if the inverse function u(v) is 
defined by u = Vt7(u), then v satisfies system (|2.2I) . with A(u, 0) symmetric and given by 

= Df(v(u))(D 2 V (v(u)) 



A(u,0) = Df(v(u))(D 2 T](v(u))) 1 --^(B(v(u))(D 2 r,(v(u))) 

(«) = (D^iv^y 1 



Moreover, 

E 
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is symmetric (being the inverse of an hcssian) and positive definite by convexity, while the dissipation 
condition (|2.5[) guarantees that 

B{u) = B(v(u))^D 2 n(v{u))^ 
is positive definite. Hence Hypothesis [T] is satisfied. 

2.2 The hypotheses assumed when the viscosity matrix is not invertible 

The aim of this section is to introduce the hypotheses that will be needed in Section |4] to study system 

v t + A(v, v x )v x = B(v)v xx (2.6) 
in the case the viscosity matrix is singular. 

Hypothesis 2. There exists an (invertible) smooth change of coordinates v = v(u) such that (|2.6p is 
equivalent to system 

E(u)u t + A(u, u x )u x = B(u)u xx , (2.7) 

where 

1. the matrix B(u) has constant rank r < N and admits the block decomposition 

*(«)=( S ,M ( 2 - 8 ) 



b(u 

with b(u) £ M rxr . Moreover, for every u there exists a constant q,(u) > such that 

V£eR r , (b(u% >c b (u)\d\ 2 . 

2. for any u the matrix A{u, 0) is symmetric. Moreover, the block decomposition of A(u, u x ) correspond- 
ing to (|2.8p takes the form 

' A n (u) A 12 {u) 



A 21 (u) A 22 (u,u x ) y (2J) 

with An e M( N - r ^ x ( N - r \ A 12 e M^-^^, A 21 e M rx(N - r *> and A 22 G MT xr . Namely, only in the 
block A 22 depends on u x . 

3. a condition of block linear degeneracy holds. More precisely, for any <rel the dimension of the kernel 
ker[An(u) — aEn(u)] does not depend on u, but only on a. 

This holds in particular for the limit cases in which the dimension of ker{A\\ —aEn) is equal to N — r 
or to 0. 

4- the so called Kawashima condition holds: for any u, 

ker(B(u)) n {eigenvectors of E~ 1 (u)A(u, 0)} = 

5. for any u, the matrix E{u) is real, symmetric and positive definite: there exists a constant ce{u) such 
that 



In the following, we will denote by 

TPf«A — ( 

E 2 \{u) E 22 {u) 
the block decomposition of E{u) corresponding to 



m = ( S±! SfS ) (2.io) 
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The change of variable v — v(u) guarantees that the initial-boundary value problem (jTTTJl is equivalent 

to 

E{u e )u £ t + A(u £ , £U%)u% = eB(u £ )u% x u e e R N 

fi( W e (i, 0))=g (2.11) 

U £ (0, x) = U , 

where v(uq) — vq, v(ub) — Vb, fi(u e ) = B(i>(it e )). 

Relying on the block decomposition of the viscosity matrix described in Hypothesis in Section 12.2.11 
the explicit definition of the function 6 is introduced. Moreover, in Section [4.31 such a definition is extended 
to a more general formulation. 

In order to justify the assumptions summarized in Hypothesis [5] it will be made reference to the works 
of Kawashima and Shizuta, in particular to [Mj and to [37] . 

In particular, in |37j it is assumed that the system in conservation form 

v t +f(v) x = (S(«)«x) (2.12) 
admits a convex and dissipative entropy rj which moreover satisfies 

(d 2 t)(v)) B(v) t = B(v)(d 2 t)(v)\ • 
If one performs the change of variables defined by w — V 77(11), finds that system (|2.6p is equivalent to 

E(w)w t + A(w)w x = (^B(w)w x S j , 

with A and B symmetric. 

It is then introduced the assumption 

Condition N: the kernel of B{w) does not depend on w. 

and it is proved that Condition N holds if and only if there exists a change of variable w = w(u) which 
ensures that system (|2.12p is equivalent to 

E{u)u t + A(u, u x )u x = B(u)u xx 

for some E(u) that satisfies condition 5 in Hypothesis [5] and some A(u, u x ) and B(u) as in condition 2 and 
1 respectively. 

Moreover, it is shown that Condition N is verified in the case of several systems with physical meaning. 

On the other side, the fourth assumption in Hypothesis [U is the so called Kawashima condition and was 
introduced in [34 . Roughly speaking, its meaning is to ensure that there exists an interaction between the 
parabolic and the hyperbolic component of system (|2.6|) and hence to guarantee the regularity of a solution 

Examples 12.11 and 12.21 in Section [2.2.21 show that, if the condition of block linear degeneracy is violated, 
then one encounters pathological behaviours, in the following sense. One may find a solution of (|2.6|) that 
are not continuously differentiable. This is a pathological behaviour in the sense that, when one introduces 
a parabolic approximation, one expects a regularizing effect. In Section fl. 1.1 1 we explain why it is interesting 
to look for an extension of the condition of block linear degeneracy to a more general setting. This problem 
will be tackled in [T3"] . 

Finally, Example 12.31 show that if the first condition in Hypothesis [5] is violated, then one can have 
pathological behaviours like the one described before. In other words, if the rank of B can vary, then one 
may find a solution of (|2.7|) which is not continuously differentiable. 
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2.2.1 An explicit definition of boundary datum for the parabolic problem 

Thanks to Hypothesis [5] system (II. 2[) is equivalent to 

E(u)u t + A(u, u x )u x = B{u)u xx liGl" 

Q(u)(t,0))=g (2.13) 

u(0, X) EE U , 

where fi(u £ ) = fi[i; e (u)]. In this section we define the function 6. Once fi(u) is known, one can obtain Q(v) 
exploiting the fact that the map v(u) is invertible. In Section [4.31 we extend this definition to a more general 
setting. 

Let r be, as in the statement of Hypothesis the rank of the matrix B. Decompose u as u = (u\, u 2 ) T , 
where u\ G H N ~ r and u 2 G W . then the equation 



can be written as 



E{u)u t + A(u, u x )u x = B(u)u xx 

E n u lt + E 12 u 2t + M\Ui x + A 12 u 2x = , 
E 21 u lt + E 22 u 2t + A 21 u lx + A 22 u 2x = bu 2xx 

The blocks En, An and so on are as in (|2.9[) and (|2.10p . 

In the first line of (|2.14p only first order derivatives appear, while in the second line there is a second 
order derivative u 2xx . In this sense, u% can be regarded as the hyperbolic component of (|2.14p . while u 2 is 
the parabolic component. As explained in section fl.1.11 one can impose r boundary conditions on u 2 , while 
one can impose on u\ a number of boundary conditions equal to the number of eigenvalues of E^ An with 
strictly positive real part. 

We recall that we denote by nn the number of strictly negative eigenvalues of An, by q the dimension 
of kerAn and by n the number of strictly negative eigenvalues of A. One can prove that the number of 
eigenvalues of E^An with strictly negative real part is equal to nn (see Lemma (|3.1|) in Section 13.1. 1[ 
which was actually introduced in [B]). Also, the dimension of the kernel of E^An is q. 

Let Ci(it, 0) G l w_r be an eigenvector of E^ 1 An{u) associated to an eigenvalue rji with non positive real 
part. Let Z z (u, 0) G R N be defined by 

Ci ■ (2.15) 



and finally let 

Z(u) := span(Zi(w, 0), . . . , Z„ 11+q (u, 0)}. 

Finally, let 

\ { \ { Wi \ ( W N - r - nil - q 



where e*j G W are the vectors of a basis in W and Wj G R N_r are the eigenvectors of E^An associated to 
eigenvalues with strictly positive real part. 

Since W(u) © V(u) = R N , every u G I 1 *' can be written as 

u = u w +u z , u w € W(u),u z G Z{u), (2.16) 

in a unique way. 

We define fi as follows. 



Definition 2.1. The function 6 which gives the boundary condition in (|2.1ip is 

fi :R N -> R N - n ^-i 
u i ^ u z , 

where u z is the component of u along W(u), according to the decomposition (|2.16p . 
We refer to Section |4~31 for an extension of this definition. 



(2.17) 
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Figure 1: the graph of the function 1*1(2) in Example 1 2. II 

| tti 



«i(0) 




«i(0) 



2.2.2 Examples 

Example 2.1. This example deals with the equation 

A(u, u x )u x = B(u)u XXl (2-18) 

in the case the condition of block linear degeneracy (the third in Hypothesis [2J does not hold. 

System ()2.18j) is satisfied by the steady solutions of (|2.7p . One therefore expects the solution of Q2.18J) 
to have good regularity properties: the example which is going to be discussed show that if the condition of 
block linear degeneracy does not hold a function u satisfying (|2.18j) may have the graph illustrated by Figure 
[T] and hence be not C 1 . Moreover, the figure suggests that a pathological behavior typical of the solution 
of the porous-media equation may occur, namely it may happen that a solution is different from zero on a 
interval and then vanishes identically. Let u = (u\, i*2) T and let 



B{u) = 

and 

A(u, u x ) = 




1 

iti 1 
1 0, 



Then the assumption of block linear degeneracy is not satisfied in a neighborhood of u\ = 0: we therefore 
impose the limit conditions 

lim Ui(x) = lim Ui{x) = 0. 

x — >+oo x — >+oo 

In this case equation (|2. 18[) writes 

uAu lx + l) = 



(2.19) 

iti = u 2x 

Since the equation satisfied by u\ admits more solutions in a neighborhood of u\ = 0, to introduce a selection 
principle we proceed as follows. 

We consider a solution u v = (u", u^) T of (|2.18p such that 

lim Ui(x) = v. 

x — >+oo 

The parameter v is positive and we study the limit v — > + . Fix the initial datum u\§ > 0. The component 
i*i satisfies the Cauchy problem 



lx u'i (2.20) 
«?(0) = u?o 
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If v is sufficiently small then uio > v and hence u\ x is always negative. In other words, u\ is a monotone 
non increasing function which satisfies v < u^(x) < u\q for every x. Also, note that if v\ < v 2 then, for 
every x, 

<(ar) < v%{x). (2.21) 

This can be deduced by a comparison argument applied to the Cauchy problem (|2.20|) . Indeed, if u\ > 
then 

V — H, 2/ — Mi 
77T- < TIT-- 



In particular, from (|2.21[) we deduce that for every x u\{x) is monotone decreasing with respect to v and 
hence it admits limit v — > + . 

Denote by u\ the pointwise limit of u\ for v — > + : we claim that ui satisfies 

01(1*13 + 1) = 0. 

Indeed, let 

Xq '■= min{a; : u (x) = 0}. 

If x < xq, then by monotoniticity for every v and for every y < x u\a > u"{y) > u\(x) > 0. Also, if v is 
sufficiently small then u"(x) > and by monotonicity u\(y) > u\(x) > 0. Consider the relation 

u 1 (x)=u w + dy. 



<{y) 

We take the limit v — > + and, applying Lebesgue's dominated convergence theorem, we get 

u\(x) = u%o — x. 

On the other side, if x > xq then by monotonicity U\(x) < 0. On the other side, u\(x) > for every x. We 
conclude that Ui(x) — if x > xq. In other words, xq — itio and 



u^x) 



uiq — x x < uio 
x > u\o- 



This function ha the graph illustrated in Figure [T] and it is not continuously differentiablc. In this sense, we 
encounter a pathological behaviour. 

Remark 2.1. An alternative interpretation of the previous considerations is the following. 
Consider the family of systems 

UlUlx + u 2x — 

Ulx = U 2 xx (2.22) 
iii(0) = v lim^^+oo m(x) = y lim x ^ +00 u 2 (x) = 

parametrized by the limit value y G M. Let U\(y) be the set of the values v such that (|2.22p admits a 
solution. The ODE satisfied by m is 

y-u\ 

U\x = 

Ml 

when u\ ^ 0. From a standard analysis it turns out that when y > 0, then U\(y) =]0, +00 [. When y = 0, 
Wi(0) = [0, +00 [, while when y < 0, Wi(y) = {y}. 

In other words, IA\ is a manifold of dimension 1 when y > 0, a manifold with boundary and with dimension 
1 when y = 0, while when y < the dimension of the manifold drops to zero, i.e. the manifold reduces to a 
point. 
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Figure 2: the graph of the function u\{x) in Example 12.21 

"1 



ui(x) 



Example 2.2. Example 12.11 shows that, if the condition of blcok linear degeneracy (the third in Hypothesis 
^ is not satisfied, then one can find a steady solution of 

E(u)u t + A(u, u x )u x = B(u)u xx (2.23) 

which is not continuously differentiable. 

As we see in the following sections, steady solutions are important when studying the limit of the parabolic 
approximation (|2. 13(1 . Other solutions that play an important role are travelling wave solutions, i.e U such 
that 

\A(U, U') - aE(U)j V = B(U)U". (2.24) 

In the previous expression, <r is a real parameter and it is the speed of the travelling wave. More precisely, 
in the following sections we study travelling waves such that a is close to an eigenvalue of the matrix A. 
Let 






A(u, u x ) := | 1 1 | B{u) := ( 1 | , (2.25) 

y 1 

which are obtained from a system in conservation form 

Ut + f(u) x = (b(u)u x J 

where 

f(u) = (u\/2 + u 2l ui +u 2 , o) 

and B(u) is the constant matrix defined by (|2.25p . 

The matrix A defined by (f2~23j) has an eigenvalue identically equal to zero: we will focus on the travelling 
waves with speed a = 0. 

System (|2.24p can then be rewritten as 

U\Ui x + u 2x = 

uix + u 2x = u 2xx (2.26) 
= u 3xx . 

As in the case considered in Example 12. 1[ the condition of block linear degeneracy is not satisfied in a 
neighborhood of u\ = and therefore we study travelling wave solutions such that 

lim ui(x) = lim u 2 (x) = lim ^3(2;) = 0. 

x — >+oo x — >+oo X — >+oc 
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From system (I2.26[) and from the previous condition one obtains 



u\/2 + u 2 = 

Ul+U 2 = u 2x 

u 3 = 0. 

Hence in the following we ignore the third component of the solution and we study only the first two lines 
of the system, which can be rewritten as 



u 2 = —u\/2 
u-\U\ x — u\/2 — u 



2/0 .. ( 2 - 27 ) 



If u\ 7^ 0, the second line is equivalent to 

Ulx = Ml/2 - 1. 

Fix uia such that < u±q < 2. We impose ui(0) = u±o and we obtain 

ui(a;) =2+{u w -2)e x ' 2 , 

This solution can be extended as far as u\ > 0: when u\ reaches the value zero, (|2.27[) admits more that one 
solution. 

In order introduce a selection principle, we proceed as in Example 12.11 and we introduce a family of 
solutions u v = (ul, u%, u%) of ([2~24]) such that 

lim u"(x) = v lim u 2 (x) = lim 1*3(0;) = 0. (2.28) 

x — >+oo x — >+oo X — >+oc 

We also impose Mi(0) = uya- In the previous expression, v is a small and positive parameter and we study 
limit v — > + . One can repeat the same considerations as in the previous example and conlude that when 
v — » + the solution u\ converges pointwise to 

f u x {x)=2+ MO) - 2)e x / 2 x<2\og (2/(2 - u 10 )) 
ui(x) = { V , ; {, (2.29) 

[ x>21og(2/(2- Ul o)J 

This function has the graph illustrated in Figure [5] and it is not continuously differentiable. 



Example 2.3. The aim of this example is to show that if in a system of the form (|2.7j) the rank of the matrix 
B is not constant, thus contradicting the first assumption of Hypothesis fJl then pathological behaviors of the 
same kind described before may appear. More precisely, we find a steady solution which is not continuously 
differentiable. 

We consider the system in conservation form 

u t + f{u) x = (b(u)u x ^J 

with 

/(«!.«») = ( S) B( Ul ,u 2 )=^f" I 

for some 7 > 3. In this case, the rank of B(u) drops from 2 to 1 when u\ reaches 0: to find out pathological 
behaviors it seems therefore natural to study the solution in a neighborhood of u\ — 0. More precisely, the 
attention is focused on steady solutions 



such that 



f(u) x = (-B(uK) (2.30) 



lim u^x) = mix) > lim u 2 (x) = (2.31) 

X — » + oo X—^ + OO 
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and it will be shown that the graph of the first component u\ has the shape illustrated in Figure [3] Hence 
i*i is not continuously differentiable and presents a behaviour like the one typical of the solutions of the 
porous-media equation: it is different from zero on a interval and then vanishes identically. 
In this case equation (12.30|) writes 



«2; 



Ml = U 2x 



and hence after some computations one obtains 



Mil 



"2 



,(3-7) 



7(7 + 1) 1 
ui(y)dy. 



-t-oo 



The equation satisfied by u\ admits more than one solution in a neighborhood of u\ = 0. 
To introduce a selection principle we consider the matrix 



B v (u v ) 



and the equation 



f + 7K) 7 ~ 1 o 
1 



(2.32) 



(2.33) 



(2.34) 



In the previous expression, v is a positive parameter and we study the limit v — > + . We keep the limit 
conditions (|2.3ip fixed. In other words, this time we do not perturb the limit condition (as in the previous 
examples), but the equation itself. Note that now B" is invertible. The solution of (12.341) with limit 
conditions (|2.31[) satisfies 
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K) 



!/\7 + l 



Indeed, from 



7+1 



(2.35) 



one obtains integrating 



^L+7(«i) 7_1 «i 



Multiplying the first line by u\ — u 2x and then integrating again one obtains (|2.35|) . From (|2.35|) we get 

(2^ + 7K) 7 fe 



>2.v 



27 
7 + 1 



^7+1 



Taking = one eventually gets 



"2^T7K)^ ' 



Fix ttio > and consider the Cauchy problem 



<(o) = «r 
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Figure 3: the graph of the function u\{x) in Example 12.31 



ui(0) 




> 



x 



One can then exploit the same considerations as in the previous examples: in particular, the fact that for 
every x Ui{x) is monotone decreasing with respect to v follows from 



if 7 < 3. One can the prove (proceeding as in the previous examples) that when v — > + , u\ converges 
pointwise to a function u\ satisfying 



This function has the graph illustrated in Figure [3] and it is not continuously differentiable. 
2.3 General hypotheses 

This section introduces the general hypotheses required in both cases, i.e. when the viscosity matrix B in 
(jl.ip is invertible and when it is singular. In the statements of the hypotheses we actually make reference to 
the formulation (|2.ip and (|2.13|) . Thus, in particular, we consider the same value uq as in ()2.ip and (|2.13| . 
First of all, we assume strict hyperbolicity: 

Hypothesis 3. There exists 6 > such that, if u belongs to a neighbourhood of uq of size 5 then all the 
eigenvalues of the matrix E~ 1 (u)A(u, 0) are real. Moreover, there exists a constant c > such that 





where 




inf {\Xi(u, 0) - Xj(u, 0)|} > c> Vi^j. 



We also introduce an hypothesis of convergence: 



Hypothesis 4. Let 



E(u £ )u £ t + A(u £ , eu%)u e x — eB(u e )u% 
Q(u £ (t,0))=g(t) 
u £ (0 7 x) = uq(x) 



'xx 
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be a parabolic initial boundary value problem such that g, Uq £ Lj oc and 

\g(0) -fi(tZ (0))|, TotVarjg}, TotVar{u } < 5 (2.36) 

For a suitable constant 5 << 1. Then we assume that there exists a time To, which depends only on the 
bound 6 and on the matrices E , A and B, such that 

TotVar{u e (t, •)} < CS Ve, t e [0, T ]. 

We also assume directly a result of convergence and of uniqueness of the limit: 

Vt S [0, T ] u £ (t, •) -> u{t) L\ oc when s -> 0+ TotVar{w(t)} < C<5. 

We point out that the uniqueness of the limit is actually implied by the next hypotheses and the unique- 
ness results for the Standard Riemann Semigroup with boundary: it is made reference to [5] for the extension 
of the definition of SRS to initial boundary value problems, while an application of this notion to prove the 
uniqueness of the limit of the vanishing viscosity solutions was introduced in [llj in the case of the Cauchy 
problem. 

It is also assumed: 

Hypothesis 5. It holds 

|fi(uo) -g\<5 

for the same constant 5 << 1 that appears in Hypothesis^ 

An hypothesis of stability with respect to L 1 perturbations in the initial and boundary data is also 
introduced: 

Hypothesis 6. There exists a constant L > such that the following holds. 
Two families of parabolic initial boundary value problems are fixed: 

uY + A{u ls ,eul s )ul s = eB{u^)ull ( u 2 - + A(u 2e , eu^)u 2 / = eB(u 2 °)u 2 x Z 

8(u^(t,0))^g\t). \ Q(u^(t 7 0))=f(t). (2.37) 

M le (0, x ) = ul(x) { u 2£ {0, x) = ul(x) 

with Uq, g 1 and Uq, g 2 in L\ oc and satisfying the assumption (|2.36p . 
Then for all t £ [0, To] and for all e > it holds 

2E/J-MI x r | ll-l -2|| ll-l -2| 



\\u^(t) - u 2 %t)\\ L1 < L(JK - u 2 \\ L1 + Us 1 - g 2 \\ L1 

From the stability of the approximating solutions and from the Lj oc convergence one can deduce the 
stability of the limit. More precisely, let u 1 and u 2 the limits of the two approximations defined above, then 



\\u\t)-u\t)\\ L ,<LUul-ul\\ L i + \\f-f\ 



L 1 



Finally, it is assumed that in the hyperbolic limit there is a finite propagation speed of the disturbances: 

Hypothesis 7. There exist constants [3, c > such that the following holds. 

Let Uq, g 1 and u 2 , g 2 in Lj oc and bounded, let u el and u e2 be the solutions of (|2.37|) and u 1 and u 2 the 
corresponding limits. If 

g 1 (t)=g 2 (t) Mt<t ul{x) = ul{x) \/x<b 

then 

\u £ \x, t)-U E2 (x, t)\ < 0^ e -o»>ta{k-»>»{0,/3(t-to)}|,|x-6+/Jt|}/e^ y x £ [ max { 0j (3(t - t a )} , b - f3t}. 

Analogously, if 

Uq(x) — u (x) V x £ [a, b] 

then 

\u £l ( Xl t)~u e2 {x, t)\ < o(e- cmia i ]x - a - 0t] ' ]x - b+0t] } /e } Vie [a + pt,b-[3i\. 
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Figure 4: the finite propagation speed in the hyperbolic limit 




0(t -t ) b-Pt b 



a a + (3t b- (it b 



Remark 2.2. Hypothesis [7] may appear a bit technical. It implies the finite propagation speed of distur- 
bances in the hyperbolic limit: an heuristic representation of this phenomenon is illustrated in Figure [5] 
Loosely speaking, the reason why we need finite propagation speed is the following. We have to be sure that 
the limit of (|2. 13|1 in the case of a generic couple of data (uq, g) can be obtained gluing together the limit 
one obtains in the case of cooked up data, namely data connected by travelling wave profiles. 

3 Characterization of the hyperbolic limit in the case of an in- 
vertible viscosity matrix 

The aim of this section is to provide a characterization of the limit of the parabolic approximation (jTTTJ) 
when the viscosity matrix B is invertible. The precise hypotheses that are assumed are listed in Section [2~T1 
and 12.31 in particular, these hypotheses guarantee that it is sufficient to study system 

E(u e )uf + A(u e , eu x )u% = eB(u £ )u% x u £ S R N 

u £ (t,0)=u b (3.1) 

U E (0, X) = Mo, 

where the matrices E, A and B satisfy Hypotheses [T] and [31 

The exposition is organized as follows. In Section T3. II we discuss some preliminary results. More precisely, 
in Section IS. 1.1 1 we recall some transversality results discussed in [^. In Section 13. 1.21 we recall the definition 
of monotone envelope of a function and we introduce some related results. In Section 13.21 we review some 
results in [7j . Namely, we give a characterization of the limit of parabolic approximation 

f E{u £ )u £ + A{u £ 2 eu%)u% = eB{u £ )u% x 

1 »'<»-)={; :;o. 

In Section f4.2.3l it is given a characterization of the limit of the parabolic approximation (|3.1[) in the case 
of a non characteristic boundary, i.e. when none of the eigenvalues of E~ 1 A can attain the value 0. The 
case of a characteristic boundary occurs when one of the eigenvalues of E~ 1 A can attain the value and it 
is definitely more complicated than the previous one. The characterization of the limit (|3.1j) in the case of 
a characteristic boundary is discussed in Section 13.2.31 

3.1 Preliminary results 
3.1.1 Transversality results 

The following lemma is discussed in [28j . However, for completeness we will repeat the proof. 



2G 



Lemma 3.1. Assume that Hypotheses^ and\3\hold. Then for any given u and u x 

1. B- l {u)A{u, u x )(= <^=> A(u, u x )^ = <=> E- 1 (u)A(w, u x )£= 

2. the number of eigenvalues of the matrix B~ 1 (u)A(u, 0) with negative (respectively positive) real part is 
equal to the number of eigenvalues of E~ 1 (u)A(u, 0) with negative (respectively positive) real part. 

Proof. The first point is actually an immediate observation. 

To simplify the notation, in the proof of the second point, we will write A, B and E instead of A(u, 0), 
B(u) and E(u) respectively. One can define the continuous path 

F : [0, 1] ^U NxN 

sh(1- s )B + sE, 

which satisfies the following condition: for every s G [0, 1], F(s) is positively definite and hence invertible. 
Indeed, 

V|*G R N (F(8)Z, > mm{c B (u), c E (u)}\£\ 2 . (3.2) 

Moreover, by classical results (see for example [38]) from the continuity of the path F(s)^ 1 A it follows the 
continuity of the corresponding paths of eigenvalues Ai(s) . . . Ajv(s)- 

Let k — 1 be the number of eigenvalues of F(0)~ 1 A = B~ 1 A with negative real part: 

i?e(Ai(0)) < . . . ite(A*_i(0)) < = ite(A fc (0)) < Re(X k+1 (0)) <... Re(X N (0)). 

Because of the continuity of Ai(s) . . . Xi(s), to prove that the number of negative eigenvalues of F(1)~ 1 A = 
E~ 1 A is (k — 1) it is sufficient to prove that for any % = 1 . . . (k — 1), Aj(s) cannot cross the imaginary axis. 
Moreover, the case Ai(s) = is excluded because otherwise the corresponding eigenvector £ should satisfy 
F(s)~ 1 A^(s) = for any s and hence Xi(s) = 0. It remains therefore to consider the possibility that Aj(s) 
is purely imaginary. 

By contradiction, assume that there exist s G [0, 1], A G K \ {o}, r G such that 

F(s)- 1 Ar = iXr : (3.3) 
just to fix the ideas, it will be supposed that A > 0. Moreover, let 



Then from (13. 3p one gets 



r = n + ir 2 , ri, r 2 G R N . 



Ar x = -XF(s)r 2 
Ar 2 = XF(s)ri 

and hence from (|3.2[) and from the symmetry of A it follows 

(Ari, r 2 ) < -Amin{cB(-u), c B (w)}|r 2 | 2 

(Ari, r 2 ) = (Ar 2 , ri) > Amin{c B (u), c_E(u)}|ri| 2 , 

which is a contradiction since A > 0. This ends the proof of Lemma I3TT1 □ 



In the following, given a matrix C we denote by V s (C) the direct sum of all eigenspaces associated with 
eigenvectors with strictly negative real part, by V C (C) the direct sum of the eigenspace associated with the 
eigenvector with zero real part, and V U (C) the the direct sum of all eigenspaces associated with eigenvectors 
with strictly positive real part. Lemma 13-11 ensures that, if Hypothesis [T] holds, then for all u, u x 

V c (E- l (u)A(u, 0)) = V c (A(u, 0)) = V c (B- l (u)A{u, 0)) 
dimV u (E- 1 {u)A(u, 0)) = dimV u (B~ 1 (u)A(u, 0)) 
dimV^E^i^Aiu, u x )) = dimV s (B _1 (u)A(u, u x )) 
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Lemma 3.2. Assume that Hypotheses^ and\^hold, then 

V u (e- 1 (u)A(v, 0)) nV B (B-\u)A(v, 0)) ={0}. 

Thanks to the previous lemma, which is a direct application of Lemma 7.1 in [6], one can conclude that 
actually 

V u (e- 1 (u)A(u, 0)) © WiT 1 ^)^ 0)) © W-B-^uJAfu, 0)) = R N . (3.4) 
3.1.2 Some results about the monotone envelope of a function 

The aim of this section is to collect some results that will be needed in Section 13.2.31 Proposition 13.11 is 
very similar to results that were introduced, in a much more general form, for example in [5] and |32j . 
However, since the situation discussed here is slightly different, for completeness we give a proof. Also, 
results analogous to Propositions 13.21 and 13.31 arc discussed in [4], but a proof is given here for completeness. 

In the following, conc[ , s ]/ will denote the concave envelope of the function / on the interval [0, s], 
namely 

(conc[ , s ]f^j (t) := inf (h(r) : h(t) > f(t) Vf £ [0, s], h is concavej. (3.5) 

We will consider only the case of a function / £ C^' 1 ([0, s]), i.e. / £ C 1 ' 1 and /' is Lipschitz continuous with 
Lipschitz constant smaller or equal to k. 

The symbol mon[Q S ]/ will denote the monotone envelope of the function / on the interval [0, s], i.e. 

^monjo, g ]/J (t) := sup |/i(t) : h(t) > f(t) V< £ [0, a], h is concave and nondecreasingj. (3.6) 

Lemma 3.3. Let f £ C^QO, s]), then 

1. conC[ j/ is continuous and satisfies ||concro )S i/||c° < |]/]]c°- 

2. 7/r £]0, s[ and 

(coiic[o, s ]/)(t) = /(r), 
£/ien concro iS ]/ is differentiable at t and 

(cone [0)S] /) (r) = /'(r). 

3. 7/r £]0, s[ and 

rconc[ , s ]/)(r) > /(r), 
then there exists a, b such that 0<a<r<b<s such that 

Vt £ [a, 6], (conc [0 ^]/J(i) = * 

^. conc [0iS ]/(0) = /(0) ; conc [0)S ]/(s) = f(s). 
Proof. By definition, 

f(r) <conc [0 , fl] /(r) Vr e [0, s] (3.7) 

and hence — ||/||e° < concm s i/(t). Moreover, the constant function ||/||co is concave and greater or equal 
to / and hence conc[ 0j s ]f( T ) < ll/llc°- Combining the two inequalities, one gets ||conC[o, s ]/||.L°° < ll/llc°- 
Moreover, conc[ s ]/ is a concave function and hence it is continuous in all the inner points of its domain, 
thus in ]0, s[. Later on we prove that it is also continuous at t — and t = s. 
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To prove the second point, assume that r G]0, s[ and that conC[ , s ]/(T) = /("?")• If h > 0, then thanks to 

fjr + h) -/(t) conC[ Q , s ]/(r + fe) - conc^ s] /(t) 
h h 

and hence 

f (r) = Um /( r + fo ) ~ /( r ) < Um conc [Q^]/( r + h)- conc [0; s] /(t) 

Note that /'(r) exists because / eCj 1 , while 

concr 0jS ]/(r + /i) - concr , s ]/(r) 

iim ; — ; ; — 

h^0+ h 

exists because the difference quotient is a monotone non increasing function. By the same reason, 
j. c onC[o, s ]/( T + h) - conc [0;Ji ] /(t) < conc [0iS] /(r + h) - conc [0; s] /(r) 

h^O- /l _ h— >0- /l 

Moreover, taking /i < one can proceed as before and get 

/'( r ) = Mm tt T + h )-tt T ) > Um conc [M /(T + /Q-conc [(M /(T) 

and hence putting all the previous considerations together one gets that (concjo s ]f)' (t) exists and it is equal 
to/'(r). 

To prove the third point, take r such that concr 0) s ]fi T ) > f( T ) an d take a and 6 as follows: 

a := mf{f < r : f(t) < conc [0 , s ]/(i)} 

and 

b := sup{i > r : f(t) < conc [0 , s ] /(*)}. < s 
To simplify the notations in the following we will also assume 

conc [0 , s ]/(a) = conc[ 0i s] f(b) = 0, 

the general case being analogous by subtracting a linear function. 

First of all, note that conc[ 0j s ]/(i) > for every t G [a, b] because conC[ 0iS ]/ is a concave function 
and hence it is greater or equal to the segment joining two of its values. We want to prove that also 
conqo, s ]f(t) < 0. More precisely, we will consider the function 

!conc [0jS ]/(i) t<a 
a<t <b 

conc [0>s] /(t) t > b. 

and we will prove that conqo, s ]f(t) < h(t). Since the function h is concave, it is enough to show that 
h(t) > f(t) for every t G [a, b]. By contradiction, assume that there exists r G [a, b] such that /(r) > 0. 
Assume that r is exactly a point at which / assumes it maximum on [a, b\. If r = a or r = 6 then there 
is nothing to prove because /(i) < for every t in [a, 6]. Moreover, let r m denote a point at which the 
maximum of coneys] / on [a, b] is assumed: since conC[ 0jS ]/ is concave, then it is non decreasing for t < r TO 
and non increasing for t > r m . Moreover, by (13. 7|) /(r) < conqo, s ]f{ T m)- Consider the function z defined as 
follows: 

!COnC[o, a ]/(*) t < t a 
/(r) a<t<b 
conc[ , s ]/(£) t>t b , 

where 

t a = min{f G [a, 6] : concp, s ]/(t) = /(r)} 
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and 

t b = ma,x{t e [a, b] : conc [0;S ]/(i) = /(r)}. 

Putting all the previous considerations together it turns out that t a < r m < tb and hence z is a concave 
function which is greater or equal to /. Thus, conc[o. s ]/ < z and hence in particular 



cone 



[o,.]/(t)<z(t) = /(t). 



This contradicts the assumption conC[ 0;S ]/(t) > /(t) on ]a, 6[. 

Note that if a > 0, then automatically /(a) = conC[ 0;S ]/(a). Now we want to prove that, even if a = 0, 
/(0) = conc[o. s ]/(0). By contradiction, assume that /(0) < conc[ 0i s ]/(0) = 0, then consider the function 

r /(o)+ii/'ii C ot t<t 

w(t) < a<t<b 
{ conc[ 0iS ]/(i) f > i, 

where f is the point at which /(0) + ||/'||co* = 0. Then w is a concave function greater or equal to / and 
hence w > conc^]/, which contradicts the assumption /(0) < conc[ 0j s ]/(0). 

In an entirely similar way one proves that /(&) = conc[ 0j s ]/(^) even if & = s. This completes the proof of 
the third point. 

To prove the fourth point we proceed as follows. From the points 2 and 3 we know that conqo.s]/ is a 
bounded function that satisfies 

(conc [0 , s] /)'(r)| < Wf'Wca Vre]0, s[ 

Thus, it has bounded total variation and hence the following limit exists: 

lim conc[ . s] /(x) := L. 

We want to prove L — /(0). Let 

S := {r G]0, s[: /(r) < conc [0! s] /(r)}. 

The set S is open and hence 

n>l 

for suitable sequences {a„} and {&„}. 
If 

inf{a„} = 0, 

there are two possibilities. If the infimum is actually a minimum, then 

/m = m_m v T£| o, 6[ 

for some b > 0. One can then exploit the same function w as before to prove that /(0) = L. If there exists 
a subsequence a nk — > + , then 

L= lim conc [0iS] /(a n J= lim /(a„J = /(0). 

Finally, if 

inf{a„} > 0, 

then 

/(t) = conc [0iS ]/(r) Vr e]0, inf{a„}[, 

then by the continuity of / L = /(0). 

Analogous considerations guarantee that conc[ s ]/ is continuous at t = s and that concp s ]/(s) = 
/(*)■ □ 
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Because of points 2 and 3 of Lemma 13.31 the function [conc^. s]f j is continuous at every point in ]0, s[. 
Also, relying on considerations analogous to those performed in the last part of the proof of Lemma 13.31 one 
can prove that ^conc^.s]/^ is continuous at t = and t = s. Note that it is actually enough to prove that 



the following limits exist: 



^lim + (coneys]/) (r) ^lim_ (conc [0i s] ./) (r). 



Indeed, one can then apply the theorem of the limit of the derivative. We conclude that the following holds 
true: 

Proposition 3.1. If f € C k ' ([0, s]) , then conc[ , s ]/ G C fc ' ([0, s\). 

The following result describes the relation between the concave and the monotone envelope. The proof 
exploits considerations similar to those used to prove Lemma l3.3l and it will be therefore omitted. 



Lemma 3.4. Let f G C\ A {[Q, s]), th 



vn. 



mOn [0, S ]/(T) 

The value tq is given by 



conc [0 , s ]/(r) if r < t 
conc[ , s ]/(r ) if r > t . 



t := max^t G [0, s] : (conC[ 0j s ]f)'(t) > oj. (3.8) 
// (conqo, s ]f)'{t) i s negative for every t, then we set tq = 0. 

Combining Lemma I3~4l and Proposition 13.11 one deduces that if / belongs to C^QO, s]) then 

mon [0 , s] / G C^QO, si]). 
The following propostion collects some estimates that will be exploited in Section f3. 2. 31 



Proposition 3.2. Let f, g G C£ ([0, s]). Then 
1. 

||conc[ , s ]/||co < ||/||c° ||(conc[o, s ]/)'|| c o < ||/'||c°- 

2. 

||conc [0>s] /-conc [ o ;S ]S'||co < \\f - g\\ c ° II (conc [0>s] / - concp^^'Hco < \\f - g'\\c° (3-9) 

Proof. To prove the first point, note that 

conc[ , s ]/(t) < conc [0 , s ](/ - g)(t) + conc [0 , s ] £/(t) Vt 

because conc[ 0i s ] (/ — g)(t) + concjo, s ]g{t) is a concave function greater or equal to / — g + g = f. Thus, by 
the first point in Lemma \3. 31 

conc [0 , s ]/(t) - conc [0!s] g(f) < ||conc [0! s] (/ - g) || c o < ||/-#|| c o Vt G [0, s] 

and by analogous considerations 

conc [0: s] g(t) - conc [0)S] /(t) < \\f - g\\ c « Vi G [0, s] 

This concludes the proof of the first point. 

To prove the second point, let us first make a preliminary consideration: 

(concjo, s ]g)\s) 7^ g'i s ) =^ 3q<s: (conc[ , s ]<?)'(i) = constant Vi G [a, s] (3.10) 
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Indeed, if there is t n — > s such that g(t n ) = conc[ 0j s ]9(tn) then g'{t n ) — (concm ^g) (t n ) and hence 

(concr 0j s ]g)'( s ) = g'( s )- I n particular, if (conc[ 0j s ]ff)'( s ) 7^ ff'( s ) then conc^^g > g on an interval [a, s] and 
thus exploiting Lemma 13.31 point 2, one conclude with (|3.10[) . 

Let t be a point at which the maximum of |(conC[ 0jS ]/ — concp, s ]3) I is assumed. Just to fix the ideas, 
suppose 

(conc[ , s ]/ - conc[ 0)S ]g) (r) > 0. 

Because of the second point in Lemma l3~3l if (concm s i/) (t) ^ /'(t), then there are a, b such that r £ [a, 6] 
and 

(conc [0 , fl] /)'(t) = f[a) b Z f a [b) Vie [a, 6]. 

Moreover, since (conc[ 0) s ]g)' is non increasing, it is not restrictive to suppose t — b. Also, because of the first 
point in the statement of Lemma [531 if b < s then (conc[ 0j s ]f)'{b) = f'(b). Thus, if (conc^.s] f)'(b) ^ f'{b) 
then b = s. If (conc[ 0; s ]<?)'( s ) = <?'( s ) then 

(conc [0 , s ]/ - conc[ , s ]5) (s) < (/ - g) (s). 

Indeed, 

f(s) = conc [0 , s ]/(s) conc [0 , s ]/(i) > /(f) Vt < b 

and hence (conc[ 0i s ]/)'(s) < f'{s). 

On the other side, if (concjo, s ]ff)'(s) ^ then combining (|3.10p and the second part of Lemma f3.3l 
one gets that there exists a < s such that 

conc [0 , s] g s = yy ' . 

s — a 

Moreover, a > a. Indeed, if a < a then (conc^ s ]f)'{ a ) > ( conc [o. s]f)'(b) and hence & would be no more a 
point at which the maximum of (concp s j/ — conc^, s ]9) 1S assumed. By Lagrange's theorem, there exists 
x e [a, s] such that 

(/ - .9) (») - — 

[/(s) ~ ~ [conc [0!s] /(a) - g(a)] 



> 



s — a 
f(s) - concp, s ]f (a) 



cone 



s — a 

(cone 



[0,s]9) (i 



[o,s]f) (x) - (conc [0 , s ]) (a 



This concludes the proof of the second point in the case (conc^ )S ]/) ( r ) > ( conc [o. s]/)'( r )- The opposite 
case (conc[ 0i s ]/)'( r ) < ( conc [o,sl/) ( T ) is entirely analogous. □ 

The following result concerns the dependence of the concave envelope from the interval: 

Proposition 3.3. Let si < s 2 , f € C k ' ([0, S2]) and assume that ||/'||c°([o. s 2 l) — Ci^o- Then there are 
constants C 2 and C3 such that 

||conc[ , Sl ]/ - conc[ , S2 ]/||co([o, Sl ]) < C 2 S Q (s 2 - si) (3.11) 

and 

||(conc[ , Sl ]/ - conc[o !S2] fy\\ C o(io, Sl ]) < C 3 S (s 2 - si). (3.12) 
Proof. We will proceed in several steps: 
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1. For every r £ [0, si], 

/(t) < conc [0 , sl ]/(T) < conc [0 , S2 ]/(r) (3.13) 

Indeed, the restriction of concro jS2 ]/(r) to the interval [0, si] is a concave and non decreasing function 

which is also greater then /. From ()3.13p one deduces that, for every r £ [0, Si], the following 
implication holds: 

/(t) = conc [0 , S2 ]/(r) conc [0 , Si ]/(t) = conc [0 , S2 ]/(t). (3.14) 

2. It is not restrictive to assume 

3s<e]si,s 2 ]: /(s) > — (conc[ 0)S1 ]/) ^ (s - s x ) + conc [0 , si ]/(si) (3.15) 
Indeed, if (|3.15[) is not satisfied then 

f conc [0 , sl ]/(T) t £ [0, si] 

^ I ^~( conc [o,si]/) T _ g [s~ s i] + c° nc [o,si]/(si) r e]«i, sj] 

is a concave function which is greater or equal to / on [0, S2]. Thus, 

ft(r) > conc [0 , S2 ]/(r) V r £ [0, s 2 ]. 

In particular, on [0, s±\ 

conc [0iSl ]/(r) > conc [0i S2 j/(t) 
and thus by (|3.13p the two coincide. 

3. Before proving (|3.1ip we introduce other preliminary observations. 
Given two arbitrary points a and b in [0, S2], let £ the line 



Then 

concjo 
Indeed, 



£ ab (s) = -^(f(a)(b-s)+f(b)(s-a 



,si](f-£ab) = (conc [0: Sl ]f^j - lab conc [0j S2 j (f - £ a = (conc [0 ^ 2] /) -4b- (3-16) 



cone 



[0, S i](/-4feJ < (conc [0 -£ q6 
because ^conC[ 0jSl ]/^ — £ a b is a concave function bigger than f — £ a b- On the other side, 

(conC[o !Sl] /^ < l ah + conC[o !Sl] (/ - 4bJ 

since 4fc + conc[ 0i Sl ] (f — £ a tj is a concave function bigger then / — £ a b + £ ao = f. 

Equalities (|3. 16|) imply, in particular, that in proving (|3.11|) and (|3.12|) one can fix two points a and b 
in [0, S2] and assume that 

f(a) = fib) = 0. (3.17) 

Indeed, instead of / one can consider / — £ a b, which satisfies (|3.17p . Thanks to (13.16p . / satisfies p. lip 
and (|3.12p if and only if / — £ a ^ does. 
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4. Define 

s* := max{a G [0, a x ] : conc [0;Sl ]/(a) = conc [0j S2 ]/(s)}. 

If a < s*, then concro )S1 ]/(a) = concm, S2 ]/( s )- Indeed, by contradiction assume that there exists a 
r e [0, s*] such that 

conc [0)Sl] /(r) < conc [0!S2 ]/(r). (3.18) 

By ()3. 13[) this implies that f(r) < concjo, S2 ]f( T ) an d hence, thanks to the third part of Lemma T3.31 
there are r a and such that r a < t < Tfc and conC[ 0jS2 ]/ on [r a , Tfc] is a linear. Moreover, f(r a ) = 
conC[o !S2 ]/(r a ) and /(r 6 ) = conc[ 0i S2 ]f{n). Because of (|3.14p . this implies /(r ) = conc [0)Sl ]/(r a ). Also, 
because of the previous observation, it is not restrictive to assume f(r a ) = /(r^). Since concm ai i/ 
is a concave function, if concm, Sl ]f( T a) = and concro )S1 i/(r) < 0, then concro )S1 i/(t) < for every 
t > t. If Tb > si, then concro lS1 ]/ < concro, S2 ]/ on ]r a , si and hence s* < r a < t, which contradicts our 
hypothesis. On the other side, if Tb < Si then from f(rb) — conc[ 0i S2 ]f{ T b) an d (|3. 14)) we deduce f(rb) = 
conc[ 0j si]/( T b)- Being a concave function, concp. Sl ]f on [ T a, Tb] is above the line joining conc[ 0i Sl ]f( T a) 
and conc[ 0i Sl \f{ T b)- Thus, in particular, conc[ . Sl ]/(r) > concr 0j S2 i/(r), which contradicts (|3 . 1 8[) . 

5. We now prove p. lip . Let s* as in the previous step. If s s]s*, si], then concp, S2 ]f > f an d hence there 
are s a , Sb such that [s*, si] C [s a , Sb] C [0, S2] and on [s a , Sb] conC[ . S2 ]/ is the linear interpolation 
between f(s a ) and f(sb)- Relying on (|3. 16|) , one can assume f(s a ) — f(sb) — 0, getting conc[ 0i S2 i/ = 
on [s a , 8b]- Thus, conc[ . Sl ]f( s *) = an d conC[ . Sl ]/ < on ]s*, s\]. This implies 

(conc [0iSl ]/)'(s*) < 0. 

Since (conC[ 0iS2 ]/)' is the derivative of a concave function, then it is non increasing and 

(conc [0iSl ]/) (s) < Vsg[s*,si]. 

In particular, conC[o iSl ]/ is non increasing on [s*, s\] and hence its minimum on [s*, s\] is assumed in 
si. Thanks to the previous step, 

max {|conc [0)Sl i/-concr 0jSl ]/|} = max {|conc [0 , Sl ]/ - concr , si i/|}, 

which is equal to — concjo, Sl ]f( s i) by the previous observations. By the forth part of Lemma 13.11 
conqo, si]/( s i) — f( s i)- To conclude, one can observe that 

-f(si) = f(s b ) - f( Sl ) < dSisb - sx) < C 1 S(s 2 - sx). 

6. Before proving (|3.12p we introduce the following observations. Let g S C^fa, b] such that g(a) — 
g{b) = 0. Then 

g(s)>^(s-a)(s-b) Vs G [a, b]. (3.19) 

Indeed, consider for simplicity the case a = 0, b = 1. Let s denote a point at which g assumes its 
minimum and suppose, just to fix the ideas, that s < 1/2. Denote by h the first derivative of g: h 
is Lipschitz continuous function with Lipschitz constant smaller than k and such that h(s) = O.Thus, 
h(s) > k(s — s) if s < s. Moreover, Let us consider the Cauchy problem for g at s = 0: 

f g'(y) = h(y) > k{s - s) 
\ 5(0) = 0. 

By a comparison argument one gets 

g(s) > ks 2 /2 - kss Vs€[0, a]. (3.20) 

Moreover, 



VaG [0, 1], 3 (a) > <?(a) > - 




(3.21) 
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To complete the proof of (|3.19[) . let us consider to cases separately. If g'(l/2) < 0, consider the Cauchy 
problem 

h'(s)=g"(s)<k 
h(l/2) < 0. 

Then h(s) < k(s — 1/2) and thus the solution of 

g '(s) = h{s) < k(s - 1/2) 
. 3(1) -0. 

satisfies 

g(s)>ks(s-l)/2 VsE [1/2, 1]. (3.22) 

Define 

{ks 2 /2 — kss s < s 
-ks 2 /2 s<s< 1/2 

fes(s-l)/2 s>l/2. 

By direct check one gets ip{ s ) > ks(s — l)/2 for every s E [0, 1]. Putting together (|3.20p . (I3.2ip and 
(|3.22|) one gets therefore 

g (s) >ip(s) > ks(s- l)/2 Vse[0, 1]. 

If g'(l/2) > then 

g(s) > ks{s - l)/2 Vs€ [1/2,1]. (3.23) 
Indeed, thanks to (|3.2ip ,g(l/2) > — fe/8. Assume by contradiction that 

0(a) < fcs(s - l)/2 (3.24) 

for some s G [1/2, 1]. Define 

s n := max{s : g(s) > fcs(s — l)/2}. 
Since g(l/2) > — k/8 and g(s n ) = ks n (s n — l)/2, then there exits r G [1/2, s n ] such that 

.g'(r) < fc(r - 1/2). 

Consider the Cauchy problem 

h'(s) = g"(s) < k 
h(r) = ff'(r) < k(r - 1/2), 

then /i(s) < fc(s — 1/2) for s > t. Then the solution of the backward Cauchy problem 

g '(s) = h{s) < k(s - 1/2) 
3(1) = 

satisfies g(s) > ks(s — l)/2 for every s E [r, 1]. This contradicts (|3.24p and hence (|3.23p holds. One 
can then exploit the same function -0 considered in the case g'(l/2) < and conclude. 

7. We now prove p,12|) . Let s* be as in the previous steps, then conC[o jS2 ]/ = conC[o iSl ]/ on [0, s*} and 
hence 

(conc [0 , Sl ]/)'(s) = (conc [0 , S2] /) (s) s E [0, s*[. 

Thus, 

||(conc [0 , si ]/)' - (conc [Q aa] /)'|| C o [0 , Sl ] = max |(conc [(Ml] /)'(s) - (conc [0;S2] /)'(s)|. (3.25) 
As in step 4, denote by s a and Sb the values such that 

(conc [0 , S2 ]/j (s) = s e\s ai s b [. 
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As before it is not restrictive to assume f(s a ) — f(sb) — and in this case 

(conc [0jSl ]/) (s) < s S [s*, si] 
Since it is a non increasing function, the maximum in (I3.25|) is given by — (conqo, S2 ]/)'( s i)- Define 

s° := sup{s S [0, si[: conc [0)Sl] /(s) = f(s)}. 
If s° = 0, then for every s 6 [0, si[ 

f(s) < conc [0 , si ]/(s) < conc [0 , S2 ]/(s). 
In particular, one has s a — and hence /(0) = 0. By the third part of Lemma 13.31 

(conc [0 , Sl ]/j (si = = ■ 

Now we can apply the previous step to / defined on [0, Sb], thus obtaining 

j~ (s~^J k . . . . A^, . 

< ^— si(s b - Si) < -{s b - Si) < -(s 2 - Si). 

si zsi 2 2 

If s° = si, then (conc[o )Sl ]/)'(si) = f'(si). Consider two subcases separately. If Sb < s 2 , then also 
(concjo, S2 ]/)'(s 2 ) = f'{s2) and hence /'(s 2 ) = because we are assuming f(s a ) = f(sb) — 0. Then 

-/'(si) = /V) - f'(si) < Hs b - Sl ) < k{s 2 - si). 

On the other side, if s& = s 2 one can proceed as follows. Since f'(si) < 0, /(si) < and /(s 2 ) = 0, 
then there exists s G [si, s 2 ] such that f'(s) — 0. Then 

-/'(si) = /'(§) - < fc(3 - si) < fc(s 2 - si). 

We are now left to deal with the case < s° < s\. One then has 

- (conc [0 , Sl] /)'(si) = f{S °l~J s l Sl) (3.26) 

Because of the implication (|3.14[) , s a < s° and hence /(s°) < 0. Moreover, /(s°) = conc[ Sl ]/(s°), 
/(s 1 ) = conc [0 , Sl ]/(si) and (conc [0;Sl] /) (s) < on [s°, s x ]. Thus, 

/(s 1 ) < < 0. 

Since /(s&) = 0, there exists s 6 [si, Sb] such that /(s) = f(s°). Applying the previous step to the 
function f — f(s°) defined on the interval [s°, s], one gets /(s 1 ) — /(s°) > k(si — s°)(si — s)/2. Inserting 
in (13.26p . one gets 

/ k k 

-(conc [0 , si] /) (si) < -(§ - si) < -(s 2 - si). 

This concludes the proof of (|3.12|) and hence of the lemma. 

□ 



Combining Lemma l3.4l and Proposition [231 one finally gets 

Proposition 3.4. Let s± < s 2 , / £ ^^([0, si]) and assume that ||/'||c°([0,s2]) — Ci^O- Then there are 
constants C 2 and C3 such that 

||mon [0 , Sl ]/ - mon[o, S2 ]/|| c o([o, sl ]) < C 2 <5 (s 2 - s x ) 

and 

||(mon[ 0jS1 ]/ - mon [0 , S2 ]/)'||co([o, Sl ]) < C 3 5o(s 2 - si). 
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3.2 The hyperbolic limit in the case of an invertible viscosity matrix 
3.2.1 The hyperbolic limit in the case of a Cauchy problem 

In this section for completeness we will give a quick review of the construction of the Riemann solver for the 
Cauchy problem. We refer to [7] for the complete analysis. The goal is to characterize the limit of 

E{u s H + A(u E , u%)u% = eB{u e )u% x 

u~ x < (3.27) 



u e (0, ,x) =- , 

[ uq x >0 

under Hypotheses [TJ OH ID and [H 

The construction works as follows. Consider a travelling profile 

B(U)U" = (A(U, U') - aE(U)\U', 

then U solves 

u = v 

B(u)v' = (A(u, v) - aE{u)\v (3.28) 
a 1 = 0. 

Let Xi(uo) be the i-th eigenvalue of E~ 1 (uq)A(uo, 0) and let r,(«o) be the corresponding eigenvalue. If one 
linearizes the previous system around the equilibrium point (uq, 0, Xi(uo)) one finds 

I N 

A(u , 0) -Xi(u )E(u ) 


The generalized eigenspace corresponding to the eigenvector is 

V c = Uu, xn(uo), cr) : u £ R N , x, a e R~\. 

It is then possible to define a center manifold which is parameterized by V c : we refer to |33j for an extensive 
analysis, here we will just recall some of the fundamental properties of a center manifold. 

Every center manifold is invariant with respect to (|3.28|) and moreover satisfies the following property: 
let (u°, p°, a ) belong to M c and denote by (u(x), p(x), a(x)) the orbit starting at (u°, p°, a ). Then 

lim (u(x), p(x), cr{x))e~ cx/2 = (0, 0, 0) lim f u(x), p(x), tr(xj) e- cx ' 2 = (0, 0, 0). 

The constant c is the separation speed defined in Hypothesis [31 A center manifold A4 C is defined in a 
neighbourhood of the equilibrium point (fii, 0, A(uj) and it is tangent to V c at (itj, 0, X(iii)). 
Also, the following holds. Define 



" := j(it , y^ y %jrj(u ), a) :u£ R N , Xj, a GR|, 



ysu ._ 

where ri(ito) . . .rj_i(uo), rj+i(uo) . . .rjv(uo) are eigenvectors of E~ 1 (uo)A(uq, 0) different than rj(uo). If 
we write R 2N+1 = V c ® then the map 

$ c : V^ c -> T/ c © V us 

that parameterizes A^ c can be chosen in such a way that if 7r c is the projection from V c © V us onto then 
7r c o $ c is the identity. 

Fix a center manifold A4 C . Putting all the previous considerations together, one gets that a point (it, v, a) 
belongs A4 C if and only if 

v = ViViiuo) + ^ipj(u, Vi, ai)rj(u ). (3.29) 
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Since all the equilibrium points (u, 0, <7j) lay on the center manifold, it turns out that when Vi = 0, then 

ipj(u, 0, o-j) = Vj, u, a t . 
Hence for all j, ^>j{u, Vi, Cj) = Vi<f>j(u, v?, ai) for a suitable regular function 0j. The relation (|3.29p can be 



therefore rewritten as 



i(n(uo) + y^ y (j)j(u, Vi, (Ti)rj(u )j :=Uifi(u, t!j, crj). 



Because of the tangency condition of the center manifold to the center space, it holds 

h(ua, 0, Xi(u ^j = ri(u ). 
Inserting the expression found for v into (13.281) . one gets 

v ix B(u)f l + vlB{u)Df t fi + ViV ix B(u)f iv = vJa{u, Vifi) - (TiE(u)\fi. 
Considering the scalar product of the previous expression with fi one obtains 

VixUh, B{u)fi) +Vi(h, B(u)f iv )] =Vi({fi, [a(u, Vifi) -aiE{u)\fi) -i»i(fj, B(u)Dfifi) 



Hence setting 



Ci(u, Vi, crj) := (f i; B{u)fi) +Vi{fi, B{u)f iv ) 

di(u, Vi, crj) := (f^ (^A(u, Vifi) ~ <?iE(u)jfi) - Vi{h, B(u)Dfifi) 



one can define 



4>i(u, Vi, <Ji) := 



a,i(u, Vi, Oi) 



Ci(u, Vi, (Ti) ' 

The fraction is well defined since Ci(uo, 0, Ai(uo)) > cb(uq) > and hence Ci is strictly positive in a small 
neighborhood. The constant cb is as in the third condition in Hypothesis [T] 
Thus, 

a, <3c* , 1 ddi 



Since 
then 



ai(ua, 0, Ai(tto)) > c B (u ) = 



do. 



(uo,0,Ai({io)) Cl 



7 (a(u , 0) - A 4 S(w )^ri(u )) - (r-i(uo), £?(u )ri(tto)> 



(ri(wo), I A(uo) 0) - XiE(u ) )f it7 ) 



(3.30) 



-ce(uo) < 0. 



In the previous computations, we have exploited the symmetry of A(uq, 0) and E(uq) and hence the fact 
that 

(ri(u , 0), (A(u , 0) - X.Eiuo)^^) = ((a(Uq, 0) - XiE(u )^n(uo), ha) = 0. 
Also, the constant ce in (|3.30|) is the same as in Hypothesis [T] 
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In conclusion, system (|3.28[) restricted to A4 C can be rewritten as 



U 1 = Vifi(u, Vi, Oi) 

v[ = (f>i(u, Vi, a t )vi (3.31) 
at = 0. 



One actually studies the following fixed point problem, denned on a interval [0, s,]: 

u(t) = u + ( fi(u(0, Vi(0, o-f(0)^ 



We have used the following notations: 



Vi(r) = /i(r, u, Vi, ai) - conc/i(r, u, u,, cr*) (3.32) 



where 



Xi[ui, Vi, fTi](0 = <^i^«i(0) °*(£)J +cb(uo)ct. 

Also, conc/i denotes the concave envelope of the function 

conc/j(r) = inf{/i(s) : h is concave, h(y) > fi(y) Vy £ [0, Sj]}. 

The link between (13.32[) and (|3.28p is the following: let (it,, Vi, Ui satisfy (|3.32|) . Assume that t>, < on }a, b[ 
and that Vk{a) = Vk{b) = 0. Define a^(r) as the solution of the Cauchy problem 

den 1 



dr v iij) 
a(a + b/2) = 

then (v,i o ckj, o q„ <Ji o ctj) is a solution to (|3 . 3 1 1) satisfying 

lim Ui o Q!i (a;) = Ui(a) lim u, o 0^(2;) = Ui(b). 



x — > — oo .r — >+oo 



Thus, iii(a) and itj (6) are connected by a travelling wave profile. 

As shown in [7], (|3.32|) admits a unique continuous solution (ui, Vi, <Ti). Also, one can show that Ufc(0) 
and Uk(si) are connected by a sequence of rarefaction and travelling waves with speed close to Aj(uo). If 
u(t, x) is the limit of 

" E(u e )ul + A(u s , u%)u% = eB{u e )u% x 
uq x < 
Ui(si) x > 



w £ (0, ,x) 



then as e — > + 

!M0 X < <7j(0)i 

Ui(s) x = ai(s)t (3.33) 

Ui(Si) X > (Ji{si)t 

The i-th curve of admissible states is defined setting 

T s 4 .u := Ui(si). 

If Si < 0, one considers a fixed problem like (|3.32|) , but instead of the concave envelope of fi one takes 
the convex envelope: 

conv/i(r) = sup{h(s) : h is convex, h(y) < f l (y) Vy}. 
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Again, one can prove the existence of a unique fixed point (itj, Vi, o~i). 

Also, in [7] it is proved that the curve T*. is Lipschitz continuos with respect to Si and it is differentiable 
at Si = with derivative given by t^uq). Moreover, the function Tj. is Lipschitz continuos with respect to 

Uq. 

Consider the composite function 

%/}(uo, s%...s N ) = T s \ o . . . T^ N u 

With the previous expression we mean that the starting point for T^ -1 is not uo but T£\uq. Thanks to 
the previous steps, the map tp is Lipschitz continous with respect to s\ . . . Sn and differentiable at si = 0, 
. . . sn = 0. The column of the jacobian are ri(uo) • ■ ■ r/v(wo)- Thus the jacobian is invertible and hence, 
exploiting the extension of the implicit function theorem discussed in [24] (page 253), the map </>(uo, •) is 
invertible in a neighbourhood of (s± . . . sn) — (0 ... 0). In other words, if u~ is fixed and is sufficiently close 
to uq, then the values of si ... sn are uniquely determined by the equation 

u~ = ip(u , si . . . s N ). 

Taking the same u~ as in (|3.27| . one obtains the parameters (s\ . . . sn) which can be used to reconstruct 
the hyperbolic limit u of (|3.27|) . Indeed, once (si . . . sjv) are known then u can be obtained gluing together 
solutions like (|3.33p . 

3.2.2 The hyperbolic limit in the non characteristic case 

The goal of this section is to provide a characterization of the limit of the parabolic approximation (|3.ip in 
the case of a non characteristic boundary, i.e. when none of the eigenvalues of E~ 1 (u)A(u, u x ) can attain 
the value 0. More precisely, we assume the following. 

Hypothesis 8. Let Ai (u) . . . Ajv (u) be the eigenvalues of the matrix E~ 1 (u)A(u, 0). Then there exists a 
constant c such that for every u 

Ai(«) < • • • < An(«) < ~\ < < I < A„ +1 (u) < • • • < X N (u). (3.34) 

Thus, in the following n will denote the number of eigenvalues with strictly negative negative real part 
and N — p the number of eigenvalues with strictly positive real part. 

To give a characterization of the limit of (|3.1[) we will proceed as follows. We wil construct a map 
4>(uq, si ... sn) which is a boundary Riemann solver in the sense that as (si . . . sn) vary, describes states 
that can be connected to uq. We will the show that the map <fi is locally invertible. Hence, given uq and Ub 
sufficiently close, the values of (si . . . sn) are uniquely determined by the equation 

Ub = <f>(u , si . .. s N ). 

Once (si ... sat) are known the limit of (|3.1[) is completely characterized. The construction of the map <fi is 
divided in some steps: 

1. Waves with positive speed 

Consider the Cauchy datum u$, fix (N — k) parameters (sjv-fc . . . sn) and consider the value 

U = 1 „„ , O . . . 1 Uq. 

The curves T Ar ~ fc . . . are, as in Section l3\2.1l the curves of admissible states introduced in 7 . The 
state uo is then connected to u by a sequence of rarefaction and travelling waves with positive speed. 

2. Boundary layers 

We have now to characterize the set of values u such that the following problem admits a solution: 

A(U,U X )U X = B{U)U XX 
U(0) = u 

lim x ^ +00 U(x) = u. 
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We have thus to study system 

U x =p 



Px = B{U)- X A{U, p)p (3 - 35) 

Consider the equilibrium point {u, 0), linearize at that point and denote by V s the stable space, i.e. 
the eigenspace associated to the eigenvalues with strictly negative real part. Thanks to Lemma [3.11 
the dimension of V s is equal to the number of negative eigenvalues of E (u)A(u, 0), i.e. to n. Also, 
V s is given by 

n n 

V s = + jrzXt(u), ^2xiXi(u)), X!...x n eM|, 

where fXi(u) . . . /i„(u) are the eigenvalues of B~ 1 (u)A((u), 0) with negative real part and )Ci(u) ■ ■ ■ Xn(u) 
are the corresponding eigenvectors. 

Denote by M s the stable manifold, which is parameterized by V s . Also, denote by <j> s a parameteriza- 
tion of M s : 



Let tt u be the projection 



<f> s : V s -» R*. 



(u, p) H-> u 

If u € ^u(4>s(si ■ ■ ■ s n )J for some si . . . s n , then system p.35[) admits a solution. Note that thanks 
to classical results about the stable manifold (see eg [33]) the map <j> 8 is differentiable and hence also 
7r u o <j) s is differentiable. In particular, the stable manifold is tangent at s\ = . . . s n = to V s and 
hence the columns of the jacobian of tt u o (j) s computed at s\ = . . . s n = are Xi ■ ■ ■ Xn{u). 

Note that the map ir u o </> s actually depends also on the point u and it does in a Lipschitz continuos 
way: 

\tt u o <f> a (ui, si... s n ) ~ir u o 4> s (u2, si . . . s„)\ < L\ui - u 2 \ 

3. Conclusion 

Define the map <fi as follows: 

4>{u , si . . . s n ) = ttu o 4>s(r^ N ~ k k o . ..T^ N u , si . . . s n J (3.36) 

From the previous steps it follows that <f> is Lipschitz continuos and that it is differentiable at s\ = 
... sn = 0. Also, the columns of the jacobian are Xi(""o) ■ ■ ■ Xn(uo), r n +i(uo) . . . rjv. Thus, thanks 
to Lemma l3.1[ the jacobian is invertible. One can thus exploit the extension of the implicit function 
theorem discussed in (24] (page 253) and conclude that the map 4>{uq, ■) is invertible in a neighbourhood 
of (si . . . sn) = (0 ... 0). In particular, if one takes Ub as in (|3.1|) and assumes that \uo — Ub\ is sufficiently 
small, then the values of si ... sn are uniquely determined by the equation 

u b = (j)(u , si . . .s n ) (3.37) 

Once the values of s\ ... sn are known, then the limit u(t, x) can be reconstructed. In particular, the 
trace of u on the axis x — is given by 

u:=T?+\o...T s N n u . (3.38) 

The self similar function u is represented in Figure 13.2.21 and can be obtained gluing together pieces 
like (13331) . 
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Theorem 3.1. Let Hypotheses]^ [TJ , and\^ hold. Then there exists S > small enough such that 
the following holds. If \uo — Uf,\ << 5, then the limit of the parabolic approximation (|3.1|) satisfies 

Ub = <t>(uo, si . . . s N ) 

for a suitable vector (s\ . . .sn)- The map <p * s defined by (|3.36p . Given uq and Ub, one can invert <f> and 
determine uniquely [s\ . . . sjv)- Once {s\ . . . sjv) are known one can determine a.e. (t, x) the value u(t, x) 
assumed by the limit function. In particular, the trace u of the hyperbolic limit in the axis x = is given by 
(13351) . 

3.2.3 The hyperbolic limit in the boundary characteristic case 

The aim of this section is to provide a characterization of the limit of the parabolic approximation (|3.1|) when 
the matrix B is invertible, but the boundary is characteristic, i.e. one eigenvalue of A(u, u x ) can attain the 
value 0. 

Let 5 be the bound on the total variation introduced in Hypothesis [H Moreover, as in the previous 
section we will denote by Xi(u) . . . Xn(u) the eigenvalues of the matrix i? _1 (u)A(u, 0). The eigenvalue 
Afe(u) is Lipschitz continuous with respect to u and hence there exists a suitable constant M such that, 
if | Afc(wo)| > MS, then |Afe(w)| > for every u in a neighbourhood of uq of size 5, \u — uq\ < S. Thus, 
in particular, if Afc(wo) > MS then Hypothesis [S] is satisfied for some c and hence the boundary is non 
characteristic. In this section we will therefore make the following assumption: 

Hypothesis 9. Let S the same constant as in Hypothesis^ and let M the constant introduced before. Then 

\\k{u k )\<MS. (3.39) 

Note that, because of strict hyperbolicity (Hypothesis [3]) the other eigenvalues of E~ 1 (u)A(u, 0) are 
well separated from zero, in the following sense: there exist a constant c > such that for all u satisfying 
\uo — u\ < S, it holds 

Xi(u) <•■■ < Xk-i(u) < -c < < c < X k+ i(u) < • • • < X N (u). (3.40) 

The notation used in this section is therefore the following: (fc — 1) denotes the number of strictly negative 
eigenvalues of i? _1 ^4(u, 0), while (N — k) is the number of strictly positive eigenvalues. 
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The study of the limit is more complicated in the boundary characteristic case than in the case considered 
in the previous section. The main ideas exploited in the analysis are described in the introduction, here 
instead we will give the technical details of the construction. As in the previous section, the characterization 
of the limit will be as follows: we will construct a map (j>(uo, s% . . . sn) which describes all the states that can 
be connected to Uo- We will then show that the map is locally invertible: thus, given uq and Ub sufficiently 
close, the values of (si . . . sn) are uniquely determined by the equation 

u b = <f>(u , si.. . s N ). 

Once (si . . . sn) are known one can proceed as in the previous section and determine a.e. (t, x) the value 
u(t, x) assumed by the limit function. The construction of the map cj> will be given in several steps: 

1. Study of the waves with uniformly positive speed 

Consider the Cauchy datum u$, fix (N — k) parameters (sj\r-fc ■ ■ • sjv) and consider the value 

rr>N — k rriN - 

u k=T SN _ k o...T SN u . 

The curves T„ JV_fc . . . are, as in Section 13.2.11 the curves of admissible states introduced in [7]. 
The state uq is then connected to Uk by a sequence of rarefaction and travelling waves with uniformly 
positive speed. The speed is uniformly positive in the sense of (|3.40p . 

2. Reduction on the center stable manifold 
Consider system 

u x = V 

Px = B{u)- 1 (A{u 1 p)~aE{u)) P (3.41) 
o x =0 

and the equilibrium point (jik, 0, X(uk))- Linearizing the system around (uk, 0, A(i2fc)) one obtains 

/ 1% \ 

B- 1 (u k )(A(u k ,0)-X k (u k E(u k )) (3.42) 
\ '0 / 

Thanks to (|3.40p , the matrix E (v,k) — \k(uk) A(iik , 0) has 1 null eigenvalue and (N — k) strictly negative 
eigenvalues. Because of Lemma 13.11 one can then conclude that the matrix (|3.42p has N + 1 null 
eigenvalues and (N — k) eigenvalues s with strictly negative real part.. 

Let V c be the kernel of (|3.42|) . V s the eigenvspace associated eigenvalues s with strictly negative real 
part, and V u the eigenspace associated to eigenvalues s with strictly positive real part. There exists a 
so called center stable manifold Ai cs which is parameterized by V S @V C : we refer to [33] for an extensive 
analysis, here we will just recall some of the fundamental properties of a center stable manifold. If we 
write R 2N+1 = V c © V s © V u , then the map 

$ cs : V s © V c -» V c © V s © V u 

that parameterizes M cs can be chosen in such a way that if tt cs is the projection from 
V c © V s © V u onto V c © V s , then n cs o $ cs is the identity. 

Every center stable manifold is invariant with respect to (|3.4ip and moreover satisfies the follow- 
ing property: let (u°, p°, <j°) belong to M cs and denote by (u(x), p(x), <j{x) the orbit starting at 
(w°, p°, cr°). Then 

lim (u(x), p(x), a(x)e- cx/2 ) = (0, 0, 0). 

x — >-{-oo V / 

The constant c is the separation speed defined in (|3.40p . We are also assuming that \k(uk) is bigger 
than — c/4: this is a consequence of Hypothesis O provided that S is sufficiently small. A center stable 



43 



manifold A4 CS is defined in a neighbourhood of the equilibrium point (uk, 0, \(u k ) and it is tangent to 
V C (BV S at (u fc , 0, A(u fc ). 

In the following, we will fix a center stable manifold Ai cs and we will focus on the solutions of (|3.4ip 
that lay on A4 CS . The center and the stable space of (|3.42[) are given respectively by 



V c = {(«,«*?*(«*), a): ue R N , v k s , a G K.| 



and 



{{U k + J2 —K(u k ),Y, vt csMu k ), Afc(«fc))}. 



Xl---Xfe-1 



In the previous expression, r k (uk) is a unit vector in the kernel of ( A(Uk, 0) — \ k {uk)j , while 

are the eigenvectors of B~ 1 (uk)(^A(uk, 0) — \k(u k )j associated to the eigenvalues [i\ . . .fJ>k-i s with 
strictly negative real part. Since <I> CS o7r cs is the identity, then a point (u, p, a) belongs to the manifold 
Ai cs is and only if 

P = yl csXt(uk) + v k cs r k + Y <j>L(u, v\ s . . . v k s , a)Xi(u k ), 

i<k i>k 

where Xk+i{uk) ■ ■ - XN(uk) are the eigenvectors of \A(uk, 0) — Xk(u k )j associated to eigen- 

values with strictly positive real part.. Since all the equilibrium points (u, 0, a) belong to the center 
stable manifold, then we must have p = when v^ s . . . v k s are all 0. Thus, since Xfe+i("fc) ■ ■ • XN(uk) 
are linearly independent, 4> l cs {u, 0, er) = for alii = k + 1 . . . N and for every u and a. This means, in 
particular, that for every i = k + 1 . . . N there exists a /c-dimensional vector 4> l cs (u, v]. s . . .p k s , cr) such 
that 

<t>cs( U , V ls--- V cs) = (<t>cs{ U , «4 • • • V cs: V): «L • • • V cs) ■ 

Define 



R, 



,(u, v\ s . ..v h cs , a) := (xi| ... \xk-i\rkj + (xk+i\- ■ ■ \xn) 



\ ^ s (u,vi...v k cs ,a) ) 



then a point (u, p, a) belongs to Ai cs if and only if 

p = R cs (u, V cs , a)V cs , (3.43) 
where V cs denotes (vl s . . . v k s ) T . Since M cs is tangent to V c © V s at (Uk, 0, X k (uk)), then 

R cs {u k , 0, \ k {u k )) = Ixfe-iln:)- 

Plugging p = R C s(u, V cs , &)V CS into system (I3.64p one gets 



((D U R CS )R CS V CS + (i? w ii cs )V^ sa: ^ V cs + RcsV£ s — B 
Moreover, the matrix 

D(u, V cs ) = (da,) d lh = J2 ^^Vj 

satisfies 

!>(«, 0) = ((A,^)^)^ = !>(«, F CS )K CS . 



(3.44) 



44 



Hence (|3.44p can be rewritten as 

\R CS +D(u, vSjVc* = \ B- 1 (A-ctE) -{D u R cs )R cs V cs )R cs V cs . (3.45) 

Since 

^R cs + D{u, V cs )^j (u k , 0, Afc(ufc)) = (xA ■ ■■ \Xk-i\rk), 

then the columns of R cs + D{u, V cs ) are linearly independent in a small enough neighbourhood. Thus 
one can find a matrix L(u, V cs , a) such that 



L(u, V CS1 a)[R cs + D(u, V cs ) 
Multiplying (13.45j) by L(u, V cs , a) one finds 

V x cs = \LB-\A - <jE)R cs - L(D U R CS )R CS V CS )V CS . (3.46) 



Define 

A cs (u, V cs , a) := LB~ X {A - aE)R cs - L{D U R CS )R CS V CS , (3.47) 

then by construction A cs (uk, 0, \k(uk)) is a (k x k) diagonal matrix with one eigenvalue equal to 
and N — k eigenvalues with strictly negative real part. Also, if Ik = (0 . . . 0, 1) and Sf. = then 



>J dK 



da 



1 4 = -k < (3.48) 



-Uk, V C 8=0, <j=X k {uk)' 

for a suitable strictly positive constant k. This is a consequence of conditions 1 and 3 in Hypothesis [T] 
In conclusion, the solutions of system (|3.4ip laying on A4 CS satisfy 

u x = R cs (u, V cs , o)V cs 

V csx - A cs (u, V cs , a)V cs (3.49) 
&x = 

where A cs (u, V cs , a) is defined by (|3.47|) . 

3. Analysis of the uniformly stable component of p.49p 

Linearizing system (|3.49l) around the equilibrium point [Uk, 0, \k{uk)) one obtains the matrix 

R cs (u k , 0, A fc (tifc)) \ 

A cs (u k , 0, A fc (w fe )) , (3.50) 
/ 

which has fc— 1 distinct eigenvalues with strictly negative real part and the eigenvalue with multiplicity 
N + 1. Also, the manifold 

E := | (it, 0, er), u G R N , a G m| 
is entirely constituted by equilibria. 

Then there exists a uniformly stable manifold M.™ which is invariant for (I3.49| and is characterized 
by the following property: given V c ° s , a ) belonging to M.™ , denote by (u(x), V cs (x), a(x)) the 
solution of (|3.49|) starting at (u°, V° s , a ). Then there exists a point (it 00 , 0, cr°°) belonging to E such 
that 



lim 

cc — >+oo 



(\u(x) - u°°\ + \V cs (x)\ + \a(x) - a°°\)e cx / 2 = 0. 
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We remark that the uniformly stable manifold Ai^ is contained but in general does not coincide with a 
given center stable manifold. Indeed, M.™ does not contain, for example, the trajectories that converge 
to a point in E with speed slower than e~ cx . 

The existence of the uniformly stable manifold is implied by Hadamard-Perron theorem, which is dis- 
cussed for example in [33 . The application of the uniformly stable manifold to the study of hyperbolic 
initial boundary values problems is derived from [3] . In our case, the existence of the uniformly stable 
manifold is loosely speaking guaranteed by the strict hyperbolicity of A (Hypothesis [3|) , namely the 
fact that the eigenvalues of A are uniformly separated one from each other. In the following, we will 
just recall some of the properties of the uniformly stable manifold. 

Let 



V s = 



{(u k +j2 v ixi(uk), J2 v &' v i--- v s 1 e R }> 



i=l 



where \i are as before the eigenvectors of B 1 (A — E) and e*j 6 R fc are the vectors of the canonical 
basis of M fc , ei = (1, . . . 0) and so on. The uniformly stable manifold MJf is parameterized by E © V s 
and hence has dimension N + k — 1. The parameterization 

<S> US : V s © E -> M fc+2 

can be chosen in such a way that the following property is satisfied. Let 
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Uu k +v k r k (u k ), v k e kl \ k {u k j), v k S r\. 



If we write R 7V+fe+1 as V c © V s © E and we denote by -k us the projection onto V s © E, then we can 
have n us o <fr us equal to the identity. 

By considerations analogous, but easier, to those performed to get (|3.43p one deduces that a point 
(u, V cs , a) belongs to Ai^ s if and only if 

V cs =R s (u, V„ a)V s , 

where V s — (w^ s . . . v k J x ) are the components of V cs along e\...e k ~\. The matrix R s belongs to 
M fex(fe-i) and 

R s (u k , 0, X k {uk)^ = (ei\ ■ ■ ■ \Sk-ij. 
We want now to compute A CS R S V S . Plugging the relation V sx = A CS V S into p.49p one gets 



R s + A 



A CS R S — D U R S R CS R S V S 



v.. 



(3.51) 



where D s is a matrix such that D s (u, 0, a) for every u and a. Thus, 

V sx = A S V S 

for a suitable matrix A s such that A s ^u k , 0, A/ c (u/ £ ^ is a diagonal matrix with all the eigenvalues with 
strctly negative real part. Plugging back into (|3.5ip one gets 



A CS R S V S 



Rs + D t 



A s + D U R S R CS R S V S W s . 



(3.52) 



4. Analysis of the center component of (|3.49p 

Linearizing system (|3.49|) around the equilibrium point (u k , 0, \ k (u k )) one obtains (|3.50[) and hence 
the center space is given by 



yc 



I (it, v k e k , a): u e l", v k , a G ffii. 
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Consider the center manifold M c , parameterized by V°: thanks to considerations similar, but easier, 
than those performed in the previous steps one gets that a point (u, V cs , a) belongs to M c if and only 
if 

V cs = f k (u, V k , <j)v kl 

where f k 6 R k and f k (u k , 0, \ k {u k )) = e k . 

Again, with considerations analogous, but easier than those performed at the previous step one gets 



A cs f k v k = 



(fk + f kv v k )(j) k + D u f k R cs f k v k 



v k 



(3.53) 



for a suitable function 4> k satisfying (j) k {u k , 0, X k (u k )j — 0. Also, thanks to (|3.48p . 



da 



= -k < 0. 



(3.54) 



u=u k ,V cs =0,<j=X k (u k 



5. Decomposition of (|3.49p in center and uniformly stable component 

In this step, we will fix a trajectory (u, V cs , tr) of (|3.49[) and we will decompose it in a uniformly stable 
and in a center component, in the following sense. We exploits the manifolds Ai c and introduced 
in the previous steps and we decompose 



V cs = R s (u, V a , o)V s + fk(u, v k , a)v k . 
Plugging this expression into the first line of (|3.49[) one gets 

u x = R CS R S V S + R cs f k v k . 
From the second line of (|3.49[) one gets 



(3.55) 



r k v kx 



D U R S (R CS R S V S + R cs f k v k ) + D V R S V S , 
D u f k (R cs R s V s + R cs f k v k ) + 



A CS R S V S + A 

cs^k^k ■ 



One can prove that 



D V R S V SX V S = D 8 (u, V a , <j)V s , 



for a suitable matrix D s which satisfies D(u, 0, it) = for every u and a (the same as in the previous 
step). Thus, exploiting (|3.52[) and (|3.53p . we get 



Rs + D t 



V s ., 



fk + fkvVk 



Vkx = A. CS R S V S + A cs f k Vk 

- D U R S (R CS R S V S + R cs f k Vk)Vs - D u f k (R cs R s V s + i? cs r fe w fc )w fe 



Rs + D s 



A s + D u RsR cs RsV s \Vs 



(fk + fkvV k )4> k + D u f k R cs f k v k 



D u Rs(R cs R s Vs + R cs fkV k )V s - D u f k (R cs R s Vs + R cs fkV k )v k 



Rs + Ds 



A S V S + (fk + fkvv k )4> k v k 



D u R s R cs fk + D u f k R cs Rs 



V s v k 



(3.56) 
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The functions R s + D s and f k + f k v k satisfy 

(R s + D s )(u k , 0, Afe(u fc )) = fe x | . . . |e fc -i 



(rfe +r fc w fe )^M fe , 0, A fe (ufc)^) = e fc 



and hence in a neighbourhood of (it^, 0, \ k (u k )) the columns of i? s + D s and 4 + f k v k are all linearly 
independent. Denote by l\ . . . i k the vectors of the dual basis and define 



L s (u, V s ,v k , u) 



Then multiplying (|3.56[) on the left by L s one gets 

V sx = A s (u, v k , V s ,a)V s 

where A s g M fc_1><fe_:L is given by 



A s (u, v k ,V s ,(j) = L & 



R s + D s 



A, - L, 



D u R s R cs f k + D u f k R cs R 6 



The matrix A s (u k , 0, 0, X k (u k )) is diagonal and all the eigenvalues have strictly negative real part. 
Also, multiplying (|3.56p on the left by £ k one gets 

Vkx — 4>kVk, 

where the real valued function <fi k (u, v kl V s , a) is given by 



D u R s R cs f k + D u f k R cs R s 



4>k{u, v k , V s ,(x) = l k {f k + f kv v k )4> k - 4 

Also, cj) k {u k , 0, 0, A fe (ufe)) = 0. 

In conclusion, system (|3.49|) can be decomposed as 

u x = R cs (u, v k f k + R S V S , <j)f k v k + R cs (u, v k f k + R S V 3 , a)R s V s 
v kx = (f>k(u, v k , V s , a)v k 
V sx = A s (u, Vk, V s , <j)V s 
a x =0 

Note that, thanks to (1X54]) . 

d<p k (u k , v k , V s , a k ) 



V 



(3.57) 



da 



-k < 0. 



(3.58) 



u=u k ,v k =0, V s =0, cr k =\ k (u k ) 



To study system (|3.57[) we will first consider the solution of (|3.57[) when u s = and V s = 0. Fix s k > 0: 
we will actually study the following fixed point problem, which is defined on the interval [0, s k ]: 



In the previous expression, 



u fc (r) = u k + r k {u k (0, MO, 0, o*(0)d£ 
Jo 

Vk(r) = f k [u k , v k , crfe](r) - mon [QtSk] f k [u k , v k , a k ]{r) 
, , 1 d 

m on[ 0iSfc ]/fe(T, Uk, v k , o k ). 



(3.59) 



; dr 



fk[uk, v k , cfc](r) = / X k [u k , v k , <J k ](£,)d(, 



(3.60) 
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where 

X k [u k , v k , <Jk}(0 = 4>k(u k (£), Vk(0> 0, ^fc(O) ( 3 - 61 ) 
The constant k is defined by (|3.58[) . 



If Sfc < 0, one considers a fixed problem like (|3.59|) but instead of the monotone concave envelope of fi 
one takes the monotone convex envelope: 

monconv/i (t) = sup{/i(s) : h is convex and monotone non decreasing, h(y) < fi{y) Vy s [0, Sk]}- 

In the following we will consider only the case Sk > 0, the case Sk < being entirely similar. 

The link between system (|3.59[) and system (|3 . 57[) in the case u s = 0, V s = is the following. Let 
(itfc, Vk, (Jk) solve (|3.59[) . assume that Vk(sk) < and define 

a(r) = - —-ds. (3.62) 

Jt V k (S) 

Let 

a := max-jr e [0, s fc ] : v k (r) = 0}, 
in the following we will prove that a(r) < +oo if and only if r > s. 

The function (ufcoa, Vk°a, 0) solves p.57p in the case u s = 0, V s = 0. If Ufc(sfc) = then (iik(sk), 0, 0) 
is also a trivial solution of (|3.57|) in the case u s = 0, V s = and the following properties are satisfied. 
Note that the function u k ° a is a steady solution of the original parabolic equation since a = 0: 

^4(wfe, u kx )u kx = B(u k )u kxx . 

Moreover, 

lim (u fe o a) (a;) = Mfe(s). (3.63) 

x — *oo 

From system p.59p we also get that u k is connected to u k (s) by a sequence of rarefaction waves and 
shocks with non negative speed. 

After finding a solution of (|3.59[) in the case u s = 0, V s = we will also find a solution of (|3.59p in the 
case u k = ilk, v k = and a = 0. We will also impose 

lim tt s (a;) = 

x — *oo 

and hence we will study the fixed point problem 

u s {x) = - [ R s (uk +u s (y), 0, V s (y))v s (y)dy 

Jx r* r / nt (3 - 64) 

V s (x) = e Ax V s (0) + e A ^[A s {u s (y)+u k ,0,V s (y))-A\V s (y)dy, 
where 

A = A s (ii k , 0, o). 

Note that again, being a — 0, u s provides a steady solution of the origin parabolic equation, 

A{lL s , U sx )u sx — B^Us^Usxj, 

Finally, we will consider a component of perturbation, due to interaction between the purely center 
component, defined by (|3.59[) and the purely stable component, which satisfies (|3.64[) . More precisely, 
we will define (U, q, p) in such a way that u = U + u k oa + u s , v k +q and V s +p is a solution of (|3.57p . 

6. Analysis of the purely center component 

The purely center component is the solution of system (|3.59p . 
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Lemma 3.5. Fix 8 > such that s k < 8 << 1. Then system (|3.59|) admits a unique solution 
(«fc, v/c, <7fc) satisfying 

(a) Uk G C ([0, Sfe]) 7 1 1 life — Ufc||co < c u 5, itfc is Lipschitz with constant Lip(uk) < 2. 

f&J «fc G C ([0, Sk)), ||«fc||c° < c jA w * s Lipschitz with constant Lip(vk) < c„<5. 

(cj cr fe G C°([0, Sfe]), ||crfe - Ajfe(«jfc)||co < c a S/r) k , a k is Lipschitz with constant Lip(a k ) < c a . 

The constants c u , c v , c v , c a) r\ k and c a do not depend on 8. The function f defined by (|3.60|) belongs 
to C^ ,:L ([0, Sk]) for a suitable constant k which again does not depend on 8. 

The symbol C^ ,:L ([0, s k }) denotes the space of functions / G C 1 ([0, s k ]) such that /' is Lipschitz con- 
tinuous and has Lipschitz constant Lip{f) < k. 

Proof. The proof of the lemma relies on a fixed point argument. Only the fundamental steps are 
sketched here. 

Define 

X uk ■= {u k G C°([0, Sfc], R N ) : \\u k - u k \\ c o < c u 8, u k is Lipschtiz, Lip(u k ) < 2}, 

X vk := {u k GC°([0, s fe ], R) : ||v fc || C o < 8, v k is Lipschtiz, Lip(v k ) < c v S}, (3.65) 

X„ k := {cr fc 6 C°([0, s fe ], R) : ||er fc - X k {uk)\\a < c a 8 2 , cr fe is Lipschtiz, Lip(a k ) <c a ). 

In the previous expression, || • || CT is defined as follows: 

H| CT :=^IH| C o, 

for a suitable constant r\ k whose exact value does not depend on 8 and will be determined in the 
following. Also the constants c u , c v , c v , c a and c CT do not depend on 8 and their exact value will be 
estimated in the following. 

If (tife, Vk, o~k) G X uk x X vk x X ak then the function / defined by (|3.60[) satisfies / G C k ' 1 for a large 
enough constant k. Moreover, exploiting (|3.39p and (|3.6ip . one gets 

WfWco < s\\Xk\\co 

4>k(uk, o, x k (u k )j + 4> k (u k , v k , o-^j - <j>k(uk, o, \k{u k )j co + X k (u k ) 



<8 



+ ||Afc(ufc) - o- k \\ c o 



< 8 ( + 0{l)(c u 8 + c v S + c a S) + MS + —8 

\ ^ , 

< ^c v 8 

if 8 is sufficiently small and c v sufficiently large. 
The same computations ensure that 



(3.66) 



2 

for a large enough c v . Finally, 



ll/'llco < bv8 (3.67) 



\f'{ x ) ~ f'(y)\ = \X k (u k (x), v k (x), a k {x)^j -X k (u k (y), v k (y), a k (y)j\ 
< 0(1)(2 + ~c v 8 + cj)\x -y\< k\x - y\ 



(3.68) 
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for a sufficiently large constant k. In the previous estimate, we have exploited the following observation: 
since by (gUgj) 

u k ,v k =0, <J k =\ k (u k ) 



da 
then 



< 0(1)6. (3.69) 

c° 



da 

Because of Q3.66p . (|3.67j) and (j3.69[> . / £ C]} and hence by Proposition ^ . 21 monjn Sk ]f 6 C\ 1 . Moreover 

i[o, a „]/llcD < Cv$ ||(mon[o (afc ]/)'||co < c v S. 



mom. 



One can then conclude that the map T defined by the right hand side of (|3.59p maps X uk x X vk X X a k 
into itself. 

Choosing r/k large enough (but independent from 6) it is possible to prove the contraction property. 
One has to exploit all the previous estimates and also properties (|3.9|) and (|3.69|) . which are needed to 
handle the second and the third component of T. 

It turns out that the Lipschitz constant of T (i.e., the constant in the contraction) is uniformly bounded 
with respect to 5. 

□ 

In the following it will be useful to know how the solution of (|3.59|) depends on the length Sk of the 
interval of definition. More precisely, in order to underline the dependence from Sk, we will denote the 
first component of the solution of (|3.59[) as u s k k . Let 

F k (u k , s k ) := u s k "(s k ), 

where the dependence from the initial point Uk is also made explicit. 
To study how F k depends on Sk we will exploit the following result: 

Lemma 3.6. Fix S « 1 and let s\, s\ such that s\ < s\ < 5. Let (u\, v\, al) and (v? s , v\, a^) be 
the corresponding solutions of (|3.59[) . Then, 

11*4- 4llc°([o,4]) + H w fc ~ v k\\c°([o,sl]) + VkS\\al -Ofc||c([o,*i]) < L k S\sl - s 2 k \, (3.70) 
where L k is a suitable constant which does not depend on S. 

Proof. In the following we will denote by T Sk and T Sfc the maps defined by the right hand side of (|3.59p 
when Sk — and Sk — s k respectively. Moreover, to simplify notations, we will denote by (u 2 a , v k , af)f 

the restriction of the fixed point of T Sk to the interval [0, si]. 
Since (u^, v\, a^) is the fixed point of T Sk , one then has 

11/ 1 1 1\ / 2 2 2\|| ^ 1 M/ 2 2 2 \ msj / 2 2 2\|| 

In the previous formula, K denotes the Lipschitz constant of T Sfc , which turns out to be uniformly 
bounded with respect to K as underlined in the proof of Lemma \'S. 5 1 Moreover, X uk , X vk and X a k 
denote the spaces defined by (|3.65|) when Sk = sj,. 
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Thus, to prove the lemma one actually reduces to prove 

11/ 2 2 2\ rrisl I 2 2 2 \ 1 1 II 2 rjisl ( 2 2 2\ II 

+ K -Tv k (u 2 s , Ufc, ^)|| c o ( [ 0)B i]) 

+ ? /feH cr fc -7^(4, 4, CT^) || c o( [0)< ,irj 
<0(1)<5|4-4|. 

In the previous formula, T u k , T„ fc and T CT fc denote respectively the first, the second and the third 
component of T Sfc . 

Since (u\, v\, <4) is the fixed point of T s %, then 



<rl( T ) = g^: mon [o,^]/fe( r ' u fe- u fe> °D- 

Hence, to prove the lemma it is sufficient to show that 

l|mon [0iS 2]/ fe [i4, v\, a%] - mon [0 , s i]/fe[^, v 2 k , cr|] || CQ[0; s i] < 0(l)S(s 2 k - 4) 

and that 

\\(mon [0> ^fklul, v 2 k , a 2 ])' - (mon [0)S i]/fe[«|, v%\, o|])'||e°[o, s ij < C>(1)<5(4 - 4)- 
This follows directly from Proposition [331 

In particular, the previous Lemma implies that 

\F(u k , sl)-F(u k , ,4)1 <3|4-4|. 
Indeed, assuming for example that 4 — 4' one can write 

14(4) - 4(4)1 < 14(4) - 4(4)1 + 14(4) - 4(4)1 

rk{ul(r), v 2 k (r), 4( T ))<* r 



□ 



(3.71) 



< L k s(4 - 4 

<(L k S + 2)(s 2 -sl) 
Take 4 = in the previous estimate: since 

\F k (u k , sl)-F k (u k ,s 2 k )\: 



one gets 



rk(u 2 k (T), vI(t), 4(r))dr 
f k (u k , 0,0) dr 

si 

+ I [f k (u 2 (r), v 2 (r), 4(r)) - f k {u k , 0, 0)}dr 
Jo 

= f k (u k ,0, 0)4 + 0(1)54, 
F k (u k , s k )=u 2 k {s k ) = u k + f k (u k , 0, 0)s k + O{\)5 2 . 



(3.72) 



(3.73) 

Thus, the function F k defined before is differentiable at s k — and the gradient is the eigenvector 
rk(u k ) = f k {u k , 0, 0). 

The function F k depends Lipschitz continuously on u k too: 



52 



Lemma 3.7. Fix S such that < 6 << 1. Take u\ u\ £ M. N such that \u\ — u\\ < aS for a certain 
constant a. If \sk\ < S, then 

\F(ul 8 k ) - F(ul s k )\ < L k \u\ - t*2l (3-74) 

for a suitably large constant Lk . 

Proof. The value F(u\, Sk) is the fixed point of the application T defined by ()3.59|) . To underline the 
dependence of T on Uk we will write 

T : l w x X u k x X v k x X ak — » X u k x X v k x X ak 

The map T depends Lipschitz continuously on u k , but one cannot apply directly the classical results 
(see e.g. [IT] ) on the dependence of the fixed point of a contraction from a parameter because of a 
technical difficulty. Indeed, the domain X u k x X v k x X a k of the contraction depends on u k . 

To overcome this problem, one can proceed as follows. Define X\ k and X 2 k the spaces defined as 
in (|3.65p with u — u\ and Uk = u 2 respectively. If (u\, v k , cr^) is the fixed point of the application 
T{u\, •) defined on X* k x X vk x X ak , then 

114 - u k \\c° < \\ul - ul\\ c o + \u\ -u\\< 2s k + aS<(2 + a)S. 

Choosing c u > (2 + a), one gets that (u\, v\, al) belongs to X 2 k x X v k x X a k- Thus, 

In the previous expression, (u k , v 2 , a 2 ) is the fixed point of the application T(u k , •). One can then 
proceed as in [17] and conclude that 

11(4) V L 4) - ( u li Vk, o-l)\\xz k xX vh xx„ k < Lk\ul-u 2 k \ (3.75) 

□ 

7. Analysis of the purely stable component 

The purely center component satisfies Q3.64[) . 

Lemma 3.8. Fix 5 such that \V S (§)\ < <5 << 1. Then system (|3.64p defines a contraction on the space 
X^ x X%, where 

X» := [u s e C°([0, +oo), R N ), \\u s \\ us < m u d} X s v := {v s 6 C°([0, +oo), R*" 1 ), < mj]. 

The constants m u and m v do not depend on 6 the norms \\ ■ \\ us and \\ ■ \\ vs are defined by 

\\u s \\ us :=r; s sup{e C2: / 2 | U s (x)|} \\V S \\ VS := sup {e cx ' 2 \V s (x)\}, (3.76) 
where i] is a small enough constant which does not depend on 5 and c is defined by (|3.40p . 

Proof. We know that A is has eigenvalues \i(u k ) ■ ■ ■ \ k -i(uk)- Relying on (|3.34j) . one obtains 

\e Kx V s (0)\ < e- cx / 2 \V s (0)\ < e- cx ' 2 5. 
Moreover, if u s £ X^ and V s £ X%, then 

\A(u k +u s ,0, V s ) -A\ < 0(1)5, 
where 0(1) denotes (here and in the following) some constant that depends neither on n nor on 5. 
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The columns of R s (u k , 0, 0) are the eigenvectors r\{uk) . . . rk-i(uk), which are all unit vectors. Thus, 
if u £ X s u and V s £ X*, then 

\R s (u k +u s ,0, V s )\ <0(1). 

From the previous observations one gets that the application defined by (|3.64[) maps X^ x X% into 
itself. 

To prove that the application is a contraction, take u ls , u 2s £ X^ and V ls , V 2s £ X%. Then 

R s (u k + u ls (y), 0, V u {y))V ls {y)dy- J R s (u k + u 2s (y) , 0, V 2s (y))V 2s (y)dy 

< 0(l)S\\u ls - u 2s \\ us + 0(l)S Vs \\V u - V 2s \\ vs + 0(l)i h \\V ls - V 2s \\ vs . 

Choosing r] s sufficiently small, one can suppose that 0(l)rj s < 1/4 in the previous expression. Moreover, 
one can suppose that S is small enough to have that 0{\)8 < 1/2, 0(l)Sr] s < 1/4. 

Finally, 



,72 



»A(x-y) 



A' 



A(x-y) 



■(ui a (y) + « fc ,0, V ls (y)) -K]v ls (y)dy 
i s (u 2s (y)+u k , 0, V 2s (y)) ~h\v 2s {y)dy 



< 0(l)-\\u ls - u 2s \\ us + 0(l)S\\V ls - V 2s \\ vs . 

Vs 

Assuming that S is small enough, 0(l)8/rj s < 1/2 and 0(1)8 < 1/2. Thus the map is contraction and 
the constant of the contraction is less or equal to 1/2 uniformly for <5 — > + . □ 

The solution of (I3.64|) depends on the parameter V s (0). The regularity of (u s , V s j with respect to 
V s (0) is discussed in the following lemma. 

Lemma 3.9. Fix 8 « 1 and let 7/(0), V 2 (0) two initial data such that (7/(0)1, \V 2 (0)\ < 8. Let 
Vg) and (u 2 , V 2 ) the corresponding solutions of (|3.64|) . Then, 



\\u\ - u 2 s \\ us + \\V ls - V 2s \\ vs < L s \V s \0) - V 2 (0)l 



(3.77) 



where L s is a suitable constant which does not depend on 8. Moreover, ifV 2 (0) = 0, then u 2 = 0. 
Also, 

ul(0) = R s (u k , 0, 0)A- 1 V ;j 1 (0) + <D(1)8 2 . (3.78) 
Equation (|3.78p guarantees, in particular, that the application 

F s {u, V s (0)) ^u s (0). 

is differentiable with respect to V s (0) when V s (0) = and that the jacobian is the matrix R s (u k , 0, 0)A _1 , 
whose columns are the eigenvectors ri(u k ) / X\(u k ) . . . r k ^i(u k ) / \k-i{u k ). 

Proof. Let T the application defined by the right hand side of (|3.64p . To underline the dependence on 
the parameter V s (0), we write 



and denote by T u and T v respectively the first and the second component of T. For every (u s , V s ) £ 



X us x X vs , 



T u [u s , 



V„ V}{Q)-T u (u Sl V s , V 2 (0) 
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and 



T, 



'„(u a , V s , V}(0) ~ T v (u s , V„ K 2 (0)) 



< \v}-v?\. 



Hence, (|3.77l) holds with L s := 1/(1 — K), where K is the Lipschitz constant of T and hence it is 
smaller than 1 because T is a contraction. Moreover, from the proof of Lemma 13.81 it follows that k 
is bounded away from uniformly with respect to <5. This concludes the proof of the first part of the 
lemma. 

To prove the second part, we observe that 

V} (x) = e^VXO) + 0{l)5 2 e- cx / 2 

and that 



\R s (u k + ul o, v; 1 ) - R s (u k , o, 0)1 <0(1)S. 



Hence, 



1(0) = 



+ CO 



R s (u k , 0, 0)V s i (y)dy 



R s {u k +ul, 0, V^ 1 ) - R s (u k + u s , 0, V 8 ) V}{y)dy 



R s {u k , 0, 0)e Ay V s 1 {0)dy + 0(1)5 2 / e~ cy/2 dy + 0(1)5 2 

J+oo 



dy 



+ 00 



= R s {u k , 0, 0)A- 1 y s 1 (0) + O(l)<5 2 . 
This completes the proof of the lemma. 



□ 



To study the dependence of F s (u k , V s (0), s k ) on u k one can proceed as follow. Fix u\, v? k £ R N close 
enough and denote by (u*, V^, 1 ) and (u 2 , V 2 ) the solution of (|3.64[) corresponding to u k = u\ and 
u k = u 2 respectively. Then classical results (see e.g. [T7]) on the dependence of the fixed point of a 
contraction on a parameter ensure that 



4\\us + \\V ls -V 2s \\ vs <L s \u k -ul\ 



for a suitable constant L s . In particular, 



\F°(ul V s (0)) - F s (u 2 , V°(0))\ < L s \u\-u 



(3.79) 



(3.80) 



The component of perturbation 

Assume that v k (s k ) < and define a as in (|3.62p : 



a(r) 



r v k( s ) 



ds . 



Let 



s := max-jr € [0, s k ] : v k (r) = 0}, 
we can now prove that a(r) < +oo if and only if r > s. 
Indeed, for s > s it holds 

and hence 



-Wfc(s) = -Ufe(s) + v k (s) < c v 5(s - s) 



a(r) 



v k (s) 



ds > 



1 



T c„(5(s-s) 



(is = 



Let 



p : [0, +00H& s fe ] 
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the inverse function of a, (3 o a(r) = r. If Vk(sk) = 0, we set 0(x) = s k for every x. 
The component of perturbation (U, q, p) is the fixed point of the transformation 



U(x) 



q(x) 



R s [u k o[3(y)+u s (y) + U(y), v k o 0(y) + q(y), V s {y)+p{y), Oj 
-R s (u k + u s (y), 0, V s (y), o)]v s (y)dy 
R s (u k of3(y)+u s (y) + U(y), v k of3(y)+q(y), V s (y)+p(y), Q}p{y)dy 
f k (u k o(3(y) + u s (y) + U(y),v k o /3(y) +q(y), V s (y) + p{y) , o) 
-f k (u k o (3{y), v k o /%), 0, o) v k o (3(y)dy 
+ [ f k (u k o[3(y)+u s (y) + U(y),v k of3(y)+ q (y),V s (y)+p(y),0)q(y)dy 

l+oo V 7 

4>k(uk°/3(y) +u s (y) + U(y), v k o (3{y) +q{y), V s (y) +p(y), a o /3J 
~4>k(uk o (3{y), v k o /3(y), 0, cr o/3J v k o [3{y)dy 
+ / 0fc(w fc o(3(y) +u s (y) + U(y), v k o (3(y) + q(y), V s (y) + p{y), o k o (3^q(y)dy 



p{x) 



A 



(u k o /%) + u s (y) + U(y), v k o /%) + V^fj/) + 
+ ^ e A( "- y) [A(u fc o /3(y) + u s (y) + U(y), v k o /3(j/) + ^(y) 
-Afu a (tf),0), K(y))l^ s (2/)d2/ . 



A 



p(y)dy 



In the equation for p we used the same notation as in p.64[) : A = A s (u k , 0, 0). 



(3.81) 



Lemma 3.10. Let (u k , v k , a k ) and (u s , V s ) be as in Lemma \3.5\ and \ S.8\ respectively. Then there 
exists a unique solution (U, q, p) of (|3.8ip belonging to Xjj x X q x X p , where the spaces Xjj , X q and 
X p are defined as follows: 

X v := [u G C°([0, +oo[, R N ) : \U(x)\ < fc^e^ 72 }, 
X p := {p G C°([0, +oo[, M^ 1 ) : |p(a:)| < fc p< 5 2 e - ra/2 }, 

and 

AT, := {q G C°([0, +oo[, R) : |g(a;)| 9 < k^e^^X 

In the previous expressions, the constants kjj, k q and k p do not depend on 5. 
The space Xjj is endowed with the norm 

\\U\\ u = r lu snv{e cx l A \U{x)\) 

where rjjj is a suitable constant which does not depend on 5. The spaces X p and X q are endowed 
respectively with the norms 

|M| p = sup{ e c */ 4 b(x)|} 

X ^ ' 

and 

\\q\\ q = S u V {e^\q{x)\). 
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Proof. It is enough to prove that the right hand side of (|3.8ip defines a contraction on Xjj x X q x X p . 
As in the proof of Lemma [3~8I one observes that A is a diagonal matrix with all the eigenvalues less or 
equal to — c. Thus, 

|e Ax £| < e- cx \t;\ V, f e R k -\ 

Moreover, from Lemma 13.81 one gets that for every x, 

\V s (x)\ < m v Se~ cx/2 \u s (x)\ < m u 8e~ cx / 2 . 
Another useful observation is the following: 

v k o (3{y)dy < s k < S. (3.82) 

+00 

Indeed, if v k (s k ) < 0, then one can exploit the change of variable r = j3{y). Recalling that (3 is the 
inverse function of a, which is defined by (|3.62p . one has 

' Vk o P(y)dy < s k = v k (r)a' '(r)dr = / dr = (s k - s) < s k . 

+oc J s_ J s_ 

If v k (s k ) = 0, then v k o (3 = and hence (|3.82p is trivial. 
One can also assume 

|fo(u*,0, 0, Xk{u k )\<K6 (3.83) 

and hence 

\<j> k (u k op(y)+u s {y) + U{y), v k c/3( y )+q(y), V s (y)+p(y), a k o 0) | < 0(1)6. 
Moreover, it holds 

K(u k op(y)+ Us (y) + U(y),v k of3(y)+q( y ), V s (y) +p(y)) - A| < 0(1)S 

Since A = A s (u k , 0, 0), this follows from the regularity of A and from the estimates 

\u k -u k \, \u s \, \v k \, \q\, \V S \, \p\<0(l)S. 

Exploiting the previous observations, by direct check one gets that the right hand side of (|3.81[) belongs 
to Xjj x X q x X p . 

To prove that the map is actually a contraction, one has to exploit again the previous remarks, combined 
with the following: for every x, 

_ U^ x ) < J_||f/1 _ U 2 \\ v e- cx /\ Iq 1 - q 2 \{x) < \\q l - q^e^ 
Vu 

and 

Ip 1 -p 2 \(x)< Wp 1 - V 2 \\ p e~ cx l\ 
One can then check that the map is a contraction, provided that r\u is sufficiently small. □ 

Since (u k , v k , a k ) and (u s , V s ) depend on the parameters s k and V s (0) respectively, then also {U, q, p) 
does. Let 

FP(u k , V s {0), s k ) := U(0). 
The following lemma concerns the regularity of F p with respect to (V s (0), s k ). 
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Lemma 3.11. The function F p is differentiable with respect to (V s (0), SfcJ at (V S (Q), Sk) — and the 
jacobian is the null matrix. 

Also, let (V}(0), s\) and (V 2 {0), s|) such that \(V}(0), s\)\, \(V 2 {0), s|) | < 5. Then, 

V s \0), 4) - FP(ul V 2 (0), 4)| < L p (|4 - 4| + 1^(0) - y s 2 (0)|). (3.84) 
Proof. By construction, 

\U(0)\ < 0{1)S 2 . (3.85) 
This implies the differentiability at (14(0), si) = 0. 

To prove the second part of the lemma, we will focus only on the dependence from V s (0) and s/.. The 
proof of the Lipschitz continuous dependence from is completely analogous. 

We can observe that (U, q, p) is the fixed point of a contraction T, defined by the right hand side of 
()3.8ip . To underline the dependence from (V s (0), si) we write 

T : X u x X q x X p x K fe " 1 x M -»■ X u x X q x X p 

and denote by Tjj, T q and T p respectively the first, the second and the third component of T. Since 
the Lispchitz constant of T is not only smaller than 1 but also bounded away from uniformly with 
respect to 5, then to prove (|3.84|l it is enough to prove that, for every (U, q, p), 

WMV^Q), si U, p, q) - Tu(V?(0), si U, p, q)^ < 0(1) (|4 - 41 + 1^(0) - V 2 (0)\) 
11^(^(0), 4, U,p, q)-T q (V 2 (0), si U,p, q)\\ q < 0(1)(|4 -4| + \V s \0)-V s 2 (Q)\) (3.86) 
\\T p (Vt{0), 4, U, p, q) - T p (V s 2 (0), s 2 , U, p, q )\\ p < o(i)(|4 - 41 + 1^(0) - V 2 (0)\). 

Let 

ai[T] ) = ~ I k ^m ds a2iT) = - V W) ds 

and (3 1 and j3 2 the inverse functions. Just to fix the ideas, let us assume that s\ < s|, then (3i(0) = 
4 < /^fe(O) — s l Exploiting Lemma I3~B1 one gets that, as far as Pi{x) < P%{x), it holds 

^^ = 4°/3 2 -4o/3 1 = 4o/3 2 -4°/3 1 + 4°/3 1 -4°/3 1 

< ||4-4llco+^(/5 2 -/3 1 ) 
<L k 5{s 2 k -sl)+c v 5{(3 2 

and hence 

(*) - /3 2 (x)| < [^(4 - 4) + (4 - 4)1 ^ 5x - ^(4 - 4) < o(i)(4 - 4)^ fe (3.87) 

HP 1 > (3 2 , then 

P l (x) - fi 2 {x) < 4* [e 2 ^* - *) - lj (4 - 4), 

where 

5 = max{y < x : /3 x (y) = /3 2 (y)}. 

Thus, estimate (|3.87[) still holds. 
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fk. 



Exploiting (|3.87[) , one gets 
\ul o (3 2 (x) - u\ o /3 1 (a;) | < \u\ c /3 2 - u\ o /3 1 1 + | U 2 o /3 1 (z) - uj o /3 1 (z) | 

< \\u\ - 4|| c o + 2|/3 2 - < L fc 5(4 - 4) + 0(1)(4 - 4)e^ fa 

<0(l)(4-4)e 5 ^ 
and, by analogous considerations, 

|4 °/? 2 04 - ^0^(1)1, \a\ o(3\ x ) - alof{x)\ < 0(l)( s 2 - S \)e^ 5x . 
To prove (|3.86|) . the most complicated terms to handle are in the form 

' (u\ o p\y) + ul(y) + U(y), 4 o (3\y) + q(y), V^y) +p(y), o) 
■■(u\ot3\y), v\op\y), 0, Q)]vlop\y)dy 
■ (ut o /3 2 (y) + u 2 s (y) + U(y), v\ o /? 2 (y) + V s 2 (y) 

- f k (ul o /3 2 (y), 4 o /? 2 (y), 0, 0)] 4 ° 

? fc (4 o /3 x (y) + 4(y) + U(y), 4 o ^(y) + q(y), V}(y) + p(y), 
- f k (ul o /3\y), 4 o ^(y), 0, 0)] (4 o fiy) - 4 ° /3 2 (2/))^ 

r fc (u\ o /3 1 (y) + 4 (y) + U(y) , 4 o /3 1 (y) + g(y) , 7/ (y) + p(y) , ' 
-r^o^y), 4.o/3 1 (y),0, 0)" 

f fc (u 2 o /3 2 (y) + ,i 2 (y ) + U(y ) , 4 o /3 2 (y ) + q(y ) , if (y ) + p(y) , 
-f k (ul o/3 2 (y), 4 o/3 2 (y), 0, 0)] |4 o/9 2 (2/)^ 



< 



Exploiting (|3.89[) the first term in the previous sum can be estimated as follows: 

(u\ o p\y) + 4(y) + U(y), 4 ° /? X (|/) + Q{V), V}{y) +p(y), 
- f k (u\ o /3 1 (y) , 4 o /3 1 (y), 0, 0) ] (4 o /3 1 (y) - v\ o /? 2 (y)) dy 

0(l)(4-4)e^^ 



-cy/2 



< , e 

'+00 

<0(l)«-.i)fc-«/ 4 



m u (5 + k u S 2 + fc g <5 2 + m„<5 + k + p5 £ 



(3.88) 



(3.89) 



(3.90) 



where we have supposed that S is small enough to obtain c v S < c/4. 
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To handle with the second term in (|3.90[) one can observe that 
k (ulof3\y) + u] (y) +U(y),vla (3 1 (y) +q(y), V} (y) +*>(»), 



D u f k > 
• 1 



dt 



{u\op\y), vlo0\y), 0, 

(u\ o f3 1 {y) + tul(y) + tU(y), v\ o p\y) + tq(y), tV}{y) + tp(y), o) [ul(y) + U(y) 

v f k (u\ o P\y) + tu\{y) + tU(y), v\ o 1 {y) + tq(y), tV}(y) + tp(y), o) [q(y)]dt 

+ J D v f k (u l k o (3\y) + tu\{y) + tU(y), v{ o f3\y) + tq(y), tV}{y) + tp(y), o) [V 3 1 (y)+p(y) 
Moreover, 

£> u f fc (4 o p\y) + tu\{y) + tU(y), v\ o /3 x (y) + fctfy), iF/(y) + tpfc,), o) [ U \{y) + U(y) 
D u f k ( u\ o /3 2 (y) + iu 2 (y) + tU(y), v 2 o /3 2 (y) 



tg(y),rt?(y)+ip(y),0 <(y) + C/(y) dtW k of3\ y )dy 



< 



oo I ^0 



u f fe (4 o ^(y) + tuiCtf) + «i ° ^(y) 



+ tq{y), tV s 1 (y)+tp(y), ^(y)-u 2 (y) diU 2 o/? 2 (y)dy 







/"I 


i ■ 



Z?„f*(uJ o ^(y) + *ui(y) + tU(y), v\ o ^(y) + i$(y), tV/fo,) + tp(y), 
(u 2 o /3 2 (y) + tu 2 (y) + tU(y), v\ o /3 2 (y) 
+ tq(y), tV s 2 (y) + tp(y), o)] [u 2 s (y) + U(y)}dt \v 2 o /? 2 (y)dy 



< 



+ cx) I JO 



< 



0(1)L.|V?(O) - V 2 (0)\e~ c y/ 2 dtjc v Sdy 

0(l)(s 2 - .sDe^ +i s |K 1 (0) - K'WIe"^/ 2 ]^-^/ 2 *^ 

0(1)^(1^(0) - K 2 (0)| + (4 - 4))e" ca; / 4 







/"I 

/+oo 


i ■ 



By analogous considerations, one can conclude that the second term in (13.90|) can be bounded by 

0(1)^(1^(0) - y s 2 (o)| + (4 - 4))e- cx / 4 . 

One can also perform similar estimates on all the other terms that appear in 

\\Tu(V}(Q), 4, U, P , q)-Tu(V 2 (0), 4, U, p, q)\\ w 

\\T q {Vt(0), si U, p, q) - r ? (y s 2 (o), 4, U, p, q)\\ q 



and 



| T p (1^(0), 4, U, p, q) - T p (V 2 (0), 4, U, p, q) 
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and conclude that (|3.86p holds. □ 

Set V s (0) = (s% ... Sfe-i) and define 

F(u k , s u ...s k ) :=F k (u k , s k ) + F s {u k , s u . . . s k -i) + F p (u k , Sl ,...s k ). (3.91) 
Combining (|3.70|) . ()3.77j) and (|3.84|) one gets that F is Lipschitz continuous with respect to (si . . . s k ): 

\F(u k , si, . . . s k ) - F(u k , §i, . . . s k )\ < L(Jsi - Si| + . . . |sfe - s k \) (3.92) 
for a suitable constant L. 



Moreover, F is differentiable at (si . . . s k ) = because of (|3.73[) , (|3.78p and (|3.85p . The colums of the 
jacobian matrix are the eigenvectors ri(u k ) . . . r k (u k ). 

Finally, F is also Lipschitz continuous with respect to u k , namely 

\F(ul Sl ... s k ) - F(u 2 k , ai ... s k )\ < L\u\ - u\\. (3.93) 

Define 

t&uo, Sl ... s N ) := F(T*£ o ■ • • o T s n Ju ), Sl ... s k ) . (3.94) 

We recall that the composite map T^+^o- • •oTj^uo) is Lipschitz continuous and differentiable at (s k+ i . . . sn) — 
(0 . . . 0), with the columns of the jacobian given by the eigenvectors r k+ i(uo) ■ ■ ■ ^(uq)- This was proved in 

m- 

Combining this with (|3.93p and p.94p one gets that 4> is Lipschitz continuous. Moreover, it is differentiable 
at (si . . . sn) and the columns of the jacobian are the eigenvectors ti(uq) . . . r^iuo). Thus, thanks to (|3.4|) . 
the jacobian is invertible and hence, exploiting the extension of the implicit function theorem discussed in 
[2"4"| (page 253), the map <fi(uo, ■) is invertible in a neighbourhood of (si . . . sjv) = (0 ... 0). 

Thus, if Ub is fixed and is sufficiently close to uo, then the values of s\ ... sjv are uniquely determined by 
the equation 

Ub = 4>{uo, si . . .s N ). (3.95) 

We will take u& equal to the boundary datum imposed on the parabolic approximation (|3.ip . Once (s\ . . . sn) 
are determined one can determine the values of the hyperbolic limit u(t, x) a.e. (t, x). This can be done 
glueing together pieces like p.33p . Here we will just make a remark concerning the trace u of the hyperbolic 
limit on the axis x — 0. Let s k determined by equation (|3.95p and let (u k , v k , <r k ) solve the fixed point 
problem (|3.59p . Define 

s := min{s e [0, s k ] : u k {s) = 0}, 
then the trace u of the hyperbolic limit on x — is 

u = u k (s) (3.96) 

The following theorem collects the results discussed in this section. 

Theorem 3.2. Let Hypotheses QJ [7J j and\Q hold. Then there exists S > small enough such that 
the following holds. If \uq — tt&| « 5, then the limit of the parabolic approximation (|3.ip satisfies 

u b = <f>(iio, si . . . s N ) 

for a suitable vector (si . . . Sjv)- The map cj> is given by 

0(«o, sl.. s N ) := F(T*+\ o ■ • • o T s n n (u ), Sl... Sfe) , 

where T^^ . . . T^ N are the curves of admissible states defined in j^jj, while the function F is defined by (|3.9ip . 
The map <f> is locally invertible: given uq and Ub, one can determine uniquely (s% . . . sn). Once (s\ . . . sn) 
are known, the value u(t, x) assumed by the limit is determined for a.e. (t, x). In particular, the trace u of 
the hyperbolic limit in the axis x = is given by (|3.96p . 
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4 The characterization of the hyperbolic limit in the case of a 
singular viscosity matrix 

The aim of this section is to describe the limit of the parabolic approximation 1)) when the viscosity matrix 
B is not invertible. The precise hypotheses assumed in this case are introduced in Sections 12.31 and 12.21 In 
particular, they guarantee that it is sufficient to study the family of systems 



E(u e )u £ f + A(u e , eu%)u% = eB(u e )u £ xx u £ € 
&{u E {t, 0)) =g u%0,x)=u o , 



pN 

(4.1) 



with E, A and B satisfying Hypothesis 6 given by Definition l2.fi 

The exposition is organized as follows. In Section 14.11 we introduce some preliminary results, while in 
Section 14.21 we give a characterization of the limit of (|4.ip . More precisely, in Section 14.2.11 we briefly recall 
the characterization of the limit in the case of the Cauchy problem: 



E{u e )u E t + A(u £ , £u £ x )u x = eB(u £ )u s x 
Uq x > 
u~ x < 



w e (0, x) 



We refer to [7] for the complete analysis. In Section P4.2.2I we introduce a lemma on the dimension of the 
stable manifold of the equation 

A(u, u x )u x = B(u)u XXl 

which is the equation satisfied by the steady solutions of (|4.1j) . Such a lemma gives an answer to question 
which was left open in [Hj. In Sections 14. 2. 31 we give a characterization of the limit in the non characteristic 
case, i.e. when none of the eigenvalues of E~ 1 A can attain the value zero. The boundary characteristic case 
occurs when an eigenvalue of E~ x A can attain the value and it is discussed in Section [4.2.41 Finally, in 
Section f4.3l we discuss a transversality lemma, a technical result which is needed to complete the analysis in 
Sections |4~!Q1 and [4~24l 



4.1 Preliminary results 

This section introduces the preliminary results needed in Section f4. 21 The exposition is organized as follows: 
Section B. 1.1 1 gives some results about the structure of the matrix A(u 7 u x ), all implied by Hypothesis [U and 
in particular by Kawashima condition. In Section [4.1.21 we rewrite equations 

E(u)u t + A{u, u x )u x = B{u)u xx (4.2) 

and 

A(u, u x )u x = B(u)u xx (4.3) 

in a form easier to handle. 

Some notations have to be introduced: given a real parameter a, let 

A(u, u x , a) := A(u, u x ) - aE{u). (4.4) 

Consequently, both systems (14. 2| and (14. 3| may be written in the form 

A(u, u x , a)u x = B(u)u xx . (4.5) 

Finally, let 



( A n (u,a) A 12 (u,a) \ = ( A n (u) - aE n (u) A 12 {u) - aE 12 (u) 
\ A 2 i(u, a) A 22 (u,u x ,a) J \ A 21 (u) - <tE 21 (u) A 22 (u, u x ) - oE 22 (u) 




and 

S = (id, z) t , w e R N ~ r , zeW (4.7) 
be the block decompositions of A(u, u x , a) and of S £ l w corresponding to (|2.8p . 
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4.1.1 Some results about the structure of A(u, u x , a) 

Even if it is not explicitly specified in the statement, all the results of this section presuppose that Hypothesis 
[H holds. 

The first observation is immediate, but useful: 

Lemma 4.1. Let Aij(u, u x , a) be the blocks of A(u, u x , o) defined by (|4.6j) . Then 

A 12 (u, a) = Aj^u, a) ^n(w, a) = A n (u, a) 

Some further notations are required: let Po(u, a) denote the projection of R w " r on ker Au(u, a). Thanks 
to the third condition in Hypothesis the dimension of this subspace is constant in u: we will denote it 
by q. The operator Pq can be identified with a matrix Po(u, a) € M. qx ( N ~ r h Similarly, the projection 
on the subspace orthogonal to ker An(u, a) is identified with a matrix P±(u, a) 6 M( JV ~' r ~' z ) >< ( jv ~' r ). The 
decomposition (|4.7j) is hence refined setting 

to = P Q w w = P±w to we R N ~ r ~ q . (4.8) 

Using the previous notations and recalling Lemma |4. 11 the product A(u, u x , <r)S writes 

{ A 21 \u, J) A 22 (u, uJ,a))[ Z ) = \ J 4 , , {Aifiu, a U , (4.9) 

V A 2l( U ) A 2l(«, CT ) A 22(U, <7, U x ) J \ Z ) 

where 

A' 21 {u, a) = A 2 i(u, <r)P T (u, <r) S M rxq A^u, a) = A 21 (u, «r)Pj(u, a) 6 ^x(Af-r-<j) 

and 

An(«, a) = Pj.(«, a)Au(u, a)Pj(u, a) e M( JV - r -^ x ( JV -'"-^ 
It is now possible to state the following result, a consequence of Kawashima condition: 
Lemma 4.2. The matrix A 2 i(u, a) € M rX9 has rank q. Thus, in particular, q < r. 
Proof. First of all, one observes that 

E~ l A(u, 0)3 = o-S A(u, 0, cr)S = 
and hence Kawashima condition writes 

ker(A(u, 0, a)) Hker(P(u)) = {0} Vu, a. (4.10) 

We claim that this implies 

A^w = => to = 0. (4.11) 
Indeed, suppose by contradiction that there exists to ^ such that A^to = 0. Then 

/ (Ai^ \ / w\ ( \ / \/»\ 

A n (A'I 21 ) T = = , 
V A\ x A T I 21 A 22 J \ J \0 J V 6 / V / 

which contradicts (|4.10| . Hence (|4.1ip holds. 

From (|4.1ip one then deduces that A 2l has q independent columns and hence it admits a submatrix of 
rank q. □ 
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In the following we will isolate an invertible submatrix of A T 2l (u, a) and we will denoted it by «u(u, a). 
More precisely, an is obtained from A^^u, a) selecting q independent rows. By continuity, one can suppose 
that the rows selected are the same for every u. 

With considerations analogous to the ones that lead to (14. 9p . one can write 



P T (a)z + P T (a)z, 



r — q)xr 



where P € M 9Xr is the projection such that an = A\ X P T and P G 
subspace orthogonal to the image of W through P. 

Thanks to the previous considerations, one obtains that A(u, u x , <r)3 writes 

/ afx(u, a) a^(u, a) \ ( w ^ 

An (u, u) a^iu, a) o,2 2 (u, a) w 

On(u, a) au(u, a) au(u, u x , a) a^u, u x , a) z 

\ a 2 i(u) a 22 (u, cr) Qt2i (w, u x , a) a 22 (u, u x , a) J \ z, J 



is the projection on the 



(4.12) 



with 



ai2 e 
an G 



]X(N — r — q) 



, a 2 2 G 



ir(r-g)x(JV-r-g) 



?x «, aai GlV 
Z G 

The corresponding decomposition of -B(u) is 



r(V-g)xg 



, «22 G 



X (r-g) 



Z E 



with 



/ 



6u(u) 6 ia («) 

V o o 02 i(m) 6 22 ( M ) y 



6xi G M« x<z , 6i2 G M 9X ( r -«), 6 21 G M^"^* 9 , 6 22 e M< r -«> x < r - 9 >. 



(4.13) 



It is worth underling explicitly that formulations (I4.12p and (|4.13|) include in particular the limit cases q = 
and q = N — r: indeed, it is enough to assume w = w, z — z and w — w, z = z respectively. 

Lemma 4.3. The matrix 

a2iari 1 (frii( a ii)~ la 2i ~ &12) - &2x(axx)~ la 2X + b 22 

is invertible. 

Proof. By contradiction, assume that there exists an (r — q)-dimensional vector £ ^ such that 
^a 2 ia n 1 ^6n(af 1 )~ 1 a 2 r L - o i2 ) - 6 2 i(an)~ la 2X + ^22^1*= 0. 



It follows that 



Define 



(a210.11 , —I, 



r-q) 



Oil bi2 
621 622 



l ll) 



'21 



I r — q 



(4.14) 



D= ( Kl ) a 21 



r—q 



then (14. 14[) implies 

Since D £ 7^ 0, this is a contradiction because of the second condition in Hypothesis [H 



□ 
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4.1.2 The explicit form of system (|4. 5[) 

The following result guarantees that system (|4.5|) . and hence systems (|4.2p and (I4.3[) can be reduced to an 
explicit form like (|4.15p 

Lemma 4.4. There exists a sufficiently small constant 5 > such that if \u — uq\ < S and \u x \ < $, then 
the following holds. System (|4.5|) can be rewritten as 

f u x = (w(u, z, a), w(u, z, a), z{u, z, a), zj {A-.lS) 
{ z x = f(u, z, a) 

for suitable functions 

w : R N x R r - q xl^r w:R N x R r - q xR^ M. N ~ r ~ q z :R N x M. r ~ q xl^M' 

/ : m N x m r - q xr^ m r - q . 

For every a and for every u the following holds. 

w(u, 0, a) = w(u, 0, cr) = z(u, 0, ct) = /(u, 0, cr) = (4.16) 
Also, these functions satisfy 

Duz\ u=aQ; - =0>u = Drz]^ = D s z | ^ = -(^"Xi (4-17) 



and 



D u w\ _ . . =0 D(jW _ . . =0 

^H^un 5=0 <r = (Ol^-Ds^ + a 22) = -^U ( - 4( a u)~ la 21 + a : 



22 ; 



Moreover, 



D u w\ _ - n =0 D a w\ - - n =0 

= <hl (a 12 A^ (aj 2 - af 2 (af 1 ) _1 a 2 r 1 ) + anfan)"^ - ai 2 ) + a^ 1 ^6i 2 
- &ii(an) _1 a 2 i^ ^21%/ (&ii( n) _1 °2i _ & is) - & 
+ &22^ ^lan^a^nO 3 ^ - af 2 (af 1 ) _1 a 2 " 1 ) + an(a^ 1 ) _1 a 21 - ai 2 ) 
+ a22^ri 1 ( a i2( a n) _la 2i - a 22 - «2i(au) _1 a2i + 022)^ 



(4.18) 



/ T \ — 1 T (4.19) 

21 (an) <% y ' 
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and 



D u f 



I U—Uq , 2—0, (7 

D sf\ - - n 



D a f 



u—uq , z—0, a 

b2iDiZ + b 2 2 - a 2 i 



- a21a.il [aizDzw + anD=,z 



iz + a 11 1 &i2 
azzDzw + ct2iDzW + a 22 ^ 

^a 2 iaii (b 1 i(a r [i)~ 1 a2i - b 12 ^j - b2i(aj i y 1 a 21 + &22^ (a21a.11 (a^Af/ (a 

) +a 22 i n 1 (a? 2 (an)" 1 



1„T 
21 



a 12 (a n ) a 21 j +au(a n ) a 

Q 2l(o[l)" 1 fflM + "22 



«21 



. T 
£»21 



T 
a 22 



(4.20) 



Proof. Using notation (|4.12p and (|4.13[) . equation rewrites as 





An(u, a) 



aJi{u, a) 
aj 2 {u, a) 



ali(u, a) 
a22i.ua) 



an(u, a) ai 2 (w, a) a\i{u, w, w, z, z, a) (^(it, w, w, z, z, a) 
\ a 2 i(u, a) a 22 (u, a) a 2 i(u, w, w, z, z, a) a 22 (u : w, w, z, za) j 



\ 


( 


w \ 






w 






z 


) 


\ 


* I 



(° 











\ 
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w x \ 




















w x 








bn(u) 


b\2{u) 










V 





b 2 i{u) 


b 2 2{u) 


J 




\ 


Zx ) 



Hence from the first line it follows 
and therefore 



a 11 (w, a)z + a 21 (u, a)z = 
,-1 t 



z{u, z) = -(a n (u, a)) a 2 i(u, a)z. 
The second line then reads 

An (it, cr)w — aj 2 (u) (aii ( u > <T )) la 2i( u > a )z + a 22( u ) a )z = 
and hence thanks to the invertibility of An(u) 

w(u, z) = Aii(u, a)(al 2 (u)(ali(u, a)) 1 a 21 (u, a) + a 22 (u, afjz. 

Moreover, one has 

(z(u, z)) x = -(a,ii(u, a)) 1 a 21 (u, a)z x - Atf (^(n, a)) 1 a 21 (u, a)ju x = z x (u, w, z, z x ) 

and hence the third line reads 

aii(u)w + ai2(u)w(u, z) + au(u, w, z)z{u, z) + a 21 (u, w, z)z = bn(u)z x (u, w, z, z x ) + b\2{u)z x 

and therefore, thanks to the invertibility of an, in a small enough neighborhood of 
(u = uq, w = 0, z = 0, z x = 0) it is implicitly defined a map w — w(u, z, z x ) such that 



D u w 
Dtw\ 



U=Uq, 2 — 0, Z x =0 
U—Uq, 2—0, Z x —0 



= 



u=So , 2=0 , 5ie= o = «n(&i2-&ii(an) ^21) 



hi \a12An {a 22 - a 12 (a u J a 21 ) + au(o n ) a 2i - a 



GG 



Then the fourth line reads 

a 2 i(u)w(u, z, z x ) + a 22 (u)w(u, z) + a 2 i(u, z, z x )z x (u, z, z x ) + a 22 (u, z, z x )z = b 2 i{u)z x (u, z, z x ) + b 22 {u)z x 

Hence thanks to Lemma 14.31 it is implicitly defined a map z x — /(it, z), which satisfies the hypotheses 
described in the statement of the lemma. □ 

In order to simplify the notations, we set 

a(u, z, cr) := a 2i a^ (ai 2 A^ (a 22 - a( 2 (a{ 1 )' 1 a 21 ) + a u( a ii)~ la 2i ~ a 2i) + a 22^n (^("n) -1 ^ 

- a T 22 - a^rXi + « 22 ) € M^-«) x ^-«) 
b(u, a) := a 21 a^(b 1 i(a{ 1 )~ 1 a 2 i - 612 J - ^i(«u) _1 «2i + ^22 € 



[(r— 9) X (r-g) 

(4.21) 

With this notations, 



= 6 (ito, 0, cr)a(uo, 0, cr). 



£*/ 

Also, if we consider the jacobian 

Dg(w(u, z, cr), w(u, z, cr), z(u, z, cr), z(u, z, a)j 

and we compute it at the point I u — uq, z — 0, cr), we get 



A11 f a i2( a n) la 2i "22 



a x la 12 A^{al 2 - af 2 (afi) 1 aJ 1 )| + a u 1 oni(af 1 ) - a 11 1 a 2 r L |*+ (a!^^ - a 11 1 6ii(a 11 1 ) T a 2 r L )6 ^ \ 

-1 T 

(4.22) 

Finally, the proof of Lemma l4~4l implies the following result: 

Lemma 4.5. Given £ 6 R r_9 , i/ie condition a(u, 0, cr)£ = A&(n, cr)£ is equivalent to A(u, 0, cr) S = XB(u) S, 
where 3 6 R-^ is given by 

( a^a^in ("22 ~ a i2( a ii) _la 2i)C + an^iiCan) -1 ^!*- a n a 2iC + (a n 1& i2 - Ci^nfau ) Ta 2i)F 1 af 

( a i2 (an)" 1 a 2i - a 22 ) I* 
-(afi)- 1 ^! 4 

v s 

(4.23) 

TTte case A = and hence at; = is, in particular, included in the previous formulation. 

Remark 4.1. For simplicity, in the following we assume that 6 _1 a is diagonalizable. Hence, in particular, 
one can find (r — q) independent vectors Si, . . . S r _ 9 such that 

(A - XiB)Ei = 0. 

To handle the general case, one should work with generalized eigenspaces. 
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4.2 The hyperbolic limit in the case of non invertible viscosity matrix 
4.2.1 The hyperbolic limit in the case of the Cauchy problem 

In this section for completeness we review the construction of the Riemann solver for the Cauchy problem 
in the case the viscosity matrix B is not invertible. We refer to [7] for the complete analysis. Also, some of 
the steps in the construction are the same as in the case of an invertible viscosity matrix, which is discussed 
in Section [3.2.11 Thus, here we will focus only on the points where the singularity of the viscosity matrix 
plays an important role. 

The goal is the characterization of the limit of 

E(u-)ul + A{u°, ul)u% = eB(u s )u% x 

u- x < (4.24) 



u E (0, , x) — i 

y 1 ' > \ u x > 

under Hypotheses El El El and El 

The construction is as follows. Thanks to Lemma 14^41 the equation of travelling waves 



B(U)U" = (A(U, U') - aE(ufj U', (4.25) 



is equivalent to the system 



u' = (w(u, z, tr), w(u, z, a), z(u, z, a), z^j 

z' = f(u, z, a) ' ( 4 - 26 ) 

a' = 0. 

When z = then for every u and a f(u, 0, a) — 0, w(u, 0, a) — 0, w(u, 0, a) — 0, z(u, 0, a) — 0. Thus, 
(uq, 0, Ai(wo)) is an equilibrium for (|3.28|) . where Xi(uo) is an eigenvalue of E~ 1 (uo)A(uq, 0). Linearizing 
around such an equilibrium point one obtains 

c 

b^aiuo, 0, a,) ) , (4.27) 


where the exact expression of the matrix c G R Nx ( r ~ q ) is not important here. Let Si(uo) denote the 
eigenvector corresponding to Xi(uo), then ^(uq is in the form given by Lemma 14.51 for some vector £(mo) £ 
M. r ~ q . Then the eigenspace of (|4.27p associated to the eigenvalue is 



yc 



{(m, Vi£i(u ), <t) : u £ R N , v u a G r\. 



One can then proceed as in the case of an invertible viscosity matrix (see Section [3.2.ip . Fix a center manifold 
Ai c , then one can verify that a point(u, z, a) belongs to M. c if and only if 

z = Vi£i(u, Vi, a), 

where ^ G M. r ~ q is a suitable vector valued function such that 

ii(u , 0, Ai(u )) = ii(u ). 

Thus, on the center manifold 

u' = (w(u, z, ai), w(u, z, ai), z(u, z, o-j), zl = Sj(u, Vi, Ui)vi 
for some vector valued function Hj G M. N such that 

Zi(u , 0, Xi(uo)j = 5i(u ). 
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Plugging the relation v! = EiVi into the equation 

B{u)u" = [A(u, v!) - aE{u)]u' 

one obtains 

[BE t + E iv Vi]v ix = ([A - aE]Ei - D u E t EiV^jVi. 
Taking both members dot product Ei one gets 

C-iVix — QiVi 

for suitable functions Ci{u 1 Vi, a) and a,i(u, Vi, a). Thanks to Lemma l4.5[ Ei(uo) is not in the kernel of B 
and hence 

c^uq, 0, Xi{u )j = (Ei(u ), B(u )Ei(u )) > 
Thus, Ci > in a neighborhood and hence one can introduce 4>i '■= CLi/c-i- Also, 

= -(Si(u ), ESi(fio)) < 0. (4.28) 



d(j)j 
da 



(« , 0, Ai(tio)) 

Thus, system the solutions of (|3.28|) laying on M. c satisfy 

it- = ViEi(u, Vi, (Ji) 

v[ = ^(u, Vi, Oi) (4.29) 

with cj> satisfying (|4.28p . One can therefore apply the results in [7]. Proceeding as in Section [3.2.11 one 
eventually defines the i-th curve of admissible states Tg.Uo, which satisfies 



0THt o 



ds. 



= Sj(u ) 

Si=0 



and all the other properties listed in Section f3. 2. II 
Consider the composite function 

i/>(u , s 1 ...s N ) = o ... Tf N u 

With the previous expression we mean that the starting point for is T^uq. It turns out that the 

map 4>{uq, •) is invertible in a neighbourhood of (s± . . . sn) = (0 ... 0). In other words, if u~ is fixed and is 
sufficiently close to uq, then the values of s\ ... sjv are uniquely determined by the equation 

U~ = ip(u , si . . . sjv). (4.30) 

Taking the same u~ as in (|4.24p . one obtains the parameters (si . . . sjv) which can be used to reconstruct the 
hyperbolic limit u of (|4.24jl . Indeed, once [s% . . . sn) are known then u can be obtained in the way described 
in Section f3. 2.11 

Remark 4.2. To construct the curves Tj.ito one proceeds as follows. Suppose that s, > 0, then one considers 
the fixed point problem 



u{t) = u k + / r k (u(0, Ufc(0i a k(0)d€ 


Vk{r) = fk(r, u, Vk, cr fe ) - conc/ fc (r, u, v k , a k ) (4-31) 
a k(r) = ^-y— conc/ fc (r, u, v k , a k ). 
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Denote by (tij, Vj, <7j) the solution of (|4.3ip , whose existence is ensured by results in [7]. Then we define 



TI.Uq = Ui(Si). 

The value Ui(si) is connected by uq by a sequence of rarefaction and shocks. This is actually true only if 
system ()4.28|) is equivalent to ()4.25[) . i.e. if Lemma l4~4l holds true. The problem here is that in (|4.28|) a is a 
parameter, while in (|4.31| can vary. Actually, it turns out that we can still apply Lemma l4~4l provided that 
the dimension of the kernel of A\\(u) — uE\\(u) does not vary in a small enough neighborhood of a — A,(iio), 
i.e. if the kernel of An(u) — Xi(uo)En(u) is trivial. This is certainly a restrictive condition. In [13) the 
authors will address the problem of finding less restrictive conditions. 

4.2.2 A lemma on the dimension of the stable manifold 



The aim of this section is to introduce Lemma 14. 7[ which determines the dimension of the stable manifold 
of the system 

B{U)U" = A{U, U')U', (4.32) 

satisfied by the steady solutions of 

E(u)u t + A(u, u x )u x = u xx . 

In this section we thus give an answer to a question introduced in [44) : more precisely, it is shown that the 
condition called there Hypothesis (H5) is actually a consequence of Kawashima condition. 
Thanks to Lemma (14.41) . system ([4.32p is equivalent to 



u x = (w(u, z, 0), w(u, z, 0), z(u, z, 0), 
z x = f(u, z, 0). 



(4.33) 



Hence, linearizing around the equilibrium point (it, 0), one finds that the Jacobian is represented by the 
matrix 

TO 

b^aiu, 0, 0) 

where a and b are defined by (|4.2ip and m is a N x (r — q) dimensional matrix whose exact expression is 
not important here. Hence the stable manifold of system (|4.33[) is tangent in (u, 0) to 

d d 

V S = { { U + Yl v i™$i{u), Y, : ".£ R }c R N+r - q , (4.34) 

»=i »=i 

In the previous expression 9\, . . . £ are the eigenvectors of & _1 a associated to eigenvalues with 

strictly negative real part. What we want to determine in Lemma (j4.4[) is the number d, i.e. the number of 
eigenvalues of 6 _1 a with strictly negative real part. 

Before stating the lemma, we recall some notations: k— 1 denotes the number of eigenvalues of E~ 1 (u)A(u, 0), 
that satisfy 

Ai(u) < X 2 (u) < • • • < Afc_i(u) < -c < 

for some constant c > 0. The number (k — 1) does not depend on u and u x because of the strict hyperbolicity 
(Hypothesis [3]) . Moreover, let 

A(u,u x )=( A /f\ /f (U \) m=(l ll{ ?\ ffl) (4-35) 

^ ' x > \ A 21 (u) A 22 (u, u x ) J V £ 2iH E 22 {u) J y ' 

be the block decompositions of A and E corresponding to the block decomposition 

*<»> - ( S ><°) 
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Figure 6: the behavior of the roots of (|4.38|) when the perturbation a — > 



k — 1 roots: Render) < N — k roots: Reni(a) > 



+9 eigenvalues :i?e/ii < 0, 



\/j,(a)\ -> +oo 



r — q eigenvalues: /Ki(cr) — > "(0) 



^(a) = 



Also, nil denotes the number of negative eigenvalues of An. Since E is positive definite, also En is and 
hence by Lemma [37T]rjii is equal to the number of strictly negative eigenvalues of E^(u)A\\(u). Hence, by 
the third assumption in Hypothesis [H ri\\ does not depend on u. Finally, q denotes the dimension of the 
kernel of An(u), which as well does not depend on u. 

Before introducing Lemma 14.71 it is necessary to state a preliminary result, which is exploited in the 
proof of the lemma: 

Lemma 4.6. Let Hypothesis^ hold. Then there exists S > such that if u belongs to a neighborhood of uq 
of size S then the following holds. The polynomial 



where the matrices a, be M( r ^'?) x ( r ^«) are defined by (|4.21|) and b is invertible by Lemma |4~S1 The symbol 
I r - q denotes the identity matrix of dimension (r — q) . 

Lemma 4.7. Let Hypotheses^ and\3\hold. Then the dimension d of the stable space (|4.34[) of system (|4.32[) 
is d = k — 1 — nn — q- 

Proof. Since the state u is fixed, in the following for simplicity the matrices A[u, 0) and B{u) will be denoted 
by A and B respectively. To prove the lemma one has to show that the equation 




(4.36) 



has degree (r — q) with respect to \x. 

The proof follows immediately from Lemma 14.51 which guarantees that 





(4.37) 



admits exactly (As — 1 — n\\ — q) roots with negative real part. 
The proof is organized into six steps: 
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1. First, it is introduced a perturbation technique. More precisely, it is considered the problem 

det [A - fj,i(a)(B + crl)) =0. a > 0, (4.38) 

where a is a positive parameter which is allowed to go to 0. The symbol / denotes the TV-dimensional 
identity matrix. Just to fix the ideas and to consider the most general case, it will be supposed that 
the matrix A is singular: 

Ai < • • • < Afe_i < = A fc < A fc+ i <■■■ < X N 

Once the value of the parameter a is fixed, the solutions //i(<r) of (|4.38|) are the roots of a polynomial 
of degree N. Moreover, the matrix (73 + a I) is positive definite and hence Lemma l3~T1 guarantees that 
(k — 1) roots have strictly negative real part, (N — k) roots have strictly positive real part and one root 
is zero. 

If one let the parameter a vary, (|4.38[) defines N algebraic functions /ii(cr), . . . //jv(c): since the coeffi- 
cients of the polynomial (14.38[) are analytic with respect to a, it is known (see for example [38.;) that 
the functions fii are analytic everywhere in the complex plane except for some so called exceptional 
points that are certainly of finite number in every compact subset. Moreover, for every fixed function 
Hi only two behaviors are possible in an exceptional point a = a: 

• Hi(cr) is continuous in a. 

• limber \Hi{o)\ = +oo. 

In the case of the algebraic functions defined by (|4. 38(1 . the point a = is an exceptional point since 
in a = the degree of the polynomial drops from N to r — q: we will denote by (Jti(0), i = 1, . . . (r — q) 
the (r - q) roots of (pj~37) . 

Thanks to the known results recalled before, however, one deduces that there are 

• r — q functions (J,i(cr) such that 

lim Hi{o) = fii(0) i = l,...r. 

a— »0 

• N — r + q functions fii (a) such that 

lim |//j(cr)| = +oo i = r + 1, . . . N. 

(7 — >Q 

The situation is summarized in Figure O 

Moreover, restricting to a small enough neighborhood, it is not restrictive to suppose that a = is the 
only singular point. 

2. As a second step we prove that to obtain the lemma it is sufficient to show that the number of functions 
fJ,i(cr) continuous in a = and such that Refii(a) < when a S K. + is equal to k — 1 — n\i — q. Since 
the number of roots of (|4.38[) with negative is equal to k — 1, this is equivalent to show that there are 
exactly rin + q functions /Ltj(<r) such that Re^i(a) < and |/ii(<r)| T +°° when a — > + . 

The only case that has to be excluded is the possibility that in the limit a function /ii(cr) with Refii(a) < 
has zero real part. 

First, one observes that equation (|4. 37[) does not admit purely imaginary roots: by contradiction, 
assume that there exists /j £ M, /i ^ and v = v r + ivi, v r , Vi G K N such that 

A(v r + iv^ = i/iB(v r + iiii). (4.39) 

Just to fix the ideas, one can assume fi > 0: as in the proof of Lemma l3~Tl it turns out that 

> -n(Bvi, Vi) = (Av r , v t ) = (Avi, v r ) = n{Bv r: v r ) > (4.40) 
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because B is positive semidefinite. Hence the only possibility is that Vi and v r both belong to kerB, 
but this contradicts (|4.39[) and hence there are no purely imaginary roots of (|4.37|) . 

Hence one is left to exclude the possibility that a function fii(a) continuous in a = and such that 
Re/j,i_(a) < converges to 0. From the strict hyperbolicity of the matrix E~ x A (Hypothesis [3]) and 
from Lemma 13.11 it follows that is an eigenvalue of A with multiplicity one and hence fik = is a 
root of (I4.37P with multiplicity one. Moreover, again from Lemma l3~Tl one deduces that the algebraic 
function //fe(cr) is identically equal to 0. Since the multiplicity of as a root of (|4.37[) is one, it cannot 
happen that other functions fii (a) , k converge to zero when a — > + . 

3. As a third step, we perform a change of variable. More precisely, the number of solutions Hi(cr) of 
(I4.38[) such that Re[ii(cr) < and |/ij(c)| j +oo when a — > + is equal to the number of roots x^(<r) of 



det \B + a I - Xi(a)A) = i=l,...N-r + q (4.41) 

such that Rexiia) < and Xi(a) — > when a — > + . This is clear if one defines Xi(a) := l/fj,i(a). 
4. In the fourth part of the proof, we study the inverse problem to (|4.41|) . i.e. the eigenvalue problem 

det ( B - xA + a-j (x)l) = j = l,... N. (4.42) 



The meaning of the notation is that in this case one fixes x and then finds a. 

More precisely, it will be studied the behavior of the N — r eigenvalues — (Tj(x) such that <Tj(x) — > 
when x — > 0. From the analysis in [T^] it follows that 

F(x) = -A n x - A% 1 b~ x A 2 ix 2 + o{x 2 ) x^O 

where F{x) denotes the projection of B — xA on the generalized subspaces converging to the kernel of 
B. The blocks An and A21 are defined by (|4.35p . Let Fq(x) be the projection of F(x) on the kernel 
of An and F±(x) be the projection on the subspace orthogonal to keiAn. 

Thus again the analysis in [H] guarantees that 

F (x) = -(Aiifb^A^iX 2 + o(x 2 ) x^O 

and 

Fj_(x) = -Aux + o(x) x^O. 
In the previous expression it has been used the notations introduced in (|4.9[) : 



Au(u, a) A 2 i{u(i) T 
A 2 i{ua) A 22 (u,u x ,(x) 



( (A^f^a) 

An(u,a) (Ai{f(u, a) 

\ A 2l( u ) A ¥l( u i CT ) ^22 (u, cr, U x ) 



From Kawashima condition it follows that the matrix {A 2X ) T b 1 A 2X is positive definite: one can proceed 
in the same way as in the proof of Lemma (|4.2p . 

Hence projecting on suitable subspaces (we refer again to [12] for the precise computations) it follows 
that the following expansions hold: 

(Tj(x) — fijX 2 + o(x) x — > 0, j = 1, . . . q (4.43) 

and 

a j (x) = XjX + o(x) x -> j = q + 1, . . . N - r. (4.44) 

In the previous expression, \ij > is an eigenvalue of {.A\i) T b^ 1 A\ x and Xj ^ is an eigenvalue of 
Ail- 
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5. As a fifth step, we analyze equation (14. 44j) : for every fixed j = q+ 1, . . . N — r, the expansion guarantees 
that the function <Jj(x) is invertiblc in a neighborhood of x = 0. The inverse function Xj(o~) satisfies 
(I4.41[) and because of (|4.44[) the condition that Rexj (a) < when a 6 M + is satisfied if and only if Xj 
is an eigenvalue of An with negative real part. This implies that there are exactly riu functions Xj(a) 
satisfying (|4.44p and such that Rexj (a) < when a G K + . 

Moreover, there are N — r — q functions Xj (a) satisfying (|4.41[1 and (|4.44|) . 

6. Since there are exactly 2q functions Xj(a) satisfying (|4.41[) and (|4.43[) . one deduces that there are q 
functions with strictly positive real part and q functions with strictly negative real part. The possibility 
of a root with zero real part is excluded on the basis of the same considerations as in step 2. 

Hence from steps 5 and 6 one obtains that there are exactly tin + q solutions Xi(a) of (|4.41[) such 
that Rexi(a) < and Xj_(a) — » when a — » + . Thanks to the considerations in steps 1, 2 and 3 this 
concludes the proof of the lemma. 

□ 

To introduce Lemma l4.8[ it is useful to introduce some further notations: let 

V s (u) := span<e nil _ g+1) . . . 6 fc _i) (4.45) 
be the space generated by the vectors that satisfy 

(A{u, 0) - /ii(u, 0)£(t2)) 0i(tt, 0) = [n(u, 0) < 0. (4.46) 

Let 

V u {u) :=span(H fc+1 ,...H JV ) (4.47) 
the subspace generated by the Si satisfying 

(A(u, 0) - X t E(u)^j E t = A, > 0. (4.48) 

When A(u, 0) is invertible the dimension of V u is equal to (JV — n), where n is the number of negative 
eigenvalues of A defined as in (|3.34p . On the other side, when A(u, 0) is singular the dimension of V u is 
equal to (N — k), where (k — 1) is the number of strictly negative eigenvalues of A as in (|3.40p . In this case, 
the following subspace is non trivial: 

V c {u) := span(e fc ), where A(u, 0) 6 fc = 0. (4.49) 

The proof of the following result is analogous to that of Lemma 7.1 in [6/. 

Lemma 4.8. Let Hypotheses^ and\^ hold. Then 

v s nv u = {o} v s nv c = {o} v c nv u = {o}, 

where V s , V c and V u are defined by (|4.45j) . (|4.49|) and (|4.47p respectively. 



4.2.3 The hyperbolic limit in the non characteristic case 

In this section we will provide a characterization of the limit of the parabolic approximation (|4.ip when the 
boundary non characteristic, i.e. when none of the eigenvalues of E~ 1 (u)A(u, u x ) can attain the value 
(Hypothesis [8]) . As in Section 15. 2. 2\ n will denote the number of eigenvalues of E~ 1 (u)A(u, u x ) with strictly 
negative negative real part and N — p the number of eigenvalues with strictly positive real part. Also, we 
recall that n\i is the number of strictly negative eigenvalues of Au(u), while the eigenvalue has multiplicity 
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To give a characterization of the limit of (|4.1[) we will follow the same steps Section 13.2.21 Thus, in 
the exposition we will focus only on the points at which the singularity of the viscosity matrix B plays an 
important role. 

The characterization of the limit works as follows. We will construct a map 4>(uq, s nil + q . . . sn) which 
describes, as (s nil+q . . . sn) vary, states that can be connected to uq. We will compose <j> with the function 
/?, which is used to assign the boundary condition in (|4.ip . We will the show that the composite map is 
locally invertible. This means that, given uo and g such that \/3(uq — g\ is sufficiently small, the values of 
(snu+q ■ ■ ■ sn) are uniquely determined by the equation 

g = 4>(u , si . . . s N ). 

Once (s nil +q . . . sn) are known the limit of (|3.1|) is completely characterized. The construction of the map 
cj> is divided in some steps: 

1. Waves with positive speed 

Consider the Cauchy datum uq, fix (N — n) parameters (s„+i . . . sn) and consider the value 

u = Tlc..T s N N u Q . 

The curves TJ^ . . -T^ are, as in Section 14.2.11 the curves of admissible states introduced in [7]. The 
state uo is then connected to u by a sequence of rarefaction and travelling waves with positive speed. 

2. Boundary layers 

We have now to characterize the set of values u such that the following problem admits a solution: 

A(U,U X )U X = B(U)U XX 

17(0) = u (4.50) 
lim^+oo U(x) = u. 

Because of Lemma 14.41 one has to study 

f U x = (w(U, z, 0), w(U, z, 0), z(U, 5, 0), z) 5 
\ z x = f(U,z,0) 

Consider the equilibrium point (u, 0), linearize at that point and denote by V s the stable space, i.e. 
the eigenspace associated to the eigenvalues with strictly negative real part. Thanks to Lemma 14.71 
the dimension of V s is equal to n — tin — q. 

n n 
i=l ^ l ^ U ' j=l 

where ni(u) ... fi n (u) are the eigenvalues of 6 _1 (u, 0, 0)o((m), 0, 0) with negative real part and £i(u) . . . £n 
are the corresponding eigenvectors. 

Denote by A4 S the stable manifold, which is parameterized by V s . Also, denote by 4> s a parameteriza- 
tion of M s : 



Let 7r u be the projection 



cj) s : V s -> R N . 



l N x R r ~ q R N 



(u, z) u 

One can then proceed as in the proof of Section 13.2.21 and conclude that system (|4.50[) admits a 
solution if u € n u (0 s (s nil + 9 +i . . . s n )J for some s„ 11+g+ i . . . s„ . Also, thanks to (|4.22[) the columns of 
the jacobian of n u o tfi s computed at s„ 11+g+ i = . . . s n = are S„ 11+g+ i . . . H n . 

Note that the map n u o (f> s actually depends also on the point u and it does in a Lipschitz continuos 
way: 

\n u o (f> s (Ul, S nil+ g + l . . . s n ) - n u o 4> s (u 2 , s„ 11+?+ i . . . s n )\ < L\ui - u 2 \. 
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3. Conclusion 

Define the map (j> as follows: 

<j>(uo, s nil+q+1 ...s n ) = tt u o (j) s (T s N N Z k k o . ..T^uo, s nil+q+1 . ..a n J (4.52) 

From the previous steps it follows that <f> is Lipschitz continuos and that it is differentiable at s„ 11+g+ i = 
. . . sn = 0. Also, the columns of the jacobian are S nil+g+ i(ito) . . . S„(uo), 0„+i(uo) . . . Ojv, where 

(a(u , 0) - Ai(uoB(uo))ei(tio) = 

for Aj(uo) > and 

\A(u , 0) - (M(ii E(vQ)jEi(uo) = ^(uoJSjf^to) 
for jUi(wo) with strictly negative real part. 

In the case of an invertible viscosity matrix (Section l3.2.2| the definition of the map <j) is the final step 
in the construction. Here, instead, one has to take into account the function fi, which is used to assign 
the boundary condition and is defined in Section [2.2.11 Consider 

fi o cj)(u , s„ 11+9+ i . . . S N ). 

Thanks to the regularity of 6 and to the previous remarks, Bo0 is Lipschitz continous and differentiable 
at s nil+q+ i = 0. . . sat = 0. Denote by 

V(u ) = span(Si(M ) . . . H„(u ), <d n+ i(u ) ■ ■ ■ Ojv(uo)- 

Lemma 1431 which is introduced in Section [431 ensures that for every V £ V(uq) 

DQ(u Q )V = => V = Q 

Thus, the jacobian of fio <fi at s nil + g +i = . . . sjy = is an invertible matrix. Thanks to the extension 
of the implicit function theorem discussed in [21] (page 253) one can conclude that the map Bo0(u o , •) 
is invertible in a neighbourhood of (snu+g+i . . . ajv) = (0 . . . 0). In particular, if one takes Uf, as in 
(13. ip and assumes that |6(uo) — g\ is sufficiently small, then the values of Snu+g+i ■ ■ • sjv are uniquely 
determined by the equation 

g = fi o (f>(u , s nil+q+ i . . . s„) (4.53) 

Once the values of s nil + 9 +i . . . sjv are known, then the limit u(t, x) can be reconstructed. In particular, 
the trace of u on the axis x — is given by 

u:=TZ\\o...T s N n u . (4.54) 

Also, the self similar function u can be obtained gluing together pieces like (|3.33p . 

Here is a summary of the results obtained in this section: 

Theorem 4.1. Let Hypotheses^ El ? and El hold. Then there exists 5 > small enough such that 
the following holds. If \R(uq) — g\ < 5, then the limit of the parabolic approximation (|3.1|) satisfies 

g = 60(wo, s nil+q+ i . . . s N ) 

for a suitable vector (s r in+<j+i • • • %)■ The map 4> is defined by (|4.ip . Given uq and g, one can invert fi o <p 
and determine uniquely {s nil + q +i . ..Sjv). Once (s T in+g+i •••Sjv) are known the value u(t, x) assumed by 
the limit function is determined a.e. (t, x). In particular, the trace u of the hyperbolic limit in the axis x = 
is given by ()4.54jl . 
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4.2.4 The hyperbolic limit in the boundary characteristic case 



This section deals with the limit of the parabolic approximation (|4.ip in the boundary characteristic case, i.e. 
when one of the eigenvalues of E~ 1 A can attain the value 0. More precisely, we will assume Hypothesis HI 
which is introduced at the beginning of Section 13.2.31 As in Section l3.2.3[ — 1 is the number of eigenvalues 
of E~ 1 (uo)A(iio, 0) that are less or equal then — c, where c is the separation speed introduced in Hypothesis 
[3J Also, Tin is the number of strictly negative eigenvalues of Ah(uq), while q is the dimension of the kernel 
ofAn(uo). 

To give a characterization of the limit of (|4.ip one can follow the same steps Section l3.2.3l Also, in Section 
I4.2.3l wc explain how to tackle the difficulties due to the fact that the viscosity matrix is not invertible. Thus, 
in the following we give only a quick overview of the key points in the characterization. 

The characterization works as follows. We construct a map 4>(uq, s nil + q . . . sjv) such that the following 
holds. As (s„ 11+g . . . sn) vary, <f>(uo, s nil + q . . . sjv) describes states that can be connected to uq. We compose 
<j) with the function /3, which is used to assign the boundary condition in (|4.ip . We then show that the 
composite map is locally invertible. This means that, given u and g such that \(3(uo—g\ is sufficiently small, 
then the values of (s nil + g . . . sjv) are uniquely determined by the equation 



Once (s nil+q . . . sn) are known the limit of (|3.1|) is completely characterized and it is obtained gluing together 
pieces like 13.331 The construction of the map 4> is divided in some steps: 

1. Waves with positive speed 

Given the Cauchy datum uq, we fix (N — k) parameters (s^+i • • ■ % and consider the value 



The curves T„ fc+1 . . . TiY are, as in Section 14.2.11 the curves of admissible states introduced in 7 . The 
state uq is then connected to u by a sequence of rarefaction and travelling waves with positive speed. 

2. Analysis of the center stable manifold 

Consider the equation satisfied by travelling waves: 



The point (v,k, 0, A&) is an equilibrium. One can than define an invariant center stable M cs manifold 
with the same properties listed in Section 13.2.31 Thanks to Lemma I4.7[ the dimension of every center 
stable manifold is k — n±i — q. Also, one can proceed again as in Section 13.2.31 and find the equation 
satisfied by the solutions of 14.2.11 laving on A4 CS . Moreover, every solution laying on M cs can be 
decomposed in a purely center component, a purely stable component and a component of perturbation, 
in the same way described in Section ^. 2. 31 Eventually, one is able to define a map F(v,k, s rill + g , . . . sk ) 
which is Lischitz continous with respect to both Uk and s nil + q , . . . sk- Also, it is differentiable at 
s nil + q = 0, . . .Sk = and the columns of the jacobian matrix are the vectors S„ 11+g (ufc) . . . Efc(ufc). 
For every i = rin + q ... K it holds 



g = fi o (j>(u , si . . . s N ). 



rrik+l rriN - 

u k=T Sk 7 +1 °---T SN u . 



B{U)U" = (A-aE)U'. 



Thanks to Lemma 14.41 this is equivalent to system 




(4.55) 



A(u k , 0) - HiE(u k Sj = S(u fe )Hj 



with the real part of \Xi less or equal to zero. 
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3. Conclusion 



Define the map <f> as 

cf>(uo, s nil+q+1 ...s n ) = tt u o (f> s (T^ N ~_ k k o . ..T^Uq, s nil+q+1 . ..s n J (4.56) 

From the previous steps it follows that <f> is Lipschitz continuos and that it is differentiable at s„ 11+g+ i = 
... sn = 0. Also, the columns of the jacobian are Si(uo) ■ • ■ Sfe(uo), Qk+i(fio) . . . 0jv, where 

(a(u , 0) - Ai(u S(u ))e i (i2o) = o 
with Xi > for every i. The vectors H are as before. 

To take into account the function /3, which is used to assign the boundary condition and is defined in 
Section 12.2. II . one considers 

fi o cj)(u , s nil+g+ i . . . s N ). 

Thanks to the regularity of 6 and to the previous remarks, Bo0 is Lipschitz continous and differentiable 
at s nil+q+ i = 0. . . sat = 0. Denote by 

V(u ) = span(Si(w ) ■ ■ •Sfc('So), fc+ i(w o ) ■ ■ • 6jv). 

Lemma l4~9l which is introduced in Section PO) ensures that for every V £ V(uo) 

DQ(u )V = V = 

Thus, the jacobian of Bo <f> at s nil+q+ \ = . . . sjy = is an invertible matrix. Thanks to the extension 
of the implicit function theorem discussed in [24 (page 253) one can conclude that the map Bo (j>(uo, •) 
is invertible in a neighbourhood of (s„ 11+9+ i . . . sn) = (0 . . . 0). In particular, if one takes as in 
([3.1)1 and assumes that \Q(uq) — g\ is sufficiently small, then the values of s nil + 9 + 1 ... sat are uniquely 
determined by the equation 

g = fi o (f>(u , s nil+q+ i . . . s n ) (4-57) 

Once the values of s nil + g +i . . . sn are known, then the limit u(t, x) can be reconstructed gluing together 
pieces like 13.331 In particular, the value of the trace u on the axis x = can be determined in the same 
way described in Section ["3. 2. 31 

Here is a summary of the results obtained in this section: 

Theorem 4.2. Let Hypotheses^ [3 > and\Q hold. Then there exists 5 > small enough such that 
the following holds. If |fi(uo) — g\ < 5, then the limit of the parabolic approximation (|3.1|) satisfies 

g = 6 o (j)(u , s nil+q+ i . . . sjv) 

for a suitable vector (s„ 11+g+ i . . . sjv). The map (f> is defined by ()4.56[) . Given uq and g, one can invert fi o <p 
and determine uniquely (s„ 11+g+ i . . . sn)- Once (s„ 11+9+ i . . . Sjy) are known, then the value u(t, x) assumed 
by the limit function is determined a.e. (t, x). 



4.3 A transversality lemma 

In the first part of this section we state and prove Lemma l4.9l It is a technical result and it guarantees that 
the map 6 o <f> that appears in both Theorems 14.11 and 14.21 is indeed locally invertible. In the second part of 
the section we introduce a new definition for the map fi which is used to assign the boundary condition in 

E{u)u t + A(u, u x )u x = B(u)u xx , s 

8{u{t, 0)) = g u(0, x) = Qq. K ' 
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This definition is an extension of Dcfiniton 12.11 and it guarantees both the local invertibily of the map 6 o <j) 
and the well posedness of the initial boundary value problem (|4.58[) . 

Before stating Lemma l4~9l we recall some notations. We denote by An, En the blocks of A, E defined 
by (14.35|) . Let Q(u) G R N ~ r be an eigenvector of fif/fwjiufu) associated to an eigenvalue with non positive 
real part. Also, let Zi(u) £ M. N be defined by 



d 
o 

and 

Z(u) := span(Zi(w), . . . , Z nil+g (u)). 
We define the complementary subspace 



(4.59) 



W(u) := span{| 



ei r • • ' ' I e r r V )'"''{ 



where e*i £ W are the vectors of a basis in W and Wj G R w_r are the eigenvectors of E^An associated to 
eigenvalue with strictly positive real part. By Definition 12.11 the function fi(u) is the component u\ of the 
vector u in the decomposition 

u = ux+u 2 , u% e W(u),u 2 e V(u). (4.60) 
Also, we define the space V{u) as 

V(u) := span(e„ 11+(Z+ i(u), . . . , Q k -i(u), 6 fe (u), E k+1 (u), . . . E N (u)), (4.61) 

when the matrix A(u, 0) can be singular (boundary characteristic case). The vectors Qi satisfy 

A(u, 0)Qi(u) = inBBt Re(fn) < 0, 

while the vectors 3j satisfy 



(A(u, 0) - Xi(u)E{u)^ Ei = Ai(u) > 



In the case of a non characteristic boundary (i.e when A(u, 0) is always invertible) the subspace V(u) is 
defined as 

V(u) := span(8„ 11+g+ i('u, 0),...,6„(u, 0), E n+ i, . . .E N (u, 0)}, 

With n we denote the number of negative eigenvalues of A, according to (|3.34p . According to the results 
in Sections 14.2.31 and 14.2.41 V(uo) is the space generated by the columns of the jacobian of the map <fi at 

(Snii+q+l ■ ■ ■ Sat) = 0. 

The following lemma guarantees that, even if W(u) and V(u) do not coincide in general, nevertheless 
they are transversal to the same subspace Z(u). This result is already known, but we repeat the proof for 
completeness. 

Lemma 4.9. The following holds: 

V(uo) © Z(uq) — R N W{u )(BZ(u )=R N 

Proof. The second part of the statement is trivial. 

In the proof of the first part for simplicity we consider only the boundary characteristic case (A(u, 0) can 
be singular), being the other case absolutely analogous. Also, we write A instead of A((uq, ). The proof is 
organized into four steps: 

1. First, we prove that if V is a non zero vector in span(Sfc + i, . . . , En), then 

(V, AV) > 0. (4.62) 
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The steps are analogous to those in the proof of Lemma 7.1 in [BJ, but for completeness we repeat 
them. 

If V £ span(Sfc + i, . . . , Sjv), then V belongs to the unstable manifold of the linear system 

u x = Au 

and hence there exists a solution of 

u x — Au 

u(0) = V lim^^^oo u(x) = 0. 
Thanks to the symmetry of A, such a solution satisfies 

— {u, Au) = 2{u x , Au) = 2\Au\ 2 > 
ax 

and hence to conclude it is enough to observe that 

lim (u(x), Au{x)) = 0. 

x^ — oo 

2. if 

nu+q 

Z = y ' XiZi, 

then 

(Z, AZ) < 0. (4.63) 

Indeed, the matrices An, En are symmetric and En is positive definite, hence they admit an orthog- 
onal basis of eigenvectors Ci • • ■ Ov-r such that (Q, AnQ) = rfcSij (&, En^). In particular, 

nu+q 



(z, az) = J2 mxi{Q> E u&<o. 



! = 1 



3. if 



k 



then 



(AO, 6) < 0. (4.64) 
Indeed, by Lemma 14.41 one deduces that the following system admits a solution: 

Bu x = Au 

u(0) = lim^+oo u(x) = 0. 

Hence, by considerations analogous to those performed in the first step, one concludes that (|4.64p holds 
true. 

4. it holds 

span(Zi, . . . , Z nil+q ,Q nil+q+ i, . . . ,0 fc ) n span(3 fc+1 , . . .,Hjv> = {0}. (4.65) 
To prove it, it is enough to show that if 



Z e span(Zi,...,Z nil+g >, 6 £ span(9„ 11+9+ i, . . . , 6 fc ), (Z + 9)^0, (4.66) 

then 



(z + Q^j ^span(3 k+1 ,...,Sjv). 
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If 

N 

6 = 2J ^Qii 

t=nn+g+l 

then 

AT JV 

<AZ, 6) = {Z, AO) == J2 Xi ^ Z ' AQ ^ = J2 IH*i(Z, BQ t ) = 

i=n 11 +q+l i=mi+g+l 

and 

(Z + Q,AZ + A6) = (Z, AZ) + 2(Z, AQ) + (6, AO), 
from (T05|) , (TO^|) and dUMt it follows that holds true. 

5. to conclude, it is enough to show that for every i = 1, . . . , (mi + (?) 

Z< £ span(e„ 11+g+ i(u, 0),...,6 fe (u, 0)>. (4.67) 

Suppose by contradiction that 

k 

.7=7111+9+1 

for suitable numbers c„ 11+9+ i, . . . c&. Because of Lemma 14. 5 [ for every j, Qj has the structure described 
by (|4.23p . where £j is an eigenvector of b a. The matrices b and a have dimension (r — g) and are 
defined by (|4.2ip . Considering the last (r — q) lines of the equality (|4.68[) one obtains 

= E c i&> 

3 

which implies Cj — for every j, because the £ 3 - are all independent. Hence (|4.68p cannot hold and 
(|4.67p is proved. 

□ 

From Lemma 1431 we deduce that the map Qocf) which appears in Theorems 14.11 and 14 . 21 is locally invertible. 
We proceed as follows. 

By construction the kernel of the jacobian DQ is Z(uq). Thus, 

DQ(u )V = 6 => VeZ(u ). 

The columns of the matrix 

DQ(u )D<p(u ) 

are D fi(i2o)6 nil + 9 +i . . . .Dfi(£to)0fc_i, DQ^Uq)^ . . .D6(uo)Sjv. To prove that the columns are all indipen- 
dent it is enough to show that 

fc-l N 

E XiDQ(u )Qi +'^2x i DQ(u )E i = => x nil+q+1 = . . . x N = 0. (4.69) 

nii+g+l £= A: 

Since O nil + g +i . . . &k-i, 3fe • • • Sjv are all indipendent, to prove (14.69[) is is enough to show that if V € V(uq) 
and 

DQ(u )V = 

then V = 0. This a consequence of Lemma 14. 9[ which states that 

V(uo) © Z{u ) = R N 



We now introduce a generalization of Definition 12.11 
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Definition 4.1. The function 

8 : R N -> l"-""- 9 

is smooth and satifies 

kernel(D8(u )) © W(u ) = R N (4.70) 

and 

kernel(r>fi(w )) © V(u ) = K W (4.71) 

Thanks to (|4.70p . the initial boundary value problem 

E(u)u t + A(u, u x )u x = B(u)u xx 
R(u(t, 0)) = g u(0, x) = uq 

is well posed. This is a consequence of the same considerations as in Section [2. 2. II 
Also, because of (|4.71l) the matrix 



is invertible. Thus, the function fio (f> which appears in Theorems 14.11 and 14.21 is locally invertible. In other 
words, the analysis in Section 2] is still valid if we use Definition 14. II instead of Definition l2.il 
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